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Section  1 
SUMMARY 


The  purpose  of  the  High-Altitude  Aircraft  Wake  Dynamics  Study  has  been  to  investigate 
the  chemically  reacting  wake  of  an  aircraft  flying  at  subsonic  and  supersonic  velocity 
in  the  upper  troposphere  and  stratosphere.  In  this  study,  the  chemical  and  fluid 
mechanical  behaviors  of  important  emission  species  were  traced  from  the  time  the 
species  exited  the  engine  exhaust  nozzle  to  the  time  aircraft-induced  perturbations 
to  the  atmospheric  environment  were  no  longer  important.  The  primary  result  of  this 
study  will  be  a wake  model  that  will  predict  the  concentration  of  engine  emission 
species  at  various  points  in  the  wake.  The  final  output  of  the  wake  model  will  be  the 
size  and  shape  of  the  wake  and  the  species  concentrations  at  the  time  when  the  atmos- 
pheric environment  dominates  the  dispersion  and  chemical  reactions  in  the  wake. 

The  High-Altitude  Wake  Dynamics  Study  has  been  organized  into  the  following  four 
major  tasks: 

• Task  I — Problem  Definition 

• Task  n — Model  Development 

• Task  HI  — Supporting  Experiments 

• Task  IV  — Model  Application  and  Update 

This  report  contains  the  technical  analyses  and  studies  performed  under  Task  I, 
Problem  Definition,  and  the  conclusions  developed  from  these  analyses.  Data  ob- 
tained under  Task  in,  Supporting  Experiments,  are  also  presented,  as  appropriate, 
to  support  the  findings  of  the  analytical  studies. 

The  relationship  of  the  wake  dynamics  study  to  other  modeling  studies  within  the 
Climatic  Impact  Assessment  Program  (CIAP)  is  illustrated  in  Figure  1-1.  The  wake 
regime  extends  from  the  time  engine  exhaust  emissions  are  deposited  in  the  atmos- 
phere  up  to  an  order  of  10  sec  thereafter.  The  wake  model,  which  will  be  developed 
under  Task  n,  Model  Development,  will  describe  fluid  mechanical  and  chemical  be- 
havior of  engine  exhaust  emission  in  the  wake  regime.  It  will  also  provide  final  wake 
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size  and  shape,  and  species  concentration  input  data  to  the  diffusion  and  transport 
regime,  which  begins  at  the  end  of  the  wake  regime  and  extends  over  a time  interval 
on  the  order  of  days.  The  diffusion  and  transport  regime  is  currently  under  study  by 
other  CIAP  contractors.  Finally,  global  circulation  models  will  extend  analysis  of 
the  interaction  of  engine  exhaust  emissions  with  the  atmosphere  up  to  a period  on  the 
order  of  years  after  their  introduction  into  the  atmosphere. 

The  basic  objective  of  Task  I was  to  examine  the  interaction  of  high-altitude  aircraft 
engine  emissions  with  the  aircraft  wake  flow  field  and  to  determine  the  extent  to  which 
this  flow  regime  must  be  considered  in  specifying  the  types  and  amounts  of  effluents 
entering  the  global  circulation  as  a result  of  high-altitude  aircraft  operations.  The 
general  approach  to  Task  I was  to  examine  the  fluid  mechanical  and  chemical  behavior 
of  the  engine  exhaust  emission  to  identify  and  evaluate,  where  possible,  the  important 
mechanisms  by  which  the  engine  emissions  can  interact  with  the  atmosphere.  For 
those  mechanisms  that  were  deemed  to  be  important,  the  dependence  of  parameters 
that  control  the  mechanism  of  interest  on  aircraft  characteristics,  flight  conditions, 
and  atmospheric  environment  were  established.  In  this  manner,  different  flight  con- 
ditions and  differences  in  the  atmospheric  environment  between  the  troposphere  and 
stratosphere  are  taken  into  account. 

The  important  features  of  the  wake  fluid  mechanical  model  were  found  to  be  the  domi- 
nance of  the  early  wake  growth  (up  to  a few  minutes)  by  the  engine  exhaust  jet  flow 
field  and  the  rolled-up  vortex  pair  flow.  At  early  times,  the  wake  growth  is  primarily 
in  the  vertical  direction  because  of  mass  detrainment  from  the  sinking  vortex  pair. 

The  vortex  pair  breaks  up  because  of  pair  instability  at  a wake  age  of  1 to  2 min  but 
continues  to  sink  for  an  additional  5 to  10  min.  At  the  time  of  vortex  pair  breakup 
(wake  age  « 1 to  2 min)  the  wake  width  is  about  twice  the  aircraft  wing  span  and  the 
wake  height  or  thickness  is  about  four  times  the  wing  span.  At  a wake  age  between  5 
and  10  min,  depending  on  the  particular  aircraft,  the  wake  height  reaches  a maximum 
of  about  seven  times  the  wing  span.  The  fluid  mechanical  behavior  beyond  this  point 
in  time  is  not  perfectly  clear.  However,  preliminary  data  from  contrail  photographs 
taken  in  the  troposphere  indicate  that  after  reaching  the  maximum  height  of  about  six 
wing  spans,  the  wake  begins  to  collapse  — because  of  gravitational  forces  — until  the 
wake  height  is  about  1/3  to  1/2  the  maximum  height.  This  collapse  occurs  in  about 
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5 to  10  min  after  the  maximum  height  is  attained.  The  gravitational  collapse  of  the 
wake  necessitates  an  increase  in  wake  width  by  a factor  of  2 to  3 to  conserve  mass. 

Additional  mechanisms  for  dispersing  the  wake  at  these  late  times  are  vertical  wind 
shear,  aircraft  wake  turbulence,  and  buoyant  rise.  The  deformation  of  wake  by  vertical 
wind  shear  is  directly  proportional  to  the  wake  height  and  hence  reaches  a maximum 
and  then  decreases.  The  turbulence  deposited  into  the  wake  by  vortex  pair  breakup 
is  expected  to  decay  to  the  atmospheric  level  in  about  2 to  3 min  after  deposition.  The 
decay  mechanisms  are  wake  mass  dilution  by  low -turbulence  atmospheric  air  and 
inviscid  dissipation  that  transfers  turbulent  kinetic  energy  from  large-scale  to  small- 
scale  eddies.  The  contribution  of  wake  and  atmospheric  turbulence  to  the  dispersion 
of  the  wake  is  relatively  small  in  the  wake  dispersion  regime. 

The  increase  in  wake  width  due  to  mass  entrainment  as  the  wake  rises  by  buoyancy  is 
also  expected  to  be  small.  The  residual  temperature  excess  in  the  wake  is  not  suf- 
ficiently large  (of  the  order  of  a degree)  to  cause  the  wake  to  rise  and  entrain  a signifi- 
cant amount  of  additional  mass. 

The  above  description  of  fluid  mechanical  behavior  of  aircraft  wakes  is  based  on  ob- 
servation of  contrails  in  the  troposphere.  It  is  currently  believed  that,  in  the  strato- 
sphere, the  wake  would  behave  qualitatively  in  the  same  fashion,  but  the  quantitative 
behavior  will  be  different  because  of  differences  in  wind  shear,  atmospheric  stability, 
and  atmospheric  turbulence.  The  wake  model  under  development  will  predict  the 
quantitative  differences  in  wake  behavior  between  tropospheric  and  stratospheric  flight 
altitudes. 

The  influence  of  supersonic  Mach  numbers  on  aircraft  wake  dynamics  is  poorly  under- 
stood at  the  present  time.  The  effects  of  shock  waves  and  differences  in  lift  distribution 
between  subsonic  and  supersonic  wings  on  vortex  roll-up  processes  are  not  known. 
Supersonic  contrail  data  are  necessary  to  provide  a basis  for  assessing  the  influence 
of  Mach  number  on  aircraft  wakes.  A flight  program  to  obtain  wake  measurements 
behind  a supersonic  aircraft  is  currently  being  planned. 


1-4 


The  possible  effect  of  shock  heating  on  the  chemical  dynamics  is  somewhat  better 
understood.  The  maximum  air  temperatures  created  by  shock  heating  by  a Mach-3 
aircraft  are  in  the  neighborhood  of  600 °K.  Temperatures  of  this  magnitude  are 
created  by  normal  shocks  that  are  located  very  close  to  the  aircraft  surfaces.  A well 
designed  supersonic  aircraft  will  not  create  large  normal  shocks  because  of  the 
excessive  drag  associated  with  such  blunt  configurations.  Most  of  the  shock-heated 
air  is  probably  swallowed  by  the  boundary  layer  where  recovery  temperatures  are 
comparable  to  the  normal  shock  temperatures.  The  hot  boundary  layer  shed  by  the 
aircraft  is  expected  to  cool  rapidly  bj^  mixing  and  is  not  expected  to  have  a significant 
influence  on  the  chemical  dynamics. 

The  important  conclusions  from  the  wake  chemical  kinetics  analyses  are: 

• In  an  engine  operating  mode,  with  relatively  low  HC  emissions,  finite  reaction 
rate  effects  are  an  important  aspect  of  the  chemistry  in  the  engine  exhaust 
expansion  nozzle  and  exhaust  jet  regime.  There  is  substantial  oxidation  of 
CO  to  CO^  in  the  nozzle  and  jet,  and  the  concomitant  production  of  H atoms 
maintains  the  concentrations  H,  O,  and  OH  well  above  their  equilibrium 
values. 

• Hydrocarbon  oxidation  has  a profound  scavenging  effect  on  the  inorganic 
radicals.  Introduction  of  large  HC  quantities  reduces  the  H,  OH,  and 

O levels  by  many  orders  of  magnitude . 

• Significant  oxidation  of  NO  to  N02  cannot  be  accomplished  in  the  ex- 
haust regime.  In  the  absence  of  HC , the  N09/N0  ratio  is  kept  very 
small  through  reduction  of  NO^  by  H atoms.  When  HC  is  abundant, 
the  effect  is  to  freeze  NO  in  the  nozzle  near  its  initial  throat  value. 

Similarly,  oxidation  of  SO?  to  SOg  appears  not  to  be  significant. 

e OH  concentrations  can  be  high  in  the  jet,  but  under  such  conditions, 
the  N02  would  be  very  low  — thereby  precluding  formation  of  HNO^ 
from  the  reaction 

no£  + oh  — hno3 
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• Even  if  OH  and  N02  were  both  initially  plentiful  in  the  hot  exhaust,  significant 
conversion  to  nitric  acid  is  forestalled  through  rapid  consumption  of  OH  by 
the  competing  reaction 


OH  + OH  + M — HO  + M 
z z 

The  wake  chemical  model,  which  will  provide  inputs  to  the  diffusion  and  transport 
regime  chemistry,  can  be  either  simple  or  complex  depending  on  requirements  of 
the  diffusion  and  transport  regime  model.  This  is  an  area  that  requires  further  study 
and  coordination  with  other  CIAP  program  elements. 

The  results  of  the  present  studies  have  indicated  that  radiative  cooling  of  the  exhaust 
jet  and  wake  buoyancy  due  to  radiative  transport  are  not  important.  The  transport  of 
particulates  in  the  wake  regime  was  also  found  to  be  unimportant  because  of  the  small 
size  of  the  emitted  particulates  that  causes  the  particulate  motion  to  coincide  with 
the  gas  motion. 
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Section  2 

FLUID  MECHANICAL  MODEL 


2.  1 OBSERVED  FEATURES  OF  AIRCRAFT  EXHAUST  TRAILS 

A study  of  the  data  gathered  in  this  program  under  Task  ni  and  the  data  reported  in 
Ref.  2-1  yields  a consistent  picture  of  the  major  features  of  aircraft  exhaust  trails 
such  as  shown  in  Figures  2-1  and  2-2. 

The  near  trail,  labeled  "Jet"  in  these  figures,  lasts  for  times  on  the  order  of  10  sec. 

In  this  region  the  exhaust  trail  occupies  a fraction  of  the  overall  aircraft  wake  since 
it  is  confined  to  turbulent  jets  issuing  from  the  engines.  The  jet  regime  is  character- 
ized by  rapid  growth  of  the  exhaust  trail,  which  starts  on  a meter  scale  at  the  engines 
and  grows  more  than  an  order  of  magnitude  to  become  comparable  to  the  aircraft  size 
at  the  end  of  the  jet  (Point  A,  Figure  2-2),  at  which  point  the  exhaust  occupies  most  of 
the  aircraft  wake.  The  end  of  the  jet  is  marked  by  a sharp  decrease  in  horizontal  rate 
of  growth,  as  shown  in  Figure  2-2. 

The  intermediate  trail  is  labeled  "Vortex,  " because  in  this  region  the  trailing  vortex 
pair  associated  with  wing  lift  dominates  the  wake.  The  striking  feature  of  this  regime 
is  a horizontal  growth  that  is  practically  nil,  while  the  vertical  growth  proceeds  at  a 
substantial  rate,  as  shown  in  the  Figures  2-1  and  2-2.  The  vortex  regime  lasts  for 
times  on  the  order  of  100  sec,  corresponding  to  distances  on  the  order  of  20  km  behind 
the  aircraft.  The  end  of  the  vortex  regime  (Points  B and  C of  Figure  2-2)  is  marked  by 
a reversal  of  the  growth  rate  pattern;  i.  e. , the  trail  grows  more  rapidly  in  the  horizon- 
tal direction  than  it  does  in  the  vertical  direction. 

After  the  vortex  pair  has  broken  down,  the  exhaust  trail  starts  dispersing  in  the  horizon- 
tal direction  by  a variety  of  mechanisms  to  be  discussed  later  and  continues  to  grow  in  the 
vertical  direction  because  of  the  downward  momentum  of  the  vortex  pair.  The  long-term 
wake  growth  behavior  is  illustrated  in  Figure  2-3  in  which  preliminary  wake  growth  data 
from  contrail  photographs  of  the  NASA  Ames  CV-990  are  presented.  The  wake  height 
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reaches  a maximum  at  a wake  age  of  about  5 min  and  then  begins  to  decrease.  This 
decrease  is  thought  to  be  due  to  gravitational  collapse,  which  occurs  when  a well  mixed, 
uniform  density  volume  of  gas  is  located  in  an  atmosphere  that  has  a finite,  stable  lapse 
rate.  The  wake  height  reaches  a minimum  between  15  and  20  min  and  then  begins  to  in- 
crease slowly.  The  methods  employed  to  obtain  the  data  presented  in  Figure  2-3  are  dis- 
cussed in  Section  2.2.4. 

The  gravitational  collapse  of  the  wake  necessitates  an  increase  in  the  wake  width  to  con- 
serve mass.  This  increase  in  wake  width  by  gravitational  collapse  occurs  primarily 
between  5 and  15  min.  Additional  mechanisms  for  dispersing  the  wake  are  vertical  wind 
shear,  aircraft  wake  turbulence,  and  buoyant  rise.  The  deformation  of  the  wake  by  verti- 
cal wind  shear  probably  becomes  effective  right  after  the  vortex  breakdown  and  is  directly 
proportional  to  the  wake  height.  The  turbulence  deposited  into  the  wake  by  vortex  break- 
down is  initially  much  higher  than  the  atmospheric  level  but  is  expected  to  decay  to  the 
atmospheric  level  in  about  2 to  3 min  after  deposition.  The  decay  mechanisms  are  wake 
mass  dilution  by  low  turbulence  atmospheric  air  and  inviscid  dissipation,  which  trans- 
fers turbulent  kinetic  energy  from  large  scale  to  small  scale  eddies.  Wake  growth  by 
wake  turbulence  diffusion  is  therefore  expected  to  be  confined  to  a few  minutes  after 
vortex  breakdown.  Thereafter  atmospheric  turbulence  will  disperse  the  wake.  The  in- 
crease in  wake  width  due  to  mass  entrainment  as  the  wake  rises  by  buoyancy  is  expected 
to  be  small.  The  residual  temperature  excess  in  the  wake  are  not  sufficiently  large  (of 
the  order  of  a degree)  to  cause  the  wake  to  rise  and  entrain  a significant  amount  of  addi- 
tional mass. 

Microdensitometer  traces  of  the  jet  flow  field  immediately  behind  the  aircraft  are  shown 
in  Figure  2-4  where  the  first  scan  shows  individual  jets  from  the  four  engine  pods,  which 
later  merge  two  by  two.  Additional  microdensitometer  traces  of  the  wake  photographs 
at  1.  5 and  7.  0 min  (Figure  2-1)  were  made  to  obtain  a qualitative  description  of  the 
mass  distribution  across  the  wake.  The  traces  are  shown  in  Figures  2-5  through  2-8. 

In  Figure  2-5  the  vortex  pair  is  clearly  visible  in  the  lower  portion  of  the  wake  and  the 
microdensitometer  traces  show  a greater  concentration  of  light  scattering  particles  in 
the  vicinity  of  the  vortex  pair.  The  scattering  particles  are  probably  ice  crystals  and  if 
their  dimensions  are  in  the  submicron  range,  their  distribution  should  be  a good  indicator 
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Figure  2-1  Shetch  of  Typical  Exhaust  Trail 
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The  photographs  show  vertical  and  lateral 
exposures  of  the  contrail  generated  by  a 
B-52  aircraft  flying  at  41,900  ft  and 
450  knots  over  San  Francisco.  Longi- 
tudinal distances  have  been  foreshortened 
in  the  drawing  to  show  the  different 
regimes. 
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Figure  2-1  Shetch  of  Typical  Exhaust  Trail 
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Figure  2-2  Time  History  of  B-52  Exhaust  Trail  Shown  in  Fig. 
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Figure  2-3  Convair  990  Contrail  Dimensions,  Preliminary  Data 
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increasing  distance  behind  aircraft 
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Figure  2-5  Vertical  Microdensitometer  Scans  - B-52  at  42  kft,  Wake  Age  = 1.5  min 
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Figure  2-6  Horizontal  Microdensitometer  Scans  — B-52  at  42  kft,  Wake  Age  = 1.  5 min 
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Figure  2-7  Vertical  Microdensitometer  Scans  — B-52  at  42  kft.  Wake  Age  = 7.0  min 
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Figure  2-8  Horizontal  Microdensitometer  Scans  — B-52  at  42  kft,  Wake  Age  = 7 min 


of  the  exhaust  gas  distribution.  Figure  2-6  is  the  corresponding  horizontal  trace,  which 
clearly  shows  the  mass  concentration  due  to  the  vortex  pair. 

Figures  2-7  and  2-8  are  vertical  and  horizontal  scans  at  7.  0 min,  well  into  the  wake 
dispersion  regime  where  organized  aerodynamic  motion  has  ceased  to  exist.  The 
microdensitometer  traces  indicate  a slightly  rounded  distribution  with  a ragged  top. 

The  concentration  of  mass  at  the  bottom  of  the  wake  evident  in  Figure  2-5  is  no  longer 
discernible.  The  breakup  of  the  vortex  pair  evidently  redistributes  the  mass  uniformly 
throughout  the  wake  cross  section. 

The  procedure  for  obtaining  the  data  shown  in  Figures  2-2  and  2-3  will  be  described  in 
detail  in  future  reports  dealing  with  Task  III  of  this  program.  For  present  purposes 
it  is  sufficient  to  indicate  that  the  photography  consists  of  still  and  motion  pictures  taken 
simultaneously  from  two  locations  — one  directly  underneath  the  flight  path  and  one  from 
the  side  at  an  elevation  angle  of  about  10  deg.  The  pictures  are  taken  in  a time  sequence, 
tracking  a particular  feature  of  the  contrail.  The  apparent  dimension  of  the  contrail  is 
measured  with  an  optical  comparator,  and  the  corresponding  height  and  width  are  ob- 
tained by  solving  a set  of  simultaneous  equations  involving  the  azimuth  and  elevation  from 
both  cameras.  In  a few  selected  cases  a microdensitometer  has  been  used  to  obtain  ac- 
curate definition  of  the  contrail  edge.  Typical  traces  are  shown  in  Figure  2-4  through 
2-8.  The  drift  velocity  of  the  contrail  is  obtained  by  geometric  addition  of  the  longi- 
tudinal drift  and  lateral  displacement,  measured  from  the  motion  picture  film.  Atmos- 
pheric conditions  are  obtained  from  the  local  sounding. 

The  exhaust  trail  is  assumed  to  coincide  with  the  observed  condensation  trail.  This 
assumption  is  based  on  the  fact  that  only  persistent  contrails  are  used  to  generate  data. 
This  type  of  contrail  is  produced  at  low  temperatures  (below  -61°  C)  and  persists  for 
hours,  which  indicates  that  it  consists  of  ice  crystals  in  a supersaturated  environment. 
Under  these  circumstances  reevaporation  would  be  negligible,  and  the  condensation 
trail  is  a good  indicator  of  the  exhaust  trail. 
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2.2  THEORETICAL  MODEL  OF  EXHAUST  TRAILS 


A theoretical  model  for  the  major  features  of  aircraft  exhaust  trails  has  been  developed 
on  the  basis  of  data  discussed  in  the  previous  section  and  from  theoretical  consideration 
of  the  physical  factors  affecting  aircraft  wakes.  For  present  purposes  the  end  of  the 
wake  is  defined  as  a point  where  the  sum  of  aircraft-induced  effects  has  become  negli- 
gible, and  the  natural  atmospheric  environment  dominates.  For  modeling  purposes  the 
wake  is  subdivided  into  jet,  vortex,  and  dispersion  regimes  (as  illustrated  in  Table  2-1), 
corresponding  to  the  observed  behavior  discussed  previously.  A summary  of  the  main 
features  of  the  model  follows,  and  a more  detailed  description  of  each  regime  is  included 
in  subsequent  subsections. 

The  jet  regime  extends  from  the  engines  to  the  point  where  the  jets  occupy  all  the 
recirculating  region  associated  with  the  trailing  vortex  pair  (indicated  as  Point  A in 
Figure  2-2).  This  corresponds  to  a time  on  the  order  of  10  sec  and  a distance  on  the 
order  of  2 km.  In  this  region  the  exhaust  gases  mix  with  the  aircraft  wake  according 
to  the  classical  laws  of  turbulent  jets. 

In  the  wake  model,  the  jet  flow  field  is  related  to  aircraft  characteristics  such  as  engine 
type,  power  setting,  and  nozzle  geometry,  and  to  flight  characteristics  such  as  flight 
altitude  and  velocity.  Given  the  aircraft  and  flight  characteristics  for  any  aircraft,  the 
jet  flow  field  can  be  predicted.  The  atmospheric  environment  in  this  regime  plays  a 
secondary  role  in  determining  the  jet  flow  field,  with  the  exception  of  the  ambient  pres- 
sure. The  ambient  pressure,  however,  maybe  classified  as  a flight  characteristic  since 
it  is  determined  by  flight  altitude. 

The  vortex  regime  extends  from  the  end  of  the  jet  regime  to  the  point  where  the  vortex 
pair  breaks  up  because  of  hydrodynamic  instability,  as  described  by  Crow  (Ref.  2-2) 
and  others  (Ref.  2-3).  This  regime  lasts  for  times  on  the  order  of  100  sec  and  extends 
over  distances  on  the  order  of  20  km.  The  unusual  feature  of  this  regime  - namely, 
vertical  growth  with  little  or  no  horizontal  growth  — can  be  explained  by  the  fact  that 
the  wake  in  this  regime  is  dominated  by  the  trailing  vortex  pair.  The  exhaust  gases, 
which  originally  fill  the  recirculating  flow  associated  with  the  sinking  vortex  pair, 
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Table  2- 
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(a)  See  Fig.  2-1  for  physical  features. 


become  detrained  from  the  pair  and  are  deposited  in  a region  of  roughly  rectangular 
cross  section,  with  the  top  at  the  original  flight  level  and  the  bottom  at  the  vortex  pair. 
Thus,  the  vortex  pair  restrains  the  exhaust  from  horizontal  growth  and  promotes  verti- 
cal growth  instead. 

The  characteristics  of  the  vortex  pair,  such  as  circulation  strength,  pair  separation 
distance,  time  to  vortex  sheet  rollup,  and  sinking  rate  are  related  to  aircraft  and  flight 
characteristics.  Some  of  the  more  important  parameters  are  aircraft  weight,  flight 
velocity,  lift  distribution,  and  wing  span.  Given  the  aircraft  and  flight  characteristics, 
the  major  features  of  the  vortex  pair  can  be  predicted. 

In  the  vortex  regime,  atmospheric  environment  again  plays  a secondary  role.  The 
structure  of  the  vortex  pair  is  determined  primarily  by  the  characteristics  of  the  air- 
craft and  flight  parameters.  The  rotational  motion  and  the  sinking  of  the  vortex  pair 
are  part  of  the  aerodynamics  of  the  aircraft. 

The  wake  dispersion  regime  lasts  from  vortex  breakup  to  a time  on  the  order  of  1,  000 
sec,  during  which  the  exhaust  changes  from  an  aircraft-dominated  condition  to  a mode 
of  natural  atmospheric  dispersal.  This  regime  has  been  generally  associated  with  pure 
atmospheric  dispersion,  and  an  unexpected  finding  of  this  program  is  that  the  aircraft 
induced  effects  last  for  times  about  an  order  of  magnitude  larger  than  the  vortex  breakup 
time.  At  the  beginning  of  the  dispersion  regime  (Points  B and  C in  Figure  2-2)  the  vortex 
pair  breaks  up,  allowing  the  horizontal  growth  to  resume  while  the  vertical  growth  rate 
decreases  markedly.  Eventually,  the  horizontal  spread  overtakes  the  height  (Point  D)  and, 
at  very  long  times,  the  exhaust  trail  assumes  the  flattened  shape  characteristic  of  atmos- 
pheric dispersion  under  stable  conditions. 

The  general  approach  in  the  wake  dispersion  regime  is  to  relate  the  four  dispersion 
mechanisms  — i.e. , wind  shear,  gravitational  collapse,  wake  and  atmospheric  turbulence, 
and  buoyant  rise  — to  the  appropriate  aircraft  parameters  and  the  atmospheric  environ- 
ment. In  this  manner,  given  the  aircraft  and  flight  characteristics  and  the  atmospheric 
environment,  the  wake  growth  in  the  wake  dispersion  regime  for  aircraft  flights  in  the 
stratosphere  can  be  estimated. 
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2.2.1  Jet  Regime 


The  jet  regime  extends  from  the  engines  to  the  point  at  which  the  exhaust  occupies  the 
entire  trailing  vortex  pair.  In  this  region  the  exhaust  is  confined  to  turbulent  jets 
issuing  from  the  engines  and  submerged  in  the  aircraft  wake.  As  the  jets  become 
immersed  in  the  vortex  field  and  are  attracted  to  the  cores  because  of  jet  buoyancy, 
the  interference  to  their  development  seems  minor.  Their  rate  of  growth  is  still  that 
of  classical  jets  immersed  in  a uniform  coflowing  stream.  The  details  of  jet-vortex 
interaction  and  appearance  of  the  exhaust  trail  are  dependent  on  the  detailed  aircraft 
configuration,  such  as  number  and  placement  of  the  engines.  However,  by  the  time 
the  end  of  the  jet  regime  has  been  reached,  the  exhaust  occupies  the  whole  vortex  re- 
circulation region,  and  the  details  of  entrainment  into  the  aircraft  wake  become  un- 
important, as  examination  of  flight  data  will  verify.  In  view  of  the  foregoing  arguments, 
the  jet  regime  is  modeled  using  the  classical  laws  of  a turbulent  jet  in  a coflowing 
stream  and  flight  data.  No  details  whose  effect  would  not  show  up  in  the  data  are  in- 
cluded in  this  model.  The  justification  of  this  simple  model  lies  in  its  ability  to  pre- 
dict experimental  results. 

The  jet  regime  is  important  primarily  because  it  determines  the  environment  in  which 
the  exhaust  may  still  be  chemically  active  as  the  result  of  the  relatively  high  temper- 
atures. As  indicated  in  Section  3,  it  has  been  found  that  the  important  chemical  reactions 
occur  either  in  the  engine  or  in  its  immediate  vicinity,  which  limits  the  need  of  detailed 
knowledge  of  the  jet  to  a region  of  the  order  of  10  exit  diameters.  A special  study  of  this 
region  has  been  made  and  is  described  in  the  following  section.  Although  the  rest  of  the 
jet  regime  (about  2 km  in  length)  is  not  of  crucial  importance  for  pollution  purposes,  it 
also  has  been  addressed  to  provide  a complete  description  of  the  exhaust  trail. 

The  Near  Field  — The  near  field  of  the  jet  regime  is  of  special  interest  because  it  is  the 
hottest  section  and  closest  to  the  engine,  and  therefore  is  a region  with  potential  for  sus- 
taining exhaust  reactions.  A detailed  study  of  the  near  field  must  take  into  account  the 
differences  in  engine  configuration,  particularly  with  respect  to  whether  or  not  the  engine 
has  a secondary  nozzle  . Most  modern  jet  engines  are  equipped  with  a variable-area 
primary  nozzle  that  sustains  sonic  flow  at  the  exit.  Thus,  under  most  operation  condi- 
tions, the  exhaust  is  underexpanded  at  the  exit  plane  and  undergoes  further  expansion  in 
the  jet.  On  the  other  hand,  some  high-performance  engines,  such  as  the  GE4  built  for 
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the  Boeing  SST,  are  equipped  with  primary  and  secondary  nozzles  able  to  expand  the 
exhaust  and  deliver  it  to  the  atmosphere  at  ambient  pressure.  This  feature  makes  an 
important  difference  because  it  gives  rise  to  two  different  cooling  mechanisms:  the 
underexpanded  jet  behind  a primary  nozzle  cools  rapidly  by  adiabatic  expansion  with 
very  little  mixing,  whereas  the  ideally  expanded  jet  behind  a secondary  nozzle  cools 
by  turbulent  mixing  with  atmospheric  air.  Both  types  of  jet  have  been  calculated,  and 
the  results  used  as  inputs  for  the  chemical  kinetics  calculations  described  in  Section  3. 

For  the  case  of  primary-nozzle,  underexpanded  jets,  calculations  are  carried  out  using 
the  method  of  characteristics  (Ref.  2-4).  This  program  has  the  capability  of  accounting 
for  external  flow  and  internal  shock-wave  effects.  Required  by  the  computations  are  a 
description  of  the  nozzle  geometry,  the  thermodynamic  properties  of  the  exhaust  at  the 
nozzle  exit  plane,  and  the  ambient  conditions.  Calculations  for  the  turbulent  shear 
layer,  where  mixing  occurs  at  the  edge  of  the  inviscid  jet,  are  carried  out  using  the  re- 
sults of  Ref.  2-5,  which  is  simple  to  apply  and  sufficient  to  provide  order-of-magnitude 
results.  The  shear  layer  involves  a minor  fraction  of  the  exhaust  up  to  four  or  five  exit 
radii,  after  which  there  is  no  need  to  calculate  the  jet  accurately  because  the  temper- 
ature has  decayed  to  values  low  enough  to  freeze  the  important  chemical  reactions.  A 
sample  calculation  using  the  theories  of  Refs.  2-4  and  2-5,  for  an  aircraft  flying  Mach 
2.26  at  19.  8 km  with  the  J85-GE-13  engine  operating  at  maximum  afterburner  setting 
is  shown  in  Figures  2-9  and  2-10.  For  modeling  purposes  a few  of  these  calculations 
will  be  made  and  the  results  correlated  to  provide  a simple  estimate  of  near -jet 
properties. 

In  the  case  of  an  ideally  expanded  (constant  pressure)  jet  behind  a secondary  nozzle, 
the  mixing  layer  at  the  edge  of  the  jet  grows  until  it  reaches  the  centerline,  after  which 
the  jet  becomes  fully  developed.  In  the  region  between  this  point  and  the  exit  plane, 
there  is  a conical  region  of  constant  properties  where  mixing  has  not  yet  taken  place. 

The  theory  of  Abramovich  (Ref.  2-6)  yields  the  following  length  for  the  conical  core: 
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Figure  2-9  J85-GE-13  Jet  Plume  Properties 
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where  u is  the  external  stream  velocity,  u.  is  the  exit  velocity  and  r is  the 
ext  igl  o 

nozzle  exit  radius.  A sample  calculation  of  the  ideally  expanded  jet  is  shown  in  Figure 
2-11  for  a GE-4  engine  at  maximum  afterburner  setting.  In  the  figure,  the  mixing  region 
width 
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where  X is  the  distance  downstream,  the  radius  of  the  core 
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and  the  velocity  profile 


are  similarly  given  by  the  theory  of  Ref.  2-6. 

The  Far  Field  — The  turbulent  jets  containing  the  exhaust  grow  in  the  developing  vortex 
field  generated  by  the  wing.  The  detailed  interaction  between  jets  and  aircraft  wake  is 
a complicated  process,  involving  the  buoyance  of  the  jets  and  the  roll-up  process  of  the 
vortex  sheet  behind  the  particular  type  of  wing  of  interest.  For  present  purposes,  these 
details  are  secondary;  the  main  objective  is  to  estimate  the  overall  rate  of  jet  growth 
and  location  of  the  point  at  which  the  exhaust  fills  the  vortex  pair.  To  this  end  the  calcu- 
lations of  Ref.  2-7  have  been  used  to  estimate  the  trajectory  of  various  points,  originally 
on  the  wing  trailing  edge,  as  they  proceed  to  roll  up  around  the  vortex  core.  The  con- 
clusion from  this  study  is  that,  for  wing-mounted  engines,  the  jets  turn  about  one  revo- 
lution around  the  vortex  core;  at  this  point  they  can  be  considered  as  fully  entrained  in 
the  vortex.  This  situation  is  confirmed  by  field  observations  and  the  study  of  photo- 
graphy of  smoking  or  contrail-producing  engines.  The  time  required  to  turn  one  revolu- 
tion around  the  core  is  as  follows: 
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Figure  2-11  Mixing  Layer  in  Jet  of  GE-4  Engine  Operating  at  Maximum  Power  at  Mach  Number  2.  7 and 
20- km  Altitude  (Ordinate  Scale  Increased  to  Facilitate  Reading  Graph) 
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where 

a 

= distance  from  the  engine  to  the  centerline  of  the  vortex 

b 

= aircraft  span 

V 

= aircraft  speed 

W 

= aircraft  weight 

P 

= ambient  density 

This  time  period  is  sensitive  to  aircraft  configuration  because  of  the  square  dependency 
on  H , but  a rough  average  is  i /b  « 0.  2,  which  yields 
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The  growth  of  the  jet  up  to  time  t^  has  been  calculated,  using  several  turbulent  jet 
theories  and,  in  particular,  that  of  Ref.  2-6  which  gives  the  jet  radius  in  the  implicit 
form 


The  parameter  p is  dependent  upon  the  ratio  of  external  velocity  to  jet  velocity  as 
given  by 


u 

PYt 


Calculations  using  these  formulas  and  photographic  data  show  that  at  the  time  t^  the 
jets  have  grown  to  approximately  the  size  of  the  trailing  vortex  pair. 
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For  modeling  purposes,  the  time  t^  is  taken  to  be  the  end  of  the  jet  regime,  at  which 
time  the  exhaust  occupies  the  entire  vortex  system  and  further  horizontal  growth  is 
curtailed.  The  time 


(2.3) 


is  a characteristic  time  for  the  jet  regime.  It  will  be  shown  in  Section  2.2.  2 that  this 
is  also  a characteristic  time  for  the  vortex  regime,  and  therefore  a convenient  time- 
scaling parameter  to  reduce  the  data  and  extrapolate  results. 

The  velocity  and  temperature  distributions  along  the  jet  are  of  interest.  Reference  2-8 
presents  results  of  an  experimental  study  widely  used  in  jet  calculations.  According 
to  this  study,  the  centerline  velocity  is  given  by 
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where  uc  is  the  centerline  velocity.  According  to  the  same  study,  the  temperature 
is  related  to  the  velocity  according  to 
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0.  7 


The  exponent  0.  7 is  somewhat  larger  than  the  theoretical  value  of  0.  5 predicted  by 
Taylor's  vorticity-transport  theory.  The  cross-sectional  distribution  of  velocity  is 
given  by 


The  concentration  and  temperature  profiles  are  related  to  the  velocity  profile  by  the 
0.  7 power  law  shown  above.  The  preceding  equations  have  been  used  to  calculate  temper- 
ature and  velocity  distributions.  Examples  of  these  are  shown  in  Figures  2-12  and  2-13. 
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Figure  2-12  GEJ85-21  Jet  Centerline  Temperature  Profile 
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Figure  2-13 GE-4  Jet  Centerline  Temperature  Profile 
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2.2.2  Vortex  Regime 


The  vortex  regime  covers  the  section  of  exhaust  trail  extending  from  the  end  of  the 
jet  regime  to  the  point  where  the  trailing  vortex  pair  breaks  up.  This  corresponds  to 
a time  of  the  order  of  100  sec.  As  indicated  in  section  2.  1,  the  vortex  pair  in  this 
regime  dominates  the  dynamics  of  the  wake.  This  vortex  pair  results  from  the  roll- 
up of  the  vortex  sheet  shed  from  lifting  surfaces,  notably  the  aircraft  wing.  It  can 
be  shown  that  the  edge  of  such  a vortex  sheet  is  unstable  and  rolls  up  because  of 
mutual  induction  between  the  vortex  filaments  that  make  up  the  vortex  sheet.  This 
mutual  induction  is  analogous  to  the  electromagnetic  phenomenon,  and,  in  fact, 
obeys  the  same  law  (Biot-Savart).  The  roll-up  of  a trailing  vortex  sheet  has  been 
studied  extensively  (e.  g.  , Ref.  2-7,  and  it  is  a well  documented  fact  that  after  a time 
of  approximately  t/10  the  edges  have  rolled  up  tightly  about  the  two  cores  that  contain 
most  of  the' vorticity  associated  with  the  wing  circulation.  Thus,  at  times  of  order  r , 
that  is, at  the  beginning  of  the  vortex  regime,  the  vortical  wake  of  a finite  wing  is  well 
represented  by  a classical  vortex  pair.  This  pair  consists  of  two  linear  counterrotating 
vortices  with  cores  that  are  small  compared  with  the  aircraft  size  and  set  apart  a 
distance  d of  the  order  of  wing  span.  [ d = (7i/4)  b for  elliptically  loaded  wings.  ] 

The  dynamics  of  a vortex  pair  as  described  above  is  well  understood  (see  Ref.  2-9). 

The  pair  cannot  be  stationary;  due  to  mutual  induction  between  the  vortices,  it  travels 
at  a steady  speed.  The  sinking  speed  is  given  by 


where  T is  the  strength  of  either  vortex.  It  can  be  shown  that  there  is  a mass  of 
recirculating  fluid,  bounded  by  a closed  stream  line,  that  travels  with  the  pair  as  shown 
in  Figure  2-14,  reproduced  from  the  latter  reference.  In  practice,  the  bounding 
streamline  may  sustain  a weak  shear  layer  due  to  velocity  mismatch  between  the 
recirculating  flow  and  the  ambient  fluid;  otherwise,  the  boundary  between  vortex 
and  ambient  fluid  is  sharp.  This  explains  the  containment  effect  that  the  vortex 
has  on  the  horizontal  spread  of  the  exhaust.  For  the  case  of  a trailing  vortex  pair 
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Figure  2-14  Streamlines  of  a Vortex  Pair  as  Seen  by  an  Observer  Moving 
With  the  Pair.  Shaded  area  marks  recirculating  flow  trapped 
in  the  vortex 
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generated  by  a wing,  the  pair  moves  downward  with  a linear  momentum  per  unit  length 

P = pTd  (2.5) 

that  balances  the  wing  lift  (which  is  related  to  T ). 

If  the  vortex  pair  were  to  behave  ideally  as  described  above,  the  aircraft  wake  would 
remain  enclosed  within  its  boundary.  In  practice,  fluid  becomes  detrained  from  the 
vortex  pair  and  deposited  in  a vertical  band.  It  is  believed  that  this  detrainment  is  due 
to  buoyancy  forces,  which  detract  from  the  downward  momentum  of  the  vortex  pair, 
causing  a decrease  of  circulation,  detrainment  of  mass,  and  increase  in  sinking  speed. 

The  vortex  is  buoyant  for  two  reasons;  it  contains  exhaust  gases,  and  it  undergoes  sub- 
sidence heating  as  it  sinks  in  a generally  stable  environment.  An  approximate  computation 
shows  that  the  source  of  buoyancy  is  about  equally  divided  between  exhaust  heat  and  subsi- 
dence, each  contributing  an  increase  in  temperature  of  the  order  of  1°C.  The  upper  central 
photograph  in  Figure  2-1  shows  the  cores  of  the  sinking  vortex  pair  and  the  detrained  ex- 
haust above  them.  This  effect,  which  is  induced  entirely  by  the  aircraft,  has  a lasting 
impact  on  the  wake  because  the  vertical  growth  exposes  it  to  wind  shear,  which  is  one  of 
the  major  mechanisms  for  further  dispersal  in  the  wake  dispersion  regime. 

The  vortex  pair  sinks  at  the  speed  U given  by  Eq.  (2.4),  which  is  relatively  insensitive 
to  buoyancy  effects  because  the  decrease  in  circulation  and  size  tend  to  compensate  each 
other.  In  a perfectly  still  environment,  this  process  would  continue  for  a long  time  because 
the  dissipative  effects  are  small.  It  is  well  known  that  frictional  losses  cannot  account 
for  the  disappearance  of  the  vortex  pair  in  a matter  of  minutes  as  is  usually  the  case.  The 
most  likely  mechanism  for  destruction  of  the  vortex  pair  seems  to  be  the  hydrodynamic 
instability  studied  in  Ref.  2-2.  According  to  this  source,  the  time  for  instabilities  to 
grow  one  e-fold  is 


tB  = 3.  7 T 

where  t is  defined  by  Eq.  (2.3).  A correlation  of  data  from  free  flights  and  laboratory 
experiments  (Ref.  2-3)  shows  that  the  time  for  the  vortex  cores  to  loop  and  start  breaking 

up  is  best  fitted  by 
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(2.6) 


tB  " 10  T 

For  modeling  purposes,  it  will  be  assumed  that  the  vortex  pair  sinks  at  a steady- 
speed  given  by  Eq.  (2.4)  or,  equivalently,  the  wake  in  the  vortex  regime  grows 
vertically  at  the  rate 


dh 

dt 


_8_  W 

3 xt  t_2 

7T  p Vb 


From  the  previous  discussion,  it  is  clear  that  t is  a scaling  time  for  the  vortex 
regime.  Thus,  if  lengths  are  referred  to  as  the  aircraft  span  b and  times  to  t , 
the  last  equation  becomes 


d (h/b)  _ 8 

d(t/r)  3 

7T 

and  a universal  description  of  the  vortex  regime  can  be  constructed  as  follows: 

tA 

Initial  time  “r  = 1.5 

Final  time  — =10 

g 

Vertical  rate  of  growth  — ^ 

7T 

It  was  previously  noted  in  section  2.2.1  that  the  jet  regime  also  scales  with  the  time  r. 
Therefore,  a simple  model  for  jet  and  vortex  regimes  can  be  assembled,  as  will  be 
shown  in  section  2.2.3. 

2.  2.  3 Scaling  Model  For  Jet  and  Vortex  Regimes 

The  jet  and  vortex  regimes  scale  with  the  same  characteristic  time  r and  can  be 
represented  in  a unified  way  as  shown  in  Figure  2-15. 
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Figure  2-15  Scaling  Model  for  Jet  and  Vortex  Regimes 
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The  jet  regime  starts  at  time  zero  and  ends  at  t / r = l.  5 according  to  Eq.  (2.2), 

.A. 

which  defines  the  abscissa  of  point  A in  the  previous  figure.  The  ordinate  of  point  A 
is  the  width  (dimensionless)  of  the  vortex  pair  which,  according  to  Ref.  2-9,  is  given 
by 


w = 2.09d  = 2.09(7r/4)b  = 1.64  b 


or 


for  elliptically  loaded  wings.  In  the  vortex  regime,  the  width  remains  constant  at 
this  value  (line  A B)  whereas  the  height  grows  at  a rate  8/tt  . (See  section  2.2.2) 

The  end  of  the  vortex  regime  is  at  t^/r  = 10  . This  sets  the  final  height  of  the  wake 

ID 

(at  vortex  breakup)  at  w/b  = 3.  83  . It  is  of  interest  to  note  that  the  height  and  width 
of  the  wake  at  vortex  breakup  depend  on  the  aircraft  span  alone.  Figures  2-16  and 
2-17  show  some  data  plotted  in  dimensionless  form  according  to  this  scaling  model. 

These  data  were  obtained  by  LMSC  and  from  Refs.  2-1  and  2-10. 

It  seems  that  in  spite  of  its  simplicity  the  model  is  reasonably  good.  One  problem 
in  testing  the  model  is  the  scarcity  of  well -documented  data.  Flight  and  aircraft 
parameters  must  be  known  to  reduce  physical  data  such  as  that  of  Ref.  2-1  to 
dimensionless  form. 

2.  2.  4 Wake  Dispersion  Regime 

Comparison  of  Data  With  Theoretical  Estimate  — The  wake  dispersion  regime  extends 
from  vortex  breakup  to  times  at  which  the  atmospheric  dispersion  is  independent  of  the 
aircraft  that  deposited  the  exhaust  trail.  Thus,  this  regime  spans  the  gap  between  the 
aircraft -dominated  early  regimes  and  the  purely  atmospheric  processes.  The  wake 
dispersion  regime  is  the  most  complex  of  the  three  presently  defined  regimes  because 
it  involves  several  modes  of  dispersal  and  because  both  aircraft-induced  and  atmospheric- 
induced  processes  coexist  and  are  of  comparable  importance. 
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NOISN3IMIQ 


2-32 


Figure  2-16  Comparison  of  Contrail  Width  Data  With  Scaling  Model 
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Figure  2-17  Comparison  of  Contrail  Height  Data  With  Scaling  Model 


It  was  an  unexpected  finding  of  the  present  study  that  aircraft-induced  effects  last  for 
relatively  long  times.  It  has  been  generally  assumed  that  aircraft  effects  die  out 
rapidly  after  vortex  breakup.  Our  acquisition  of  long  time  data  and  a closer  look  at 
the  major  processes  contributing  to  wake  dispersal  indicate  that  aircraft-induced 
effects  may  be  present  in  significant  amounts  up  to  times  one  order  of  magnitude  larger 
than  vortex  breakup,  that  is,  on  the  order  of  thousands  of  seconds. 

The  prediction  of  the  fluid  mechanical  behavior  of  aircraft  wakes  in  the  dispersion 
regime  involves  the  prediction  of  both  the  wake  height  and  width.  Unfortunately,  clear- 
cut  data  on  long-term  wake  growth  in  both  the  vertical  and  horizontal  directions,  which 
requires  simultaneous  coverage  from  two  camera  stations,  are  not  current^  available. 
However,  combined  width  and  height  data  from  a single  camera  station  have  been  ob- 
tained and  resolved  into  wake  height  and  width  by  means  of  certain  assumptions,  key 
data  points  and  consistency  requirements.  The  data  previously  presented  in  Figure 
2-3  is  a case  in  point.  In  this  case,  good  data  of  the  combined  height  and  width  were 
obtained  from  the  ground  for  times  up  to  11  min  at  a moderate  elevation  angle  (-35  deg) 
and  between  15  and  25  min  at  a low  elevation  angle  (-6  to  8 deg).  In  addition,  several 
excellent  photographs  were  taken  from  the  aircraft  — one  of  the  contrail  height  at  7 min 
and  another  of  the  width  at  33  min.  These  data  were  analyzed  in  the  following  manner. 

For  the  first  11  min,  the  wake  height  was  varied  parametrically,  and  the  wake  width 
was  computed  from  the  observed  combined  height  and  width  data.  It  was  found  that  in 
order  for  the  width  to  increase  monotonically,  the  wake  height  had  to  first  increase 
then  decrease  as  shown  in  Figure  2-3.  Another  requirement  was  that  the  wake  volume, 
computed  as  the  product  of  height  times  the  width,  could  not  decrease.  This  latter 
requirement  is  not  completely  satisfied  at  all  times  for  the  preliminary  data  shown  in 
Figure  2-3,  but  this  slight  deficiency  can  be  corrected  by  additional  massaging  of  the 
data.  The  wake  width  thus  obtained  was  then  faired  into  the  data  point  obtained  from 
the  aircraft  at  33  min.  The  faired-width  data  were  then  used  to  obtain  the  wake  height 
from  the  low  elevation  ground  data  (15  £ t ^ 25  min).  The  resulting  height  data  were 
found  to  follow  a consistent  trend  with  the  assumed  height  variation  for  the  first  11  min. 
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One  interesting  result  of  this  analysis  is  that  the  increase  in  wake  height  in  the  disper- 
sion  regime  scales  with  the  same  characteristic  length  b and  time  t (t  = pVb'Vw)  as 
in  the  vortex  regime.  This  correlation  of  the  wake  height  growth  is  illustrated  in  Figure 
2-18.  The  fact  that  the  wake  growth  correlates  with  the  same  characteristic  length  and 
time  suggest  that  the  wake  height  growth  in  the  dispersion  regime  is  due  to  the  residual 
downward  momentum  of  the  sinking  vortex  pair.  At  the  present  time  it  is  thought  that 
the  increase  in  wake  height  due  to  the  residual  downward  momentum  is  being  opposed  by 
gravitational  forces  (gravitational  collapse  mechanism),  which  exist  when  a well-mixed, 
uni  form -density  fluid  is  placed  in  a stable  atmosphere.  The  densitometer  trace  shown 
in  Figure  2-7  tends  to  support  the  contention  that  the  wake  is  well  mixed  after  vortex -pair 
breakdown.  The  gravitational  force  is  proportional  to  wake  height  while  the  downward 
momentum  per  unit  mass  is  inversely  proportional  to  wake  height,  since  the  wake  is  being 
diluted  as  it  grows.  At  some  point  (t/r  « 60),  the  gravitational  forces  overcome  the 
downward  momentum  and  the  wake  begins  to  collapse.  The  wake  collapse  continues 
until  the  wake  height  is  about  0.  4 of  its  maximum  height. 

Based  on  the  wake  height  variation  observed  for  the  Convair  990,  long-term  B-52  contrail 
data  were  also  analyzed.  These  data  along  with  a theoretical  estimate  of  the  contributions 
of  wind  shear,  buoyant  rise,  gravitational  collapse,  and  wake  and  atmospheric  turbulence 
diffusion  to  the  wake  width  are  shown  in  Figure  2-19.  Details  of  how  these  estimates  were 
obtained  are  discussed  in  subsequent  paragraphs.  Unfortunately,  time  did  not  permit 
similar  theoretical  estimates  to  be  made  for  the  Convair  990  data. 

Regarding  the  data  and  theoretical  estimate  shown  in  Figure  2-19,  it  can  be  seen  that 
wind  shear  is  the  predominant  dispersal  mechanism.  Mass  entrainment  due  to  buoyant 
rise,  which  results  in  an  increase  in  wake  width,  is  seen  to  be  relatively  small.  Wake- 
width  growth  due  to  gravitational  collapse  is  also  relatively  small  and  only  occurs  during 
the  time  the  wake  height  is  decreasing.  The  contribution  due  to  turbulence  deposited  into 
the  wake  by  vortex-pair  breakdown  is  also  small  and  is  confined  to  early  times,  since  the 
wake  turbulence  decays  to  the  atmospheric  level  in  a few  minutes.  Wake  diffusion  due  to 
atmospheric  turbulence  is  a sizable  contributor  to  the  overall  wake  dispersion,  especially 
at  late  times.  Although  there  is  considerable  scatter  in  the  data,  the  agreement  between 
data  and  the  theoretical  estimates  is  encouraging. 
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Figure  2-18  Dispersion  Regime  Wake  Height  Scaling 
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Figure  2-19  Comparison  of  B-52  Contrail  Data  With  Theoretical  Estimates 


Wind  Shear  — Wind  shear  contributes  to  dispersal  of  the  exhaust  trail  by  advection  of 
the  top  and  bottom  of  the  trail  at  different  speeds.  For  simplicity,  the  wind  velocity 
may  be  assumed  linear  through  the  trail  height.  This  assumption  seems  justified  in 
view  of  the  coarseness  and  uncertainty  associated  with  wind  data.  Under  this  assump- 
tion, the  rate  of  growth  of  the  trail  width  is  given  by 


dw 

dt 


h sin  6 


dy 

dz 


(2.7) 


where 

h = trail  weight 

6 = angle  between  the  wind  and  the  trail  axis 
v = wind  speed 
z = vertical  coordinate 

The  correction  for  wind  shear  is  straightforward  and  will  be  easily  implemented  in  the 
mathematical  model.  It  is  assumed  that  appropriate  inputs  of  wind  shear  will  be  made 
available,  a task  that  is  considerably  more  difficult.  As  an  illustration  of  some  of  the 
problems  involved,  Figure  2-20,  reproduced  from  Ref.  2-11,  shows  the  differences  in 
wind  shear  that  may  be  obtained  by  different  manipulations  of  the  same  data. 

An  indication  of  the  relative  importance  of  wind  shear  is  given  in  Figure  2-19  where 
the  contribution  of  wind  shear  to  horizontal  dispersal  was  calculated  on  the  basis  of 
the  local  atmospheric  sounding  at  the  time  the  contrail  was  produced. 

It  can  be  expected  that  stratospheric  flights  will  be  subjected  in  the  average  to  less 
wind  shear  than  the  one  shown  in  the  figure,  which  occurred  in  the  high-shear  region 
near  the  tropopause. 

Wake  Turbulence  Diffusion  — In  the  current  model,  it  is  assumed  that  wake  growth  by 
turbulent  diffusion  in  the  dispersion  regime  is  due  initially  to  turbulence  deposited  in 
the  wake  by  vortex  breakdown.  The  scale  of  turbulence  created  by  vortex  breakdown 
is  assumed  to  be  of  the  order  of  the  wake  width  at  the  time  of  breakdown.  The  turbu- 
lent energy  in  the  wake  at  wavelengths  comparable  to  the  wake  width  is  initially  much 


2-38 


Figure  2-20  Rawin  Measurements,  Washington,  D.  C.  , 14  March  1954,  2100  GCT 
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greater  than  is  present  in  the  free  atmosphere.  As  the  wake  grows  in  width  and  height, 
mixing  with  the  atmosphere  and  viscous  dissipation  will  decrease  the  wake  turbulence 
level  continuously  until  it  attains  the  turbulence  level  in  the  free  atmosphere.  Additional 
turbulence  will  be  deposited  into  the  wake  by  the  turbulent  boundary  layer  and  by  viscous 
shear  and  mixing  at  the  edge  of  the  engine  exhaust  jets.  The  scale  of  this  turbulence  is 
probably  of  the  order  of  1 m or  less  and  therefore  should  not  play  an  important  role  in 
diffusing  the  overall  wake. 

The  kinetic  energy  in  the  wake  due  to  vortex  pair  breakdown  was  estimated  by  equating 

it  to  the  induced  drag  and  by  assuming  that  the  kinetic  energy  was  distributed  uniformly 

over  the  entire  wake  cross  section.  This  estimate  for  the  B-52  flight  of  interest  yielded 

an  average  turbulence  velocity  of  2.4  m/sec.  If  it  is  further  assumed  that  all  of  this 

energy  is  at  the  wavelength  equal  to  the  wake  width,  the  resulting  spectral  density  is 

somewhat  higher  than  that  present  in  what  is  qualitatively  termed  severe  clear  air 

turbulence.  This  estimate  can  be  considered  an  upper  bound,  since  the  turbulence 

would  be  expected  to  distribute  itself  over  some  wavelength  interval.  According  to  the 

findings  of  many  investigators,  severe  turbulence  corresponds  to  a dissipation  rate  of 
2 3 

about  e = 700  cm  /sec  . The  order  of  magnitude  of  the  initial  dissipation  rate  can 

u3 

also  be  estimated  from  the  relationship  e~-r-  > where  u is  the  turbulence  velocity 

ASt 

and  JL  is  the  scale  length,  which  is  assumed  to  be  equal  to  the  wake  width.  Based  on 

the  above  relationship  the  initial  dissipation  rate  is  again  estimated  to  be  about 

2 3 

700  cm  /sec  . 


If  the  dissipation  rate  is  known,  an  estimate  can  be  made  of  the  wake  growth  by  wake 
turbulence  diffusion.  Walton  (Ref.  2-12)  for  example,  has  obtained  a solution  to  the 
diffusion  equation  of  the  form 


<7 


2 


o 1/3  4/3 

+ 2e  (J 
o 


2/3  t2 


(2.8) 


where 

c = cloud  size  = standard  deviation  of  concentration  distribution 

cr  = initial  cloud  size 
o 

e = dissipation  rate 
t = time 
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In  this  solution,  the  diffusion  coefficient  is  taken  to  be 


K 


e 1/3  a4/3  (t) 


(2.9) 


In  the  aircraft  wake,  both  dissipation  rate  and  the  wake  size  are  functions  of  time.  In 
the  free  atmosphere  case,  dissipation  rate  is  taken  to  be  a constant.  An  estimate  of 
the  wake  growth  by  wake  turbulence  can  be  made  by  utilizing  the  above  solution  in  a 
quasi-steady  approximation.  The  approximation  assumes  that  the  rate  of  growth 
(dcr/dt)  for  the  case  e = f (t)  and  cr  = cr  (t)  can  be  approximated  by  solutions  for 
e = constant  and  cr  - a (t). 


The  wake  growth  rate  is  then  given  by  cr 


d a 
dt 


dt  + a 

o 


where  dcr/dt  is  obtained  from  Eq.  (2.8) 

It  now  remains  to  estimate  how  the  wake  dissipation  rate  e will  decay  from  the 

2 3 

initially  high  value  of  the  order  of  700  cm  /sec  to  the  atmospheric  value  which  is  of 
2 3 

the  order  of  1 cm  /sec  . There  are  three  mechanisms  for  decreasing  the  dissipation 
rate  in  the  wake:  (1.)  dilution  of  the  high-energy  air  in  the  wake  by  mixing  with  the  low- 
energy  atmospheric  air,  (2)  an  increase  in  wake  scale  due  to  wake  growth,  and  (3)  in- 
viscid  dissipation  of  large-scale  turbulent  energy. 

A 

o + ^ 

A atm 

where  A is  the  wake  cross-sectional  area  and  the  subscript  o refers  to  conditions 
at  time  of  vortex  breakdown.  The  scale  change  effect  enters  because  the  largest  scale 
size  of  the  turbulent  energy  deposited  by  vortex  breakdown  is  of  the  order  of  the  wake 
width  at  time  of  vortex  breakdown.  Since  the  wake  grows  because  of  various  mechanisms 
and  the  turbulent  scale  of  the  deposited  energy  does  not,  the  dissipation  rate  will  decrease 

inversely  as  the  wake  width. 


The  dilution  effect  can  be  estimated  by 


€ = e 

o 
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+ e 


atm 


The  inviscid  dissipation  of  large-scale  turbulent  energy  to  smaller  scales  will  also 
decrease  the  dissipation  rate.  If  we  assume  that  the  large  scale  eddies  lose  1/e  of 
their  kinetic  energy  in  the  characteristic  time  £/ u,  a rough  estimate  of  the  turbulent 
velocity  decay  can  be  obtained  by  means  of  a stepwise  evaluation  of  the  characteristic 
time.  The  dissipation  rate,  including  all  of  the  above  decay  mechanisms,  is  then 
estimated  from 


e 


e 


atm 


The  evaluation  of  the  above  estimate  of  the  dissipation  rate  Figure  2-21  indicates  that 

the  wake  turbulence  will  decay  to  the  atmospheric  value  in  about  2 to  3 min  for 

2 3 2 3 

e = 700  cm  /sec  and  in  about  1 to  2 min  for  e = 100  cm  /sec  . 
o o 


The  wake  and  atmospheric  turbulence  contribution  to  wake  growth  was  evaluated  using 
the  dissipation  rate  decay  as  estimated  above  in  the  quasi-steady  approximation 
described  previously.  The  results  of  these  calculations  are  shown  in  Figure  2-19. 

2 

The  difference  between  the  results  for  assumed  initial  values  of  e of  700  and  100  cm  / 

3 ° 

sec  is  insignificant  because  of  the  rapid  decay  of  dissipation  rate  to  the  atmospheric 

value.  Diffusion  of  the  wake  a few  minutes  after  vortex  breakdown  is  due  entirely  to 

atmospheric  turbulence. 


In  the  present  case,  the  wake  diffusion  is  scale-dependent  according  to  Eq.  (2.  9). 
Partial  justification  for  this  form  of  the  diffusion  coefficient  relies  on  the  comparison 
of  predicted  growth  rates  with  observed  growth  rates  of  nuclear  debris  clouds  (Ref. 
2-12)  and  similar  data.  These  data  may  include  the  combined  effect  of  vertical  trans- 
port and  wind  shear.  In  other  words,  the  scale  dependent  diffusion  prediction  may  in- 
clude the  effects  of  wind  shear.  If  this  is  the  case,  the  contribution  of  wake  and  atmos- 
pheric turbulence  to  the  wake  growth  may  be  less  than  that  currently  estimated. 
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DISSIPATION  RATE,  e (CM  /SECT) 


Figure  2-21  Estimate  of  Wake  Dissipation  Rate  Decay 
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Gravitational  Wake  Collapse  — When  a well-mixed  region  is  created  locally  in  a stably 
stratified  environment,  the  fluid  in  the  well-mixed  region  tends  to  acquire  neutral 
stability  while  the  surrounding  fluid  remains  stable.  As  a consequence  local  buoyancy 
forces  are  created,  with  the  top  of  the  mixed  region  becoming  negatively  buoyant  and 
the  bottom  positively  buoyant  with  respect  to  the  environment.  This  creates  a vertical 
collapse  of  the  well-mixed  region,  which  flattens  and  spreads  horizontally  as  every 
parcel  of  mixed  fluid  seeks  its  own  density  level  in  the  stratified  environment. 

A number  of  laboratory  experiments  have  been  carried  out  to  study  these  effects  in 
water  tanks  containing  stratified  solutions  (Refs.  2-13,  2-14,  and  2-15,).  It  has  been 
observed  that  mixed  regions  created  with  various  agitators  or  by  the  turbulent  wakes  of 
self-propelled  bodies  collapse  as  described  above.  The  degree  of  collapse  and  the  final 
shape  depend  on  the  initial  degree  of  stratification  and  the  thoroughness  of  the  mixing- 
process.  The  experiments  generally  involve  much  larger  stratification  than  that  of 
the  standard  stratosphere  and,  except  for  those  of  Ref.  2-15,  are  limited  to  relatively 
short  times.  Therefore,  application  of  these  results  to  aircraft  wakes  involves  sub- 
stantial extrapolation.  Nevertheless,  a certain  amount  of  wake  collapse  should  be  ex- 
pected in  stratospheric  wakes,  and  the  extrapolation  is  justified  within  the  error  bounds 
of  the  present  investigation. 

The  water-tank  experiments  mentioned  above  show  that,  in  the  wake-collapse  process, 
the  length  scale  is  given  by  the  size  of  the  body  generating  the  mixed  region,  and  the 
time  scale  is  given  by  the  Brunt- Vaisala  period 


1 


The  theoretical  estimate  of  the  increase  in  wake  width  due  to  gravitational  collapse 
shown  in  Figure  2-19  is  based  on  conservation  of  mass.  As  the  wake  height  decreases, 
the  wake  width  must  increase  to  compensate  for  the  decrease  in  wake  volume  due  to 
collapse.  During  collapse,  other  mechanisms  are  present  which  simultaneously  in- 
crease the  wake  width.  Since  the  wake  height  decrease  for  the  B-52  was  estimated 
from  preliminary  Convair  990  data,  the  estimate  of  the  wake  growth  by  gravitational 
collapse  is  somewhat  uncertain. 
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Buoyant  Rise  - Buoyant  rise  of  the  exhaust  trail  has  been  proposed  as  a growth  mecha- 
nism at  late  times  (Ref.  2-16).  The  growth  results  from  entrainment  of  ambient  air, 
through  turbulent  mixing.  The  source  of  buoyancy  in  the  wake  is  the  thermal  energy 
in  the  engine  exhaust,  and  subsidence  heating  due  to  the  vortex  motion.  A quick  calcu- 
lation shows  that  the  temperature  excess  in  the  wake  is  not  large  but,  on  the  other  hand, 
small  temperature  differences  can  produce  significant  buoyant  rise.  In  Ref.  2-16,  the 
buoyant  rise  and  associated  rate  of  growth  are  calculated  using  chimney  plume  theories. 
In  the  present  investigation,  the  problem  was  approached  from  a different  point  of  view, 
namely,  by  considering  the  exhaust  trail  as  a line  thermal  deposited  instantaneously. 
Initially  the  thermal  is  assumed  to  have  a finite,  circular  cross  section  and  a finite 
amount  of  buoyancy  per  unit  length,  represented  by  the  temperature  difference  above 
ambient.  At  time  zero  the  thermal  is  allowed  to  rise  freely  in  an  environment  of 
arbitrary  stratification.  The  theoretical  approach  of  Ref.  2-17  for  instantaneous  point 
sources  was  followed  with  the  proper  reformulation  for  the  case  of  a line  source.  The 
basic  assumptions  are  that  the  rate  of  turbulent  entrainment  is  proportional  to  the 
speed  of  rise  through  an  entrainment  coefficient  a , and  that  an  incompressible  formu- 
lation can  be  used  for  the  atmosphere  provided  potential  rather  than  physical  density 
is  used  throughout. 

The  differential  equations  expressing  conservation  of  volume,  mass,  and  momentum 


dr 


2 


dt 


2 r a u 


(2.  11) 


(2.  12) 


(2.  13) 


where  r is  the  radius,  P the  potential  wake  density,  and  p the  potential  ambient 

ci 


density,  an  arbitrary  function  of  altitude.  This  system  can  be  solved  in  closed  form 
in  terms  of  quadratures  containing  the  arbitrary  ambient  density.  For  present 
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purposes  it  was  simpler  to  calculate  numerical  solutions  to  the  differential  equations 
with  a high-speed  computer.  The  entrainment  coefficient  suggested  by  Ref.  2-17 
(Q!  « 0.  3)  was  used,  and  the  buoyant  rise  was  calculated  from  time  zero  to  the 
height  of  first  velocity  reversal.  The  initial  temperature  difference  was  given  as  an 
input  parameter.  Estimates  of  this  difference  were  made  (conservatively)  by  assign- 
ing the  entire  jet  engine  energy  to  the  temperature  rise  in  a cylindrical  wake  having 
uniform  temperature  and  the  same  cross  section  as  the  observed  wake.  In  an  effort 
to  account  for  all  sources  of  heat,  the  average  temperature  rise  due  to  subsidence 
was  also  calculated  and  added  to  the  previous  one.  The  results  are  shown  in  Figure  2-22 
in  the  form  of  growth  histories.  It  can  be  seen  that  for  the  B-52  wake  the  temperature 
increase  due  to  the  exhaust  alone  is  0.2°C,  and  including  average  subsidence  it  is 
0.  5°C.  The  maximum  is  about  1°C.  For  the  case  of  an  SST  flying  at  20  km  in  a 
standard  stratosphere,  the  temperature  difference  is  larger  (due  to  the  high-power 
output,  particularly  with  afterburner)  but  the  wake  is  smaller  and  the  ambient  more 
stable,  resulting  in  total  growth  comparable  to  that  of  the  B-52. 
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Figure  2-22  Growth  Due  to  Buoyant  Rise  of  Air  Parcels  Equivalent  to  B-52  and  SST 
Exhaust  Trails 
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Section  3 

CHEMICAL  KINETICS  MODEL 


3.1  SUMMARY 

The  present  investigation  is  directed  at  the  nature  and  extent  of  chemical  reactions  in 
the  exhaust  jets  and  wakes  of  high- altitude  aircraft.  The  object  is  to  determine  what 
reactions  might  have  an  important  effect  on  the  composition  of  hot  gases  as  they  issue 
from  the  engine  into  the  surrounding  stratosphere.  Experimentally  measured  exhaust 
compositions  at  the  exit  plane  of  an  engine  will  not  reflect  accurately  the  final  chemical 
composition  of  the  exhaust  in  the  stratosphere  if  reactions  in  the  wake  are  important. 

The  objective  of  Task  I,  Problem  Definition,  is  to  identify  essential  features  or  con- 
siderations needed  for  subsequent  development  of  a wake  flow  and  chemistry  model  in 
order  to  define  species  concentrations  and  wake  configurations  at  a time  in  the  wake 
history  when  air  craft- induced  motions  have  diminished  and  are  dominated  by  the  natural 
stratospheric  motions.  The  technical  approach  for  the  wake-chemistry  portion  of 
Task  I has  been  to  screen  the  many  possible  chemical  reactions  to  determine  their 
probable  importance  in  the  engine  exhaust  jet,  the  near  wake,  and  the  far  wake  — i.e.  , 
to  as  late  as  20  min  after  passage  of  the  aircraft.  To  examine  the  reactions  in  the 
exhaust  jet,  however,  it  has  been  necessary  for  the  calculations  to  cover  areas  even 
farther  upstream  — i.e.,  within  the  engine  nozzle.  This  has  been  necessary  to  examine 
the  rate  histories  for  many  of  the  species  before  they  reach  the  nozzle  exit  so  that  their 
concentrations  at  the  exit  plane  can  be  established.  Many  atoms  and  radicals  are  found 
to  be  present,  and  it  cannot  be  reasonably  anticipated  that  experimental  test  data  on 
representative  turbojet  engines  will  be  capable  of  quantitatively  identifying  all  of  the 
important  species.  Thus,  the  engine  nozzle  kinetics  work  is  a necessary  part  of  the 
chemical  reaction  screening  process.  In  this  approach,  simplifying  relationships 
have  been  sought,  such  as  steady  state  and  partial  equilibrium  relationships  among 
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minor  species , to  discover  when  such  relationships  can  be  used  for  making  reliable 
predictions . The  reliability  of  the  simplifying  methods  is  ascertained  by  comparison 
of  the  results  with  fully  coupled  computer  solutions  of  the  reacting  flow  problem. 

The  emphasis  in  final  development  of  the  chemistry  model  will  be  influenced  by  the 
interface  with  and  requirements  of  the  chemical  modeling  program  for  the  post-wake 
dispersion  regime.  It  is  anticipated  that  the  NO^  concentrations  and  total  NO^  will 
be  first  on  the  list  of  priorities . It  is  anticipated  that  any  significant  conversion  of 
NO  in  the  wake  to  HNO  will  be  of  major  interest.  Formation  of  sulfate  aerosol  in 

X o 

the  wake  and  diffusion  and  transport  regime  will  also  be  of  interest.  It  seems  likely 
that  the  diffusion  and  transport  regime  and  subsequent  climatic  analyses  will  be  less 
sensitive  to  uncertainties  in  wake  CO  and  hydrocarbon  levels  than  to  NO^  concentration. 
In  light  of  this,  a simplified  wake-chemistry  model  may  serve  nicely.  If,  on  the  other 
hand,  detailed  information  on  minor  wake  species,  exhaust  generated  radicals,  etc.  , 
is  required,  a much  more  general  and  comprehensive  wake  kinetics  model  will  be 
needed. 

The  scope  of  the  chemical  screening  process  is  broad;  i.e.  , thermochemical  reactions, 
photochemistry,  and  particulate  generating  reactions  are  all  considered.  Because 
oxides  of  nitrogen  are  thought  to  play  a particularly  important  role  in  the  stratosphere, 
their  reactions  are  given  special  attention,  and  analysis  of  the  photochemical  catalytic 
NO  and  ozone  history  following  passage  of  an  aircraft  through  the  stratosphere  is 

X 

carried  out  to  times  well  beyond  the  20-min  "end  to  wake"  time. 

For  studies  of  thermal  reactions  in  the  engine  exhaust  regime,  the  anatytical  approach 
involves  setting  up  a chemical  kinetics  model  using  those  elementary  reactions  that 
play  a role  in  changing  the  concentrations  of  engine  emissions  in  the  exhaust  nozzle 
and  the  near  jet  downstream  from  the  engine  exit  plane.  Numerical  methods  are  used 
to  integrate  the  resulting  governing  chemical  rate  equations  up  to  the  point  where  the 
exhaust  flow  mixes  with  the  ambient  atmosphere.  Initial  conditions  for  the  computations 
are  compatible  with  available  data  on  engine  exhaust  composition. 
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In  particular,  reactions  affecting  the  oxidation  of  CO  to  CO  , NO  to  NO  , and  SO 

LA  Li  LA 

to  SOg,  as  well  as  conversion  of  N02  to  HNO^,  are  of  interest.  A model  consisting 
of  elementary  forward  and  backward  reactions  pertinent  to  these  exhaust  chemical 
systems  has  been  used  in  this  study.  Reaction  rate  constants  are  taken  from  recent 
evaluations  available  in  the  literature.  The  fluid-dynamic  structure  of  the  exhaust  jet 
region  was  obtained  from  analytical  predictions  and  experimental  data  available  from 
the  literature. 

The  major  conclusions  reached  in  the  investigation  are  as  follows: 

• In  an  afterburning  engine,  with  relatively  low  HC  emissions,  nonequilibrium 
effects  are  an  important  aspect  of  the  chemistry  in  the  engine  exhaust  expan- 
sion nozzle  and  exhaust  jet  regime.  There  is  substantial  oxidation  of  CO  to 
C02  in  the  nozzle  and  jet,  and  the  concomitant  production  of  H atoms  main- 
tains the  concentrations  of  H,  O,  and  OH  well  above  their  equilibrium  values. 
The  partial  equilibrium  of  some  bimolecular  reactions  are  the  basis  for  use- 
ful relationships  among  the  concentrations  of  H,  O,  OH,  and  H9.  However, 
these  radicals  cannot  be  related  to  major  species  because  recombination  re- 
actions are  far  out  of  equilibrium. 

• Hydrocarbons  have  a profound  scavenging  effect  on  exhaust  inorganic  radicals. 
Introduction  of  HC  in  quantities  of  the  same  order  as  CO  reduce  the  H,  OH, 
and  O concentrations  by  many  orders  of  magnitude. 

• Oxidation  of  NO  to  NO  is  not  significant  in  the  exhaust  jet.  If  unburned  hy- 
drocarbons  are  not  present,  the  N02/N0  ratio  is  maintained  well  below  its 
equilibrium  value  by  fast  reactions  such  as 

NO  + H -*  NO  + OH 

La 

N02  + O — NO  + 02 

With  HC  present  and  inorganic  radicals  consequently  few,  the  effect  is  to 
freeze  the  N09  in  the  nozzle  near  its  throat  value. 

La 

• In  the  jet  mixing  region  both  reactions 

NO  + O + M — N02  + M 

no9  + oh  — hno3 
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are  slow  relative  to  OH  recombination  via 


OH  + OH  + M — H 09  + M 

Consequently,  there  is  no  regime  of  sufficient  duration  in  which  OH  and  N09 
are  simultaneously  plentiful,  and  stable  nitric  acid  is  not  formed  in  significant 
quantities. 

• Oxidation  of  SO  to  SO  in  the  hot  exhaust  jet  is  forestalled  by  H and  O atom 

Z U 

attack  on  S0o  in  a manner  similar  to  the  N0-N02  system. 

• In  the  jet  mixing  region  and  the  near  wake,  the  possibility  of  aerosol  forma- 
tion by  the  following  reaction  pair  can  be  ruled  out: 

S09  + O + M — S03  + M 

s°3  + h20  - h2so4 

Other  reactions  that  might  lead  to  particulate  formation  appear  to  be  too  slow 
to  contribute  significantly  in  the  wake  regime. 

• In  ground  tests  the  exhaust  gases  are  maintained  at  substantially  higher  tem- 
peratures for  longer  residence  than  at  high  altitude.  As  observed,  the  chem- 
ical kinetics  model  predicts  the  occurrence  of  significantly  more  CO  oxidation 
in  the  exhaust  jet  for  an  engine  under  static  test  conditions  at  sea  level,  than 
it  predicts  for  stratospheric  flight. 

• The  subject  of  chemical  kinetics  in  turbulent  flows  is  considered.  To  a large 
extent,  chemistry  is  found  to  be  "diffusion  limited"  in  the  exhaust  jet  core 
and  "reaction  rate  limited"  in  the  wake.  The  simple,  "well-mixed"  kinetics 
model,  which  ignores  diffusion  and  turbulence  in  the  exhaust  jet,  appears  to 
agree  reasonably  well  with  ground-test  data  on  CO  profiles.  This  may  mean 
that  the  net  influence  of  turbulence  on  the  chemistry  is  small. 

• Addition  of  unburned  exhaust  hydrocarbons  to  the  stratosphere  does  not  appear 
to  be  significant.  Catalytic  chain  effects  of  hydrocarbons  on  ozone/NO^,  cycles 
degrade  the  hydrocarbons  and  are  unimportant.  Removal  of  NO  r by 
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hydrocarbons  in  the  stratosphere  is  an  important  area  for  study,  but  the 
contribution  from  engine  exhaust  hydrocarbons  is  dwarfed  by  the  ambient 
methane. 

• The  release  of  chemical  enthalpy  due  to  reaction  of  exhaust  species  is  too 
small  to  have  an  important  effect  on  the  fluid  mechanics.  Thus,  the  chem- 
istry may  be  uncoupled  energetically  from  the  flow  field  calculations. 

• The  species  of  primary  interest  in  evaluating  possible  climatic  effects  of 
high-altitude  aircraft  — NO  - is  not  removed  or  depleted  by  gas-phase 

X 

chemical  reactions  in  the  aircraft  wake  regime. 


3.2  CARBON  MONOXIDE  OXIDATION  - EQUILIBRIUM  OH  TREATMENT 


3.2.1  Major  Oxidation  Reaction 

A two- reaction  zone  model  is  generally  accepted  as  the  most  useful  for  kinetic  interpre- 
tation (Ref.  3-1).  The  gas  turbine  combustor  exhibits  a two-reaction  zone  behavior 
wherein  CO  is  formed  in  the  regions  near  the  fuel  nozzle  and  further  oxidation  to  C09 
takes  place  downstream  of  this  point.  The  degree  of  conversion  to  C09  is  strongly 
dependent  upon  equivalence  ratio. 

The  accepted  rate  controlling  step  in  CO  oxidation  in  hydrocarbon  systems  is  the 
reaction  (Refs.  3-1  and  3-2): 


CO  + OH  — C09  + H (3.1) 

The  reactions  CO  + O — C09  and  CO  + 09  — C09  + O are  much  too  slow  (Ref.  3-1). 
3.2.2  Kinetic  Details 

In  general,  it  is  necessary  to  consider  the  reaction  mechanism  including  all  reactions 
affecting  OH  and  H to  arrive  at  the  proper  OH  concentration.  This  reaction  scheme 
(Ref.  3-3)  is  as  follows: 


CO  + OH  ~ 

H + °2  ^ 
OH  + OH  ^ 

° +H2  - 
H2  + OH 

H2  + M ^ 

°2  + M ^ 
H + OH  — 


co2  + H 

OH  + O 

h9o  + O 

H + OH 
H20  + H 
H + H + M 
O + O + M 
H20  + M 
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3.2.3  Overall  Kinetics 


In  those  situations  in  which  one  slow  reaction  determines  the  overall  reaction  rate, 
the  overall  or  global  rate  approach  is  often  quite  useful,  and  can  circumvent  the  rigors 
of  integrating  a complete  reaction  scheme  (Ref.  3-1).  The  oxidation  of  CO  controlled 
by  Eq.  (3. 1)  provides  just  such  a situation,  and  many  authors  have  given  overall  ex- 
pressions for  CO  oxidation.  These  various  expressions  are  reviewed  by  Singh  and 
Sawyer  (Ref.  3-2)  and  will  not  be  discussed  here. 


The  major  problem  with  empirically  determined  overall  rate  equations  is  that  they 


may  only  be  valid  for  the  particular  conditions  under  which  they  were  measured.  One 
approach  to  overall  CO  oxidation  is  to  consider  partial  equilibration  of  some  of  the 
reactions  above. 

In  post-combustion  zone  mixtures  of  hot  gases  with  excess  O 2 and  H90  , the  assump- 
tion of  OH  equilibrium  is  quite  good  (Ref.  3-1).  Gas-turbine  exhaust  streams  typically 
meet  these  conditions,  and  the  time  involved  in  passing  from  the  combustor  to  the  exit 
will  allow  OH  equilibration.  Using  the  equilibrium  relations 


2 


O + H20  = OH  + OH  , K 


(3.2) 


0 + 0 = 0,  Kp  = 

2 2 (PQ) 


2 


(3.3) 


2 


the  OH  concentration  may  be  written  as 


[OH] 


(3.4) 
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l/4 

[The  (RT)  factor  results  from  the  use  of  equilibrium  constants  in  terms  of 
pressure. ] 

The  equilibrium  constants  Kp^  and  Kp9  are,  of  course,  well  known  from  thermo- 
chemical tables  and  can  be  fitted  to  an  equation  of  the  form 


K 


P 


AT 


a 


-E/RT 

e 


The  latest  recommended  value  for  in  Eq.  (3. 1)  is  (Ref.  3-4) 


k 


1 


4.2  x 10U 


-1080/RT  3 . -1  -1 

e cm  -mole  -sec 


Using  the  foregoing  k^  expression,  the  fitted  forms  for  in  Eq.  (3.4),  combining 

with  Eq.  (3.  5)  for  the  disappearance  of  CO  according  to  Eq.  (3. 1),  and  eliminating  a 
“0  218 

weak  T ‘ dependence  through  its  evaluation  at  T = 1500 °K  , in  the  range  of  in- 
terest, one  obtains  Eq.  (3.6). 

= -k1  [CO]  [OH]  (3.5) 


11221  = _io12-6  e 39800/Rr  [CO]  [H90] 1/2  [O.J  1/4  (3.6) 

The  H.;0  and  09  concentrations  are  given  by  their  equilibrium  values.  The  example 
of  CO  oxidation  for  given  profiles  follows. 

To  compute  the  relative  change  in  the  CO  mole  fraction  when  jet  engine  exhaust  gases 
expand  over  specified  temperature  and  pressure  profiles,  one  may  recast  Eq.  (3.6) 
and  integrate  stepwise  according  to  Eq.  (3.7): 
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At 


/ 


e 


-39800/RT  0.75  0. 25  v0.25 
p XH2o  xo2 


(3.7) 


o 


where 


X = mole  fraction 

p = gas  density 

M = average  molecular  weight  (assumed  constant) 

At  = time  step  taken  during  integration 

( ) = initial  value 

v ' f\ 


o 


Additional  simplification  usually  results  because  the  mole  fractions  of  HgO  and  O 2 
are  frequently  almost  constant  in  the  nozzle,  and  these  terms  have  small  exponents 
in  Eq.  (3.7). 

For  the  GE-4  engine  at  maximum  power,  Eq.  (3.7)  has  been  integrated  numerically 
for  temperature,  pressure,  time,  and  distance  relationships  obtained  from  fluid 
mechanics  analysis.  These  relationships  appear  in  Figure  3-1.  Equilibrium  fractions 
of  H20  and  02  at  the  secondary  nozzle  throat  were  computed  for  the  overall  equiva- 
lence ratios  given  by  the  following  engine  manufacturer's  data  for  the  corresponding 
power  setting; 


The  results  appear  in  Figure  3-3.  Clearly,  the  centerline  CO  mole  fraction  stops  de- 
creasing almost  immediately  after  the  nozzle  expansion  begins.  Examination  of  the 
strong  temperature  dependence  of  the  CO  reaction  rate  in  Eq.  (3.  7)  shows  why  this 
result  is  obtained,  in  view  of  the  temperature  histories  shown  in  Fig.  3-1.  As  will 
be  seen,  however,  this  very  severe  freezing  of  CO  is  not,  in  fact,  obtained.  Thus, 


Altitude  Mach  Power 
Condition  (ft)  Number  Setting 


Maximum 

Power  65,000 


2.7 


1 
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£ 

w 

o 
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Q 
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Figure  3-1  GE-4  Secondary  Nozzle  Flow  Properties 


CO 

^co^, 


Figure  3-2  CO  Concentration  Ratio  Profile  for  GE-4  at  Maximum  Power 
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this  equilibrium  OH  assumption  is  not  justified  for  the  engine  exhaust  analysis,  in 
light  of  results  from  detailed  chemical  kinetics  numerical  analysis.  Actually,  the  OH 
is  maintained  far  in  excess  of  its  equilibrium  concentration  by  a chain  sequence  initi- 
ated by  H atoms  generated  in  CO  oxidation  by  the  reaction  shown  in  Eq.  (3. 1). 
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3.3  NITRIC  OXIDE 


3.3.1  General 

Nitric  oxide  emissions  from  gas  turbine  combustors  are,  in  general,  much  less  than 
those  predicted  by  thermochemical  equilibrium  considerations.  Experimental  meas- 
urements (Ref.  3-5)  and  theoretical  kinetic  models  (Ref.  3-6)  have  established  that 
the  NO  is  formed  in  the  high-temperature  post-flame  gases  in  the  primary  zone  of 
the  combustor.  The  residence  time  in  this  zone  is  not  long  enough  for  the  NO  to 
reach  its  equilibrium  concentration,  and  rapid  quenching  in  the  secondary  zone  freezes 
the  NO  concentration  at  its  subequilibrium  level  (Ref.  3-6).  Both  Heywood  et  al. 

(Ref.  3-7)  and  Caretto  et  al.  (Ref.  3-8)  have  pointed  out  that  the  NO  should  not  de- 
compose when  its  concentration  is  less  than  the  local  equilibrium  value. 


In  a gas  turbine  engine  with  an  afterburner,  there  is  the  possibility  of  additional  NO 
formation  in  the  afterburning  region.  However,  it  is  expected  that  the  short  residence 
times  and  lower  temperatures  in  this  region,  as  compared  with  the  primary  zone,  will 
prevent  significant  NO  formation  in  the  afterburner. 

3.3.2  Kinetics 

The  accepted  mechanism  for  NO  formation  in  a high-temperature  combustion  system 
is  (Refs.  3-7  and  3-8): 


N + NO  — N2  + O 


(3.8) 


N + O, 


2 


NO  + O 


(3.9) 


N + OH  ^ NO  + H 


(3.10) 


H + N90  — N2  + OH 


(3.11) 
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(3.12) 


° + N2°  ^ N2  + °2 


o + N20  ^ 2 NO 


(3.13) 


The  first  two  reactions  represent  the  well  known  Zeldovich  (Ref.  3-9)  mechanism  and 
are  generally  sufficient  to  account  for  NO  production  in  high-temperature  combustion 
systems.  The  last  three  reactions  can  be  important  at  lower  temperatures  in  lean 
mixtures.  The  [N2J  , [O]  , [02J  , [OH],  and  [H]  are  given  by  their  equilibrium 
values,  and  if  [N]  and  [N90]  are  given  by  their  steady- state  values,  then  the  mole 
fraction  of  NO  can  be  shown  to  be  given  by  (Ref.  3-6) 


dX^^/dt  | is  evaluated  for  initial  conditions  in  the  GE-13  jet  expansion  case 
previously  given, 


Thus,  for  the  time  scale  of  the  jet  expansion,  there  will  be  no  change  in  NO  mole 
fraction  due  to  the  reaction  scheme  (3. 8)  to  (3.  13). 

3.3.3  Other  Reactions 

The  question  arises  as  to  the  possibility  of  decomposition  of  NO  in  the  outer  cooler 
portion  of  the  jet  where  [NO]  exceeds  its  equilibrium  value. 


dX. 


= 0. 113  x 10  5/sec 


{a  « 1) 


o 
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Laurendeau  and  Sawyer  (Ref.  3-10)  have  considered  NO  decomposition  due  to  the 
reaction  scheme: 


NO  + N — ■ 0 + N2 

(3.  14) 

NO  + O — N + 02 

(3.15) 

no2  + m ^no  + o + m 

(3.  16) 

NO  + 02  ^ O + N02 

(3.  17) 

NO  + NO  — N20  + 0 

(3.18) 

n2o  + m — n2  + o + m 

(3.19) 

Starting  with  an  initial  mixture  of  10%  NO  and  90%  inert,  they  computed  the  NO  de- 
composition at  1700°K  and  0.7  atm.  Their  results  show  that  the  time  scale  for  the 
beginning  of  significant  NO  decomposition  is  on  the  order  of  10  sec  under  these  con- 
ditions. N09  levels  resulting  from  this  reaction  scheme  are  3 to  4 orders  of  magni- 
tude less  than  the  NO  levels.  Thus,  in  the  GE-13  exhaust  jet  considered  previously, 
for  example,  where  « 0.1  and  T < 1700°K,  the  decomposition  of  NO  due  to 

reaction  (3.  14)  to  (3.  19)  will  be  insignificant. 

Possible  reactions  involved  in  converting  NO  to  NO.-,  in  the  outer  jet  are: 

NO  + O + M — N02  + M 
NO  + O — N09  + hy 

(3.20) 

(3.21) 

NO  + NO  + 02  — N09  + N02 

(3.22) 

no2  + no2  — no  + no  + o2 

(3.  23) 

NO  + NO  — N02  + N 

(3.24) 
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A discussion  of  these  reactions  and  their  respective  rate  constants  is  given  in 
Refs.  3-11  and  3-12.  Reaction  Eq.  (3.23)  is  insignificant  at  temperatures  less  than 
5000 JK  and  can  be  eliminated  immediately.  The  rate  constants  for  reactions  (3.  20) 
to  (3. 23)  are: 


k3.21  = 3'9  X 10' 


/1870 

exp 

\ RT 

( T n 

-1.94 

\ 

\296, 

) 

3 , -1  -1 

cm  -mole  -sec 


k3  99  = 2.4  x 109  e+1046/RT  Cm6-mole  2-sec  1 


, . 1A12  -26900/RT  cm2-mole  ^-sec  1 

k3.23  - 4 x 10  e 


Considering  the  low  temperatures  and  very  small  equilibrium  O atom  concentrations  in 
the  outer  jet  region,  reaction  (3.22)  will  dominate  N09  formation  in  this  region.  The  NO 
decomposition  rate  under  this  scheme  is  thus  given  by 


d [NO] 
dt 


..  0 ..9  1046/RT 

(4.  8 x 10  e 


) [NO]"[09] 


or  in  terms  of  mole  fractions 


dX 


NO 


dt 


= (4.  8 x 


1A9  1046/RT 

10  e ) 


ifi Lx2  X 
M XNOXOr 


Taking  the  last  data  point  in  the  GE-13  jet  expansion  as  typical  outer  jet  conditions 

gives: 

p = 0.37  x 10  4 g/cm2 
T = 950  °K 


dX 


NQ  _ i 4 x 10-2  X2  X 
dt  1,4  10  XNOXOr 
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Considering  that  typical  magnitudes  for  the  mole  fractions  are 


XNO  < °'001 


< 0.2 


“9 

dX^/dt  will  be  of  order  10  /SEC.  Thus,  there  will  be  no  significant  change  in  NO 
mole  fraction  in  the  jet  due  to  formation  of  NO„. 

The  overall  conclusion  is  that  based  on  the  equilibrium  assumption,  the  mole  fraction  of 
NO  remains  constant  in  the  exhaust  jet.  A complete  analysis  presented  in  the  following 
sections  demonstrates  that  equilibrium  is  not  attained  and  thus  invalidates  this  analysis. 
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3.4  CARBON  MONOXIDE  OXIDATION  AND  NOx  KINETICS  - COMPLETE  KINETIC 
TREATMENT 

3.4.1  Rationale 

The  equilibrium  analysis  utilized  in  section  3.2  is  questionable  in  the  light  of  observed 
high  OH  concentrations  in  jet  engine  exhausts  (Ref.  3-13).  Therefore,  to  further  in- 
vestigate changes  in  CO  and  NO  concentrations  in  the  wake  region  of  high-altitude 

X 

aircraft,  simplified  treatments  had  to  be  abandoned.  It  was  necessary  to  set  up  a 
chemical  kinetics  model  including  all  elementary  reactions  which  might  conceivably 
contribute  to  the  alteration  of  emission  levels. 

3.4.2  Reaction  System 

• CO  oxidation  to  C09  has  only  one  important  elementary  reaction  (Refs.  3-1, 

3-2,  3-3: 


CO  + OH  ~ C09  + H 


• This  necessitates  inclusion  of  all  reaction  controlling  OH  and  H levels: 


h + o2 

OH  + OH 

° + H2 
H2  + OH 


OH  + O ' 


h2o  + O 

H + OH 
H90  + H 


fast  bimolecular 
radical  shuffle 
reactions 


H + H + M 
O + O + M 
H + OH  + M 


H2  + M 
09  + M ► 
Ho0  + M 

Z 


3 body  association 
reactions 
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• Reactions  involving  H09  also  affect  the  pool  of  radical  species: 


H + C>2  + M 
H + HO„ 


OH  + HO_ 


O + HOr 


H02  + M 
OH  + OH 
H2°  + o2 
OH  + 0„ 


H + HO, 


H + HO, 


H2  + °2 
H2°  + ° 


• H209  Reactions: 


OH  + OH 

+ M 

H2°2  + M 

ho2  + 

H°2 

— 

H2°2  + °S 

H2°2 

+ OH 

H20  + HO. 

H + 

H2°2 

- 

H2  + H02 

Nitric  acid  reactions: 

N02 

+ OH 

- 

HNO, 

Sulfur  reactions: 

so2  + o 

+ M 

SO,  + M 

H + 

S03 

SO,  + OH 

0 + 

so3 

— 

so2  + o2 
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• Principal  reactions  exchanging  NO  and  N02: 


NO  + O + M ZT  N09  + M 
H + N02  ^ NO  + OH 

O + N09  H NO  + 0.? 

NO  + H09  rz  N09  + OH 

• Other  reactions  involving  NO,  NO  N , and  N O: 

La  Lt  La 


0 + N2 

NO  + N 

Zellovitch  NO  formation, 
important  at  high 

M + 02 

— 

NO  + 0 

temperature  only 

NO  + NO 

n2o  + 0 

H + NO 

£a 

— 

N2  + OH 

0 + N20 

N2  + 02 

n2  + no9  r;  no  + n9o 
n9o  + m t:  n9  + o 
n2  + no9  r:  no  + n9o 

• Ozone  reactions: 

O + 02  + M 

O +°3 
H + 03 
NO  + Oa 


= °3  + M 

- °2  + °2 
= OH  + 02 

= N02  + °2 
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• Methane  Oxidation: 


CH4  + H 

* 

CH3  + H2 

CH4  + OH 

CHg  + H20 

ch4  + 0 

CtL  + OH 

CH4  + °2 

t; 

CH3  ♦ H02 

CH3  + 0 

t; 

CH  0 + H 

CH3  + °2 

— ► 

CH  0 + OH 

CH20  + H 

HCO  + H 

CH20  + OH 

HCO  + H20 

ch2o  + 0 

HCO  + OH 

HCO  + O 

CO  + OH 

HCO  + OH 

- 

CO  + h2o 

HCO  + H 

co  + h2 

The  above  76  reactions  (forward  and  reverse),  excluding  the  methane  oxidation  re- 
actions, can  be  taken  to  describe  the  fuel-lean  chemistry  in  the  engine  exhaust  nozzle 
and  plume  up  to  the  point  where  mixing  between  the  exhaust  plume  and  the  surrounding 
atmosphere  becomes  important.  A necessary  condition  for  their  sufficiency  is  that 
the  hydrocarbon  emissions  be  very  low  so  that  carbon  monoxide  oxidation  is  dominant. 
For  local  conditions  where  this  is  not  the  case,  the  above  well-known  methane  oxid- 
ation scheme  is  convenient  to  represent  HC  oxidation.  Since  the  number  of  possible 
species  and  reactions  is  exceedingly  large,  the  assembly  of  a complete  mechanism  with 
rate  constants  to  represent  higher  HC  (e.g.  , CgH^g)  oxidation  is  probably  not  possible 
at  this  time. 

The  reactions  and  their  respective  rates  are  summarized  in  Table  3-1.  Reverse 
reaction-rate  constants  are  calculated  from  equilibrium  constants.  References  given 
in  the  table  refer  to  the  special  reference  list  given  in  section  3.7. 


Table  3-1 


CHEMICAL  REACTIONS  IN  THE  CO  /NO  /SO  /CH  /AIR  SYSTEM 

X X X 


Reaction 

Number 

Reaction 

Reaction  Rate  Variables^) 

Ref 

(b) 

A 

n 

E 

1 

H + N20  = N2  + OH 

3.01  x 1013 

0 

10,800 

19 

2 

N + OH  = NO  + H 

4.21 

X 

1013 

0 

0 

19 

3 

0 + N20  = N2  + O2 

3.61 

X 

1013 

0 

24,000 

13 

4 

N2  + N02  — NO  + N2O 

1.41 

X 

1014 

0 

83,000 

13 

5 

H + N02  = NO  + OH 

7.25 

X 

1014 

0 

1,930 

6 

6 

H02  + NO  = OH  + N02 

6.00 

X 

1011 

0 

0 

18 

7 

0 + N2  = NO  + N 

1.44 

X 

1014 

0 

75,580 

13 

8 

NO  + 0 = N + O2 

4. 10 

X 

109 

0 

38,340 

13 

9 

NO  + 0 + M = N02  + M 

1.05 

X 

1015 

0 

-1,870 

3 

10 

NO2  + 0 = NO  + O2 

1.0  x 1013 

0 

600 

3 

11 

NO  + NO  = N20  + 0 

7.05 

X 

ion 

0 

65,000 

13 

12 

M + N2O  = N2  + O + M 

6.30 

X 

io14 

0 

56,800 

13 

13 

CO  + OH  = C02  + H 

4.  20 

X 

1011 

0 

1,080 

1 

14 

H + 02  = OH  + 0 

2.20 

X 

1014 

0 

16,800 

7 

15 

OH  + OH  = H20  + 0 

6.30 

X 

I—1 

0 

to 

0 

1,100 

7 

16 

0 + H2  = H + OH 

1.80 

X 

1010 

1 

8,900 

7 

17 

H2  + OH  = H20  + H 

2.20 

X 

1013 

0 

5,150 

7 

18 

M + H2  = H + H + M 

2.  20 

X 

1014 

0 

96,000 

7 

19 

M + H + OH  = H20  + M 

1.40 

X 

1023 

-2 

0 

7 

20 

M + 02  = O + 0 + M 

2.50 

X 

1018 

-1 

118,500 

6 

21 

M + H + 02  = HO2  + M 

1.50 

X 

1015 

0 

-1,000 

7 

22 

H + H02  = OH  + OH 

2.50 

X 

1014 

0 

1,900 

7 

23 

OH  + HO2  = H20  + 02 

1.00 

X 

1043 

0 

1,000 

7 

24 

0 + ho2  = OH  + o2 

5.00 

X 

1013 

0 

1,000 

7 

25 

H + H02  = H2  + O2 

2.50 

X 

1013 

0 

700 

7 

(a)  Reaction  rate  constant  k = ATn  exp  (E/RT).  Units  are  moles,  cm,  sec, 
°K,  cal. 

(b)  References  in  this  table  appear  as  a special  reference  list  in  Section  3.7. 
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Table  3-1  (Cont.) 


Reaction 

Number 

Reaction 

Reaction  Rate  Variables 

Ref 

A 

n 

E 

26 

H + HO£  = H2O  + 0 

1.00  x 1013 

0 

1,000 

7 

27 

N02  + OH  = HNO3 

2.90  x 1012 

0 

0 

8 

28 

H2O2  + OH  = H2O  + HO2 

1.00  x 1013 

0 

1,800 

7 

29 

OH  + OH  + M = H202  + M 

7.1  x 1014 

0 

-5,100 

7 

30 

HO2  + H02  = H2O2  + O2 

6.5  x 1013 

0 

0 

7 

31 

S02  + O + M = SO3  + M 

4.5  x 1014 

0 

0 

12 

32 

SO3  + H = S02  + OH 

6.5  x 1014 

0 

10,800 

12 

33 

SO3  + 0 = SO2  + O2 

6.5  x 1014 

0 

10,800 

12 

34 

C H4  + H = C H3  + H2 

5. 1 x 1043 

0 

12,900 

9 

35 

ch4  + OH  = ch3  + h2o 

7.9  x 1013 

0 

5,800 

9 

36 

CH4  + 0 = CH3  + OH 

2.0  x IO13 

0 

9,200 

15 

37 

CH4  + o2  = ch3  + ho2 

1.0  x 1013 

0 

63,000 

15 

38 

CH3  + 0 = CH20  + H 

1.9  x 1013 

0 

0 

10 

39 

CH3  + 02  = CH20  + OH 

7.5  x 1010 

0 

0 

10 

40 

CH20  + H = HCO  + H2 

1.0  x 1013 

0 

2,000 

10 

41 

CH20  + OH  = HCO  + H20 

7.0  x 1010 

0.7 

1 , 000 

15 

42 

CH20  + O = HCO  + OH 

4.0  x 1011 

0.6 

4,000 

10 

43 

HCO  + O = CO  + OH 

1.8  x 1011 

0.5 

0 

10 

44 

HCO  + OH  = CO  + H20 

1.1  x 1011 

0.5 

0 

10 

45 

HCO  + H = CO  + H2 

1.5  x 1012 

0.5 

0 

10 

46 

0 + O2  + M = O3  + M 

7.2  x 1012 

0 

-1,050 

17 

47 

0 + O3  = O2  + O2 

1.2  x 1013 

0 

4,800 

17 

48 

H + O3  = OH  + O2 

1.6  x 1013 

0 

0 

6 

49 

NO  + O3  = NO2  + O2 

6.7  x 1011 

0 

2,450 

6 

50 

H + H202  = H2  + H02 

1.7  x 1012 

0 

3,800 

7 
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It  turns  out  that  many  of  these  reactions  are  not  important  in  the  nozzle  and  wake 
regimes.  A condensed  list  of  essential  reactions  is  given  in  a later  section. 

3.4.3  Analysis 

The  analysis  of  the  effect  of  the  chemical  reactions  in  the  engine  exhaust  nozzle  and 
exhaust  plume  was  carried  out  using  the  NASA  Lewis  general  kinetics  program 
(Ref3-14)and  other  Lockheed  programs.  The  NASA  program  is  capable  of  treating 
reacting  flows  with  the  fluid  mechanical  equations  coupled  to  the  chemical  rate 
equations. 

The  GE-4  engine  operating  at  maximum  power  at  65,000-ft  altitude  and  flight  Mach 
number  2.7  was  taken  as  the  most  conservative  case  for  evaluating  exhaust  reactions. 
This  operating  condition  gives  the  highest  engine  temperatures  and  hence  the  greatest 
possibility  of  significant  chemical  reaction  in  the  wake. 

• The  engine  exhaust  nozzle  conditions  as  determined  from  engine  manufactur- 
er's data  are  as  follows: 

- Nozzle  throat  radius  57.9  cm 

- Throat  static  temperature  1730° 

— Throat  static  pressure  0.839  atm 

• Equivalent  fuel  C0H10,  0 = 0.715 

O lo 

• The  nozzle  shape  is  approximately  conical  with  the  radius  given  by 

r - r*  = 0 . 26  x 

where  r is  the  nozzle  radius,  x the  length,  and  * indicates  the  throat 

• The  nozzle  is  designed  to  expand  the  flow  to  ambient  pressure  - 0.0562  atm 
at  65,000  ft. 

With  the  pressure-matched  condition  at  the  nozzle  exit  plane,  the  exhaust  flow  continues 
at  constant  temperature  and  pressure  until  the  shear  layer  generated  at  the  exhaust  jet/ 
ambient-atmosphere  interface  closes  in  and  mixes  with  the  exhaust  jet. 
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NOZZLE  EXIT  PLANE 


SHEAR  LAYER  - 


Figure  3-3  Near  Jet  Exhaust 

After  leaving  the  nozzle,  the  exhaust  forms  a turbulent  jet  with  characteristics 
illustrated  by  Fig.  3-3.  The  boundaries  of  this  jet  grow  as  free  turbulent  shear  layers 
between  the  jet  core  and  the  supersonic  freestream.  These  shear  layers  grow  as  the 
jet  proceeds  downstream  until  they  meet  at  the  center,  after  which  the  exhaust  becomes 
a fully  developed  turbulent  jet.  In  the  region  between  this  meeting  point  and  the  exit 
plane  of  the  exhaust  nozzle  there  exists  a conical  shaped  region  in  the  center  of  the 
jet  which  is  unaffected  by  the  external  stream,  and  assuming  that  the  nozzle 
expansion  is  matched  to  ambient  pressure,  the  flow  properties  are  constant  in  this 
region.  Thus,  the  exhaust  gas  temperature  on  the  jet  centerline  will  be  constant, 
equal  to  its  value  at  the  nozzle  exit,  downstream  to  the  point  where  the  shear  layers 
meet. 

The  jet  exhaust  forms  a compressible  turbulent  jet  flow  in  a co-flowing  external  stream. 
For  the  purposes  of  conservative  analysis,  it  is  necessary  to  determine  the  longest 
possible  initial  constant  temperature  core.  For  this  case,  compressibility  effects 
which  tend  to  shorten  this  initial  region  are  not  included.  The  conditions  at  the  engine 
nozzle  exit  plane  were  calculated  from  published  data  for  the  GE-4  engine  at  power 
setting  1,  that  is, full  power  with  afterburner.  The  flow  was  isentropically  expanded 
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with  y = 1.28  in  the  nozzle  to  the  exit  plane  where  the  static  pressure  was  matched 
with  ambient  for  65, 000 -ft  altitude. 

Using  the  theory  of  Abramovich  (Ref.  3-15)  for  these  conditions  and  an  aircraft 
Mach  number  of  2. 7,  the  length  of  the  initial  core  region  of  constant  temperature  was 
found  to  be  22.15  m.  Experiments  pertaining  to  this  situation  were  carried  out  by 
Forstall  and  Shapiro  (Ref.  3-16)  and  an  expression  was  derived  for  the  length  of  this 
core  by  curve  fitting  data.  Use  of  this  expression  yields  a value  of  22. 0 m for  the 
length  of  the  constant  temperature  core. 

Therefore,  for  the  conditions  mentioned  above-  full  power  operation  of  the  GE-4  with 
afterburner  at  Mach  2. 7 at  65,  000  ft— the  centerline  jet  temperature  wil]  be  954°K 
and  this  value  will  be  constant  for,  at  most,  22  m downstream  of  the  exhaust  plane. 

The  complete  chemical  kinetic  model  analysis  then  involves  two  steps.  First,  the 
model  is  run  subject  to  the  constraint  of  the  conical  nozzle  flow,  and  second,  the 
calculations  are  continued  at  exit  plane  temperature  and  pressure,  out  to  the  end  of 
the  unmixed  region.  Possible  reactions  in  the  mixing  zone  and  beyond  are  considered 

in  later  sections. 

The  initial  composition  for  the  nozzle  expansion  was  taken  to  be  that  determined  by 
chemical  equilibrium  at  the  nozzle  throat,  with  the  exception  of  CO,  NO,  and  NOl? 
(equilibrium  calculations  were  performed  using  a Lockheed  Thermo-chemical  equi- 
librium program).  The  species  CO  and  NO  are  well  known  to  be  out  of  equilibrium  in 
jet  exhaust  because  of  chemical  kinetic  freezing  effects  in  the  combustor  region.  CO 
and  NO  concentrations  have  been  estimated  by  a review  of  available  aircraft  exhaust 
emission  data.  NO concentration  has  been  arbitrarily  given  a large  value  — 10  percent 
of  the  NO  concentration  — for  the  purpose  of  illustrating  the  difficulty  of  preserving  N02* 
Actual  N09  levels  in  jet  exhausts  are  highly  uncertain  at  present. 

Initial  throat  concentrations  of  all  species  considered  in  the  first  runs  are  given  in 
Table  3-2.  Hydrocarbons  and  H209  were  not  taken  into  consideration  for  this  particular 
case. 
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Table  3-2 


INITIAL  THROAT  CONCENTRATIONS 


Species 

Concentration 

(moles/cm3) 

Mole  Fraction 

H 

4.47  E-12 

7.57  E-07 

n2° 

4.79  E-13 

8.11  E-08 

N2 

4.41  E-06 

7.47  E-01 

OH 

1.97  E-09 

3.33  E-04 

N 

5.56  E-17 

9.41  E-12 

NO 

1.77  E-09 

3.00  E-04 

O 

9.25  E-ll 

1.56  E-05 

°2 

3.30  E-07 

5.59  E-02 

N°2 

1.77  E-10 

3.00  E-05 

H°2 

8.09  E-13 

1.36  E-07 

CO 

1. 77  E-08 

3.00  E-03 

C°2 

5.37  E-07 

9.10  E-02 

h2° 

6.04  E-07 

1.02  E-01 

H2 

7.58  E-ll 

1.28  E-05 

KNO 

3.80  E-17 

6.43  E-12 

5.91  E-09 

1.00  E-04 

3.4.4  Discussion  of  Kinetic  Results  Without  Excess  Hydrocarbons 

Under  these  starting  conditions,  the  model  calculations  yielded  the  concentration  pro- 
files shown  in  Figure  3-4.  Excluded  from  Figure  3-4  are  the  trace  species  SO  , N, 

O 

and  HNO  , whose  concentration  levels  remain  at  all  times  less  than  1 ppb,  and  the 

O 

major  species  CO^,  O^,  N?,  and  H90.  Concentrations  calculated  from  thermodynamic 
equilibrium  considerations  are  indicated  on  the  far  right  of  the  figure  for  comparison 

The  results  in  Figure  3-4  show  that  a nonequilibrium,  but  only  slowly  changing, 
chemical  balance  exists  in  the  exhaust  plume.  Evidently,  at  the  low  density  cor- 
responding to  65,000-ft  altitude,  recombination  is  very  slow.  All  of  the  active  in- 
organic radical  and  atomic  species  are  found  to  exist  in  super-equilibrium  concentrations. 
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Figure  3-4  Centerline  Concentration  Profiles,  GE-4  Secondary 
Nozzle  and  Exhaust  (Zero  CH^) 
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EQUILIBRIUM 


The  origin  and  maintenance  of  these  radicals  is  a consequence  of  carbon  monoxide 
oxidation.  At  throat  temperature,  CO  + OH  — C02  + H is  essentially  a branching 
step  because  the  H atoms  produced  react  via  H + 09  — O + OH,  followed  by 
0 + H20  — OH  + OH  and  H + H20  — H2  + OH.  This  chain  branching  effect  is  con- 
firmed by  the  numerical  results  in  Table  3-3.  There  are  given  reaction  rates  and 
net  radical  species  production  for  the  six  reactions  mainly  governing  H,  O,  and  OH 
levels.  It  can  be  seen  that  while  the  CO  oxidation  reaction  rapidly  consumes  OH,  the 
net  effect  of  all  six  reactions  is  a rapid  generation  of  additional  OH  as  well  as  O and 
H radicals. 


Table  3-3 

RADICAL  FORMATION  AT  STATION  10-cm 
DOWNSTREAM  OF  NOZZLE  THROAT 


Reaction 

Net  Reaction  Rate,  r -► 
(moles/cm^-sec)  n 

5 

H + N02 

^ NO  + OH 

0.46  E-6 

13 

CO  + OH 

= co2  + h 

7.74  E-6 

14 

H+°2 

^ OH  + O 

5.03  E-6 

15 

OH  + OH 

^ h2o  + O 

-2.64  E-6 

16 

° + H2 

■z.  H + OH 

-0.06  E-6 

17 

H9  + OH 

^ h2o+h 

-1.49  E-6 

Net  Effect  on  Radicals 


- net  dt^  = -r5  + r13  ' r14  + r16  + r17  = °'7  E-6* 

“et  ^t21"  r14+r15'r16  = 2-5E-6» 

net  ijOHl  = rg  - r13  + ru  -2rlg  + rlg  - r1?  = 4.5  E-6* 

*(moles/ cm^-sec) 

The  initial  fast  reaction  only  serves  to  demonstrate  that  the  assumption  that  the  radicals 
were  in  equilibrium  with  major  species  in  the  throat  is  unrealistic.  It  is,  of  course, 
generally  not  possible  for  one  species  (e.g.  , H)  to  be  in  equilibrium  and  another 
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(e.g.  , CO)  to  be  out  of  equilibrium,  when  there  is  a fast  reaction  (e.g.  , CO  + OH 
CO?  + H)  available  to  exchange  them. 

In  the  cooler  exhaust  jet,  CO  oxidation  is  still  a chain  mechanism  since  now  CO  + OH 
CO  + H is  followed  mainly  by  H + O + M — HO  + M and  then  H + HO  — OH  + OH. 

u Z Z Z 

This  notion  that  the  high  radical  concentrations  are  a result  of  carbon  monoxide  oxi- 
dation as  the  driving  reaction  is  confirmed  by  computations  in  which  the  CO  oxidation 
rate  is  arbitrarily  reduced  to  negligible  values.  The  output  then  shows  about  an  order 
of  magnitude  lower  concentration  of  O,  H,  and  OH  radicals  in  the  nozzle. 

It  is  encouraging  to  note  that  a recent  optical  measurement  (Ref.  3-13)  of  OH  concen- 
tration in  a relevant  afterburning  situation  under  altitude  condition  gives  an  OH  con- 
centration in  good  agreement  with  the  super-equilibrium  value  computed  here. 

Understanding  of  the  details  of  the  various  chemical  balances  in  the  exhaust  jet  is 
aided  by  Table  3-4.  There  are  listed  rate  constants,  toward  reaction  rates,  and  net 
rates  for  the  various  reactions  at  a position  2 m downstream  from  the  nozzle  exit 
plane.  Table  3-4  can  then,  in  turn,  be  used  to  identify  the  reactions  that  are  of 
importance  to  individual  species  rates.  For  example,  the  OH  balance  shown  in 
Table  3-5  is  instructive.  This  gives  all  the  important  reactions  involving  OH  as  well 
as  their  rates  and  the  net  effect.  As  can  be  seen  in  the  table,  the  aforementioned  CO 
oxidation  chain  has  the  consequence  that  the  net  disappearance  of  OH  is  more  than  an 
order  of  magnitude  less  than  its  rate  of  participation  in  CO  + OH  — C09  + H.  Hence, 
it  would  be  a good  approximation  to  assume  the  OH  concentration  constant  in  com- 
puting the  rate  of  CO  oxidation.  The  very  sluggish  nature  of  the  net  recombination 
means  that  the  radicals  will  persist  out  into  the  region  where  the  plume  mixes  with 
ambient  air. 

3.4.5  Exhaust  Chemistry  With  Significant  Hydrocarbons 

Test  data  are  somewhat  scarce  for  SST  engines,  but  existing  data  show  that  quantities 
of  unburned  hydrocarbons  (measured  as  total  C)  may  be  large  along  certain  streamlines 
with  afterburning  (Refs.  3-17  and  3-18).  To  take  some  account  of  the  effect  of  hydrocarbor 
on  the  chemistry,  reactions  34  through  45  (Table  3-1),  a methane  oxidation  scheme, 
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Table  3-4 

CHEMICAL  KINETIC  SITUATION  TWO  METERS  DOWNSTREAM 
OF  NOZZLE  EXIT  PLANE 


REACTION 

FORWARD  RATE 
CONSTANT 

k* 

FORWARD  REACTION 
RATE 

rf  =kf  7T  [A  ]Vi 
i 1 

NET  REACTION  RATE 

r*  - kr  7r  (A.)'u 

i 3 

COMMENTS  * 

1 

h + n2o 

N2  + OH 

1.01 

E+ll 

1.57  E-12 

1.57  E-12 

F 

2 

N + OH 

NO  + H 

4.21 

E+13 

3.21  E-17 

-3.63  E-17 

B 

3 

o + n2o 

— 

N2  + °2 

1.14 

E+8 

1.39  E-15 

1.39  E-15 

F 

4 

n2  + no2 

— 

NO  + N20 

1.35 

E-5 

1.94  E-26 

-7.46  E-21 

R 

5 

H + N02 

— 

NO  + OH 

2.62 

E+14 

1. 86  E-10 

1.86  E-10 

F 

6 

H02  + NO 

— 

OH  + N02 

1.25 

E+ll 

1.51  E-ll 

1.51  E-ll 

F 

7 

o + n2 

NO  + N 

6 93 

E-4 

7.69  E-20 

-3. 66  E-17 

R 

8 

NO  + O 

= 

N + 02 

6.73 

E+0 

3.30  E-19 

1.78  E-19 

B 

9 

NOz  + M 

NO  + O + 

M 

1.41 

E+l 

2 72  E-20 

-1.75  E-10 

R 

10 

NO  + 02 

o + no2 

4.20 

E+l 

3 93  E-16 

-4.01  E-12 

R 

11 

NO  + NO 

n2o  + 0 

9 01 

E-4 

5.  06  E-23 

-1  94  E-16 

R 

12 

NzO  + M 

- 

n2  + o + 

M 

6.09 

E+l 

2. 56  E-18 

-4  42  E-14 

R 

13 

CO  + OH 

- 

C°2  + H 

3.17 

E+ll 

2. 32  E-8 

1.85  E-8 

B 

14 

H + 0o 

— 

OH  + O 

3.11 

E+10 

3.26  E-7 

-2.30  E-9 

FASTEST  REACTION 

2 

BUT  VERY  NEAR  PE 

15 

OH  + OH 

= 

H20  + 0 

3.53 

E+l  2 

5.91  E-8 

-8.48  E-12 

PE 

16 

° + h2 

H + OH 

1.57 

E+ll 

8.06  E-9 

9. 74  E-10 

NEAR  PE 

17 

H2  + OH 

zt 

h2°  + h 

1.45 

E+l  2 

4.65  E-8 

5.63  E-9 

NEAR  PE 

18 

H2  + M 

= 

H + H + M 

2.22 

E-8 

3.96  E-24 

-2. 26  E-10 

R 

19 

H + OH  + M 

HO  + M 

1. 54 

E+l  7 

3. 80  E-9 

3.78  E-9 

F 

20 

°2  + M 

= 

O + O + M 

1.85 

E-12 

5.  26  E-26 

-3.46  E-12 

R 

21 

M + H + 02 

zl 

H02  + M 

2.54 

E+l  5 

1.91  E-8 

1.91  E-8 

F 

22 

H + H02 

zr 

OH  + OH 

9.175 

E+13 

1. 24  E-8 

1.24  E-8 

F 

23 

OH  + H02 

— 

h2°  + o2 

5.90 

E+l  2 

3.90  E-10 

3.90  E-10 

F 

24 

o + ho2 

OH  + 02 

2.95 

E+13 

3.12  E-9 

3.12  E-9 

F 

25 

H + H02 

H2  + °2 

1.73 

E+13 

2.34  E-9 

2.34  E-9 

F 

28 

H + H02 

H2°  + ° 

5.90 

E+l  2 

7.98  E-10 

7.98  E-10 

F 

27 

NOz  + OH 

HN°3 

2.90 

E+l  2 

1.00  E-12 

-1.62  E-17 

PE 

*F  - FORWARD  REACTION  DOMINATES. 

R - REVERSE  REACTION  DOMINATES. 

PE  - FORWARD  RATE  = REVERSE  RATE,  REACTION  IN  PARTIAL  EQUILIBRIUM. 
B - BOTH  FORWARD  AND  BACK  RATES  SIGNIFICANT  AND  NONEQUAL. 
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Table  3-5 


PREDOMINANT  CONTRIBUTIONS  TO  OH  BALANCE  TWO  METERS 
DOWNSTREAM  OF  EXHAUST  NOZZLE  EXIT  PLANE 


NO. 


22f  H + 
24f  O + 


REACTION 


H02  — OH  + OH 

ho2  — OH  + O2 


16  O + H2 

5f  H + N02 


— OH  + H 

— OH  + NO 


13 


CO  + OH 


- CO2  + H 


17 

19f 

14 

23f 


H,  + OH 
H + OH  + M 
O + OH 
H02  + OH 


NET  CHANGE  OH 


r 

r 

r 

r 

- r 

- r 

- r 

- r 

- r 


f 

22 

f 

24 

16 

f 

5 

13 

17 

f 

19 

14 
f 

23 


z:  h2o  + H 

— H20  + M 

zr  h + o2 

— H20  + O9 


= 2.48  E-8 
= 0.31  E-8 
= 0.10  E-8 
= 0.02  E-8 
= - 1.85  E-8 
= - 0.56  E-8 
= - 0.38  E-8 
= - 0.23  E-8 
= - 0.04  E-8 


COMMENT 
OH  FORMATIONS 
PRIMARY 
SECONDARY 

| SMALL  CONTRIBUTIONS 

OH  DESTRUCTIONS 
PRIMARY 

SECONDARY 

SMALL  CONTRIBUTION 


(MOLE/CM3-SEC) 

(MOLE/CM3-SEC) 

(MOLE/CM3-SEC) 

(MOLE/CM3-SEC) 

(MOLE/CM3-SEC) 

(MOLE/CM3-SEC) 

(MOLE/CM3-SEC) 

(MOLE/CM3-SEC) 

(MOLE/CM3-SEC) 


d [QHl  _ 
dt 


-0.15  E-8  (MOLE/CM3-SEC) 


d [CO] 

TO  BE  COMPARED  WITH  —7—  = -r,0 

at  lo 


-1.85  E-8  (MOLE/CM3-SEC) 
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have  been  included  at  this  point.  As  we  shall  see,  sacrifice  of  p (defined  as  p = H, 

OH,  and  O together)  in  favor  of  CH3,  0H2O,  HCO,  and  CO  is  the  net  effect  on  the 
exhaust  chemistry . The  magnitude  of  p radical  reduction  is  highly  dependent  on  the 
HC  quantities. 

Two  sets  of  concentration  profiles  are  shown  in  Figures  3-5  and  3-6.  Both  are  for 

the  same  operating  condition  as  Figure  3-4.  Starting  throat  concentrations  of  all  species 

are  taken  as  before,  except  that  the  portion  of  CH^  is  assumed  initially  at  500  ppm  in 

Figure  3-5  and  3,000  ppm  in  Figure  3-6.  In  the  first  case  most  of  the  CH^  is  oxidized 

in  the  nozzle  and  jet,  with  only  50  ppm  total  HC  (including  CHO  and  CH  O)  remaining 

£ 

at  the  exit  plane  and  subsequently  disappearing  in  the  near  wake.  Total  p is  forced 
down  to  about  one-fifth  of  its  zero  HC  level  at  the  exit  plane,  but  builds  up  somewhat 
in  the  wake  as  the  last  of  the  HC  vanishes. 

When  the  CH^  is  set  at  3,000  ppm  throat  concentration,  very  little  is  completely  oxi- 
dized; i.e. , the  wake  HC  level  is  frozen  at  3,000  ppm  with  about  half  of  it  present  as 
CH^  and  the  rest  as  CH90  and  CHO.  This  has  the  consequence  of  nearly  eliminating 
p from  the  wake,  as  shown  in  Figure  3-6. 

Although  it  may  seem  peculiar  at  first  that  the  addition  of  CH^  reduces  the  levels  of 
H,  OH,  and  O,  it  should  be  emphasized  that  this  is  a highly  nonequilibrium  situation, 
where  super-equilibrium  amounts  of  methane  are  introduced  and  in  turn  produce 
super-equilibrium  amounts  of  hydrocarbon  radicals  (CH  , CH00,  and  HCO)  at  the 
expense  of  H,  OH,  and  O.  In  fact,  Figure  3-5  shows  that  H,  OH,  and  O begin  to  return 
as  the  added  hydrocarbons  disappear.  For  higher  levels  of  added  CH^,  as  in  Fig- 
ure 3-6,  the  hydrocarbon  oxidation  is  incomplete  and  the  H,  OH,  and  O remain  bound 
in  the  frozen  hydrocarbon  radicals.  The  levels  for  the  three  cases  are  shown  plotted 
together  in  Figure  3-7  for  comparison. 

These  calculations  serve  to  point  out  the  profound  difference  of  carbon  monoxide  and 
hydrocarbon  oxidation  at  nozzle  expansion  temperatures.  The  former,  CO  + OH  — 

CO  + H,  builds  up  p through  activity  of  the  H atom,  whereas  the  hydrocarbon  oxi- 
dation cannot  sustain  production  of  p , and  freezing  occurs  as  soon  as  the  available 
radicals  are  scavenged.  The  general  conclusion  drawn  here  is  that,  if  the  HC  levels 
are  high  the  radical  H,  OH,  or  O cannot  be  present  in  significantly  large  concentrations. 
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Figure  3-5  Centerline  Concentration  Profiles,  GE-4  Secondary  Nozzle 
and  Exhaust  (500  ppm  CH^  at  Throat) 


3-34 


1000  ppm  10 


CH, 

4 

ch2o 

CHO 


NO 


1 ppm  10 


1 

H 

O 

2 

w 

o 

2 


.-6 


Figure  3-6 
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Centerline  Concentration  Profiles,  GE-4  Secondary  Nozzle 
and  Exhaust  (3,000  ppm  CH4  at  Throat) 
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Figure  3-7  Effect  of  Hydrocarbon  on  N02  in  GE-4  Secondary 
Nozzle  and  Jet  Exhaust 
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3.4.6  Radical  Partial  Equilibrium 


Another  simplifying  condition  evident  from  Table  3-4  is  that  the  bimolecular  shuffle 
reactions  14  through  17  are  in  partial  equilibrium.  That  is,  their  forward  reaction 
rates  are  approximately  equal  to  their  reverse  rates.  This  condition  is  exhibited  in 
the  table  by  the  fact  that  the  forward  rate  is  much  larger  than  the  net  rate.  Since  the 
ratio  of  forward  to  back  rate  constant  is  the  equilibrium  constant,  an  algebraic  equation 
can  be  written  relating  the  concentration  of  species  in  each  of  these  reactions,  as 
follows: 


where  K represents  the  equilibrium  constants  and  the  brackets  signify  concentrations. 
Since  only  three  of  these  equations  are  independent,  a complete  solution  is  not  ob- 
tainable, but  the  concentration  of  the  variable  species  H,  OH,H2>  and  O can  be  ex- 
pressed as  a function  of  just  one.  If  OH  is  selected  as  the  independent  quantity, 
we  have 


= [OH]  [01 
14  [H]  [02  ] 


[h2o]  [Ol 


[OH2] 


= [H]  [OH] 
16  [O]  [H2] 


[h2o]  [H] 

K17  = [H2]  [OH]  = K15K16 


[H] 


2 


[h2] 


[H]  = [OH] 


[o]  k14  [o2] 
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These  relations  are  expected  to  be  of  value  in  setting  initial  nozzle  exit  conditions  from 
test-cell  measurements  of  OH.  All  the  inorganic  radical  concentrations  must  be  known 
if  meaningful  downstream  calculations  are  to  be  made.  It  is  clear  from  the  results 
in  the  last  two  sections  that  probably  none  of  these  radicals  are  in  absolute  chemical 
equilibrium  anywhere  in  the  nozzle.  Some  experimental  measurements  are  therefore 
essential. 

3.4.7  Discussion  of  Axial  CO  Oxidation 

The  initial  CO  concentration  far  exceeds  its  local  equilibrium  value  and  therefore 
tends  to  oxidize  to  C02  in  the  nozzle  and  exhaust  jet.  Oxidation  of  CO  is  governed 
by  its  reaction  with  OH  radicals,  CO  + OH  — • C02  + H.  The  H atoms  produced  by 
this  reaction  then  enter  the  fast  bimolecular  shuffle  reactions  14  through  17  (Table  3-1) 
and  OH  is  regenerated.  These  fast  reactions,  coupled  with  the  slow  recombination,  act 
to  maintain  O,  H,  and  OH  at  super-equilibrium  concentrations. 

The  phenomenological  aspects  of  CO  oxidation  have  already  been  discussed  earlier. 
Numerical  results  show  that  a considerable  decrease  in  CO  can  occur  in  the  nozzle 
and  jet  under  low  HC  conditions.  Calculations  for  the  zero  HC  run,  presented  in  Fig- 
ure 3-4  show  that  the  proportion  of  CO  in  the  exhaust  is  decreased  from  3,000  ppm  in 
the  nozzle  throat  to  1,780  ppm  at  the  exit  and  then  is  further  reduced  to  1,150  ppm  at 
the  tip  of  the  jet  core  region. 

The  corresponding  levels  for  the  runs  with  500  ppm  CH^  initially,  where  CO  was  again 
3,000  ppm  at  the  throat,  are  2,810  ppm  at  the  exit,  and  2,180  ppm  finally.  When  the 
portion  of  CH^  was  set  at  3,000  ppm,  the  CO  level  remained  essentially  frozen 
throughout. 

The  kinetic  model  is  run  at  constant  pressure  for  the  jet  calculations  so  a chemically 
induced  temperature  change  is  computed  between  the  nozzle  exit  and  the  jet  cone  tip. 

In  the  zero  HC  case,  which  had  the  largest  amount  of  CO  oxidized,  the  magnitude  of 
this  change  was  still  only  5°  K (in  950°  K) . This  tends  to  confirm  our  intuitively  held 
notion  that  the  exhaust  species  chemistry  can  have  no  meaningful  effect  on  the  fluid 
mechanic  properties. 
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3.4.8  Integrated  CO  Oxidation  in  Exhaust  Jet  Core 


The  above  results  without  hydrocarbons  show  oxidation  of  35%  of  the  available  CO  occurs 
between  the  nozzle  exit  plane  and  the  tip  of  the  constant  properties  region.  However, 
this  result  is  for  the  centerline  of  the  wake;  streamlines  nearer  the  jet  boundary  en- 
counter the  reaction  quenching  inward  propagating  shear  layer  earlier,  so  the  total  in- 
tegrated CO  oxidation  is  not  so  great. 

If  the  wake  zone  of  constant  conditions  is  assumed  to  be  exactly  conical,  the  total 
reduction  can  be  estimated  in  a simple  manner . Let  L be  the  length  of  the  undisturbed 
exhaust  cone;  x be  the  streamline  distance;  CO  (x)  be  the  axial  CO  profile  calculated 
assuming  constant  properties;  r be  the  local  radius  for  nozzle  centerline;  and  subscript 
E refer  to  the  exit  plane.  Clearly  the  CO  contained  in  the  annulus  of  area  dA  = 2n  rdr 
is  frozen  at  the  position 


Therefore,  the  portion  of  the  total  CO  remaining 


T 


Total  CO 


contributed  by  this  element 


Total  CO  in  exit  plane 


is 


_ [CO(x)l  2t r rdr  = [COQg)]  2_r_  d 

" lcolEAE  [c°1E  4 

E 


and  with  r = r£  (1  - x/L)  and  dr  = -j-  dx 
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so  that  the  total  remaining  is  given  by 


0 


0 


A numerical  integration  of  this  equation  using  the  [CO(x)]  profile  data  presented  in 
Figure 3-4  yields  the  value  T =0.87  corresponding  to  a 13%  total  CO  oxidation  beyond 
the  nozzle  exit. 

The  extent  of  CO  oxidation  in  the  far  wake  beyond  the  jet  core  depends  almost  entirely  on 
the  rate  at  which  the  exhaust  becomes  diluted,  the  rate  falling  off  as  the  square  of  the 
dilution.  The  significance  of  the  CO  oxidation  is  not  the  final  wake  concentration  of 
CO  or  the  total  CO  added  to  the  stratosphere.  The  estimated  total  impact  of  the  CO 
added  by  an  SST  fleet  is  to  increase  ambient  stratospheric  levels  by  approximately 
0.2  ppb/30  ppb  =0.7%  (Ref.  3-19). 

Carbon  monoxide  is  not  a "pollutant"  in  this  sense.  The  significance  of  CO  oxidation 
is  the  resultant  production  of  hydrogen  atoms . These  contribute  substantially  to 
reactions  which  tend  to  prevent  oxidation  of  pollutants  like  NO  and  SO^  in  the  exhaust 
jet.  This  point  will  be  discussed  more  fully  in  later  paragraphs. 

3.4.9  Discussion  of  NO  Results 


Oxides  of  nitrogon  are  the  species  of  major  interest  with  regard  to  environmental 
impact.  The  portion  of  NO  , as  N09  in  the  hot  exhaust  is  of  importance  because  of  the 

X z 

possibility  of  nitric  acid  formation  to  provide  an  NO  "sink"  in  the  wake, 

X 


Experimental  evidence  (Refs.  3-20,  3-21)  indicates  that  the  afterburning  region  in 
augmented  turbojet  engines  produces  very  little  NO  so  that  NO  levels  in  the  exhaust 
are  associated  with  the  NO  production  in  the  main  combustors.  Typical  combustor 


x 


N02  + OH  — HNOs 
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NO  exit  levels  range  from  100  to  500  ppm  (Refs.  3-22 , 3-23).  The  initial  NO 
concentration  is  taken  here  as  300  ppm. 

Emission  levels  ofNO^  are  more  uncertain  than  the  NO  levels.  Both  Fletcher  (Ref.  3-23) 
and  Caretto  (Ref.  3-24)  indicate  that  at  full  power  the  oxides  of  nitrogen  will  be  essen- 
tially all  nitric  oxide,  NO.  However,  recent  test  data  (Ref.  3-25)  seem  to  indicate  that 
significant  N09  levels  exist  at  all  engine  power  levels. 

In  the  measurement  of  NO^  levels,  one  method  is  to  use  a chemilumescent  detector 
with  a converter  to  convert  N09  to  NO  so  that  the  detector  measures  total  NO  . NO 
is  then  determined  by  the  difference  between  NO  and  NO.  Mellor  (Ref.  3-26)  has  re- 

X 

cently  reported  problems  in  using  the  NOg-NO  converters  and  has  attributed  lack  of 
reproducibility  of  results  to  catalytic  reactions  in  the  converter.  This  illustrates  the 
point  that  the  presence  of  catalytic  walls,  such  as  sampling  probes  in  the  flow,  can 
promote  reactions  between  NO  and  oxygen  species  which  might  not  take  place  in  the 
gas  phase.  Recent  test  data  (Ref.  3-25)  utilized  an  electrochemical  device  to  measure 
N09  and  NO  . 

Z X 

Thus,  at  the  present  time  NO^  emission  levels  are  quite  uncertain.  For  this  investi- 
gations we  have  taken  30  ppm  for  the  NO  concentration  at  the  nozzle  throat. 

Thermodynamic  equilibrium  considerations  give  an  NO^/NO  ratio  of  0.01  at  exit  plane 
conditions.  The  results  of  the  kinetic  calculations  without  HC,  as  shown  in  Figure  3-4, 
reveal  that  conversion  of  NO  to  N09  actually  does  not  occur.  In  fact,  any  NO^  assumed 

present  at  the  throat  is  quickly  destroyed.  At  the  exit  plane  the  concentration  has  fal- 

-5 

len  to  4 ppb,  corresponding  to  an  N02/NO  ratio  of  only  10 

NO  is  kept  far  below  its  equilibrium  level  by  fast  reaction  such  as 

Li 

91  NO  + O + M — N09  + M 

6f  NO  + H09  — N02  + OH 

5f  H + N02  — NO  + OH 

10r  O + N02  — NO  + 02 


3-41 


and  the  fact  that  the  H and  the  O levels  are  maintained  far  above  their  equilibrium 
concentrations.  The  second  and  fourth  reactions  above  become  unimportant  as  the 
temperature  drops  and  are  already  much  slower  than  the  other  two  reactions  at  the 
exit  plane. 

The  contributions  of  the  NO^  balance  at  the  2-m  position  can  be  tabulated  as  follows. 
Net  Change  NO  , No  Excess  HC: 

-10  ^ 

= 1.86  x 10  (moles/cm  -sec) 

-10  , 8 

= 0.15  x 10  (moles/cm  -sec) 

-10  .8 

= - 1.97  x 10  (moles/cm  -sec) 

-10  8 

= - 0.04  x 10  (moles/cm  -sec) 


-10  Q 

= Z r < 0.01  x 10  (moles/cm  -sec) 


- r. 


- r 


10 


d [ NOr 
dt 


Thus , in  the  absence  of  HC  the  N09  level  is  largely  controlled  by  a steady  state 
between  its  formation  by  reaction  with  O atoms  and  destruction  by  H atoms.  Given 
such  a balance  between  reactions  9 and  5,  an  expression  for  the  steady-state  N02/NO 
ratio  is  obtained  as  follows: 


O 


d [N02] 
dt 


k9  [NO]  [O]  [ M]  - k5  [H]  [N021 
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Therefore 


'SS 


kq  [Mj  [Oj 


9 


kg  [H] 


The  O/H  ratio  is  thus  the  controlling  quantity. 

Substituting  the  zero  HC  exit  plane  concentrations  of  M,  O,  and  H,  and  kr  and  kn 

5 9 

evaluated  at  the  exit  plane  temperature,  in  the  expression  yields 


This  result  compares  favorably  with  the  results  of  the  full  computer  solution,  i.e.  , 


For  minor  species  resulting  from  rapid  formation  and  even  faster  destruction  proc- 
esses, a simple  local  steady  state  approximation  can  often  furnish  a useful  order-of- 
magnitude  estimate  of  the  ratio  between  the  oxidized  and  reduced  forms  of  the  species 
in  question. 

When  sufficient  hydrocarbons  are  present,  the  steady-state  NO  position  is  not 
reached.  The  expression  above  is  still  a good  approximation  in  the  case  of  500  ppm 
CH  , but  with  3,000  ppm  CH  the  H is  driven  so  low  that  NO  is  no  longer  destroyed 
and  ( [ NO  ] / [NO]  ) = 0.05.  Of  course,  any  formation  would  also  be  slowed  since 
O is  driven  equally  as  low  as  H. 


(no2) 

^N0Wit,  steady  state 


0.62  x 10 


-5 


(no2) 

^N^exit,  finite  rate 


1.1  x 10 


-5 
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3.4.10  H,  O,  and  OH  Recombination  in  Cooling  Wake 


An  assessment  of  predominant  radical  recombination  paths,  as  air  is  entrained  and  the 
wake  cools,  can  be  based  on  calculations  such  as  presented  in  Table  3-6.  There  are 
listed  rate  constants  and  initial  reaction  rates  for  recombination  and  contributing  bi- 
molecular  reactions,  assuming  the  exit  plane  species  are  cooled  instantaneously  to 
220°K.  A ten-fold  air  dilution  has  also  been  assumed  for  these  sample  calculations. 
The  original  exit  plane  concentrations  are  those  corresponding  to  Figure  3-4. 

It  is  evident  from  the  tabulated  results  that  each  radical  has  one  prevailing  recombina- 
tion reaction  that  is  at  least  an  order  of  magnitude  faster  than  its  competitors.  For 
H atoms  this  reaction  is  (4)*  H + O + M — HO  + M,  for  OH  radicals  it  is  reaction 

Z Z 

(6)  OH  + OH  + M—H  O + M,  and  for  O atoms  reaction  (9)  O + O + M — O + M 

Z Z Z o 

dominates.  , 

The  half  time  for  O atom  recombination  via  (9)  at  220  °K  is  given  by: 


ln2 

1/2  ' k9[02][M] 


= 4 x 10 


sec  at  20  km 


H and  OH  recombination  are  similarly  fast. 

Reaction  (7)  SO  + O + M — SO  + M and  (8)  NO  + O + M — NO  + M are  slow  com- 

Z O Z 

pared  to  loss  of  O atoms  in  (9).  Keeping  in  mind  that  (9)  becomes  relatively  even  faster 
with  more  air  dilution,  it  is  apparent  that  the  extent  of  SO.-,  and  NO  oxidation  in  the  cool- 
ing wake  by  (7)  and  (8)  will  be  small. 


*Retum  number  in  this  discussion  refer  to  those  used  in  Table  3-6. 
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Table  3-6 

REACTION  RATES  FOR  SUDDEN  COOLING  OF  EXHAUST  TO  220  K 


Reaction 

k 

(mole,  cm,  sec) 

r "*  3 

(male/cm  -sec) 

(1) 

H + OH 

M 

— 

H£0  + M 

2.  9 E18 

3.  2 E-8 

(2) 

H + H 

+ 

M 

- 

H2  + u 

3.  0 E15 

4.2  E-ll 

(3) 

O + O 

+ 

M 

— 

°2  + M 

1.  7 E15 

2.  3 E-ll 

(4) 

H + °2 

+ 

M 

- 

H02  + M 

1.5  E16 

2.2  E-6 

Prevailing  H sink 

(5) 

H + 0 

+ 

M 

— 

OH  + M 

5.  0 E15 

7.3  E-ll 

(6) 

OH  + OH 

+ 

M 

— 

H2°2  + M 

1.  1 E20 

8.  0 E-7 

Prevailing  OH  sink 

(7) 

S02  + o 

+ 

M 

— 

S03  + M 

4.  0 E15 

2.  6 E-ll 

(8) 

NO  + O 

+ 

M 

- 

N02  + M 

7.  6 E16 

1.  6 E-9 

(9) 

0 +02 

+ 

M 

— 

°3  + M 

9.  0 E14 

1.2  E-7 

Prevailing  O sink 

(10) 

OH 

+ 

0 

— 

°2  + M 

1.  5 E13 

5.  0 E-8 

(11) 

OH  + OH 

— ■ 

h2°  + o 

5.  1 Ell 

1.2  E-9 

3.4.11  NO  in  the  Near  Wake 
x 

Another  possible  oxidation  process  for  NO  is  the  Bodenstein  mechanism 


NO  + NO  + 02  — N02  + N02 


This  has  the  rate  equation 


d (NO) 
dt 


2 kb  (NO)2  (02) 


and  hence  a half  time  for  NO  disappearance  of 


1 

Tl/2  2 kb  <NO)0  <02) 
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This  relationship  between  an<^  (NO)^  *S  s^own  Pitted  in  Figure  3-8  for  some 

altitudes  of  interest.  The  rate  constant  is  k,  = 2.4  x 10^  exp  (+  1046/RT) 

6 —2  —1  ^ 
cm  -mole  -sec  (Ref.  3-27).  Figure  3-8  shows  that  the  Bodenstein  oxidation  is 

insignificant  in  the  wake. 


Upon  mixing  of  the  exhaust  with  ambient  air,  ozone  is  introduced  and  the  subsequent 
effect  of  the  reaction 


k 

NO  + 03  - N02  + 02 

must  be  considered  from  the  standpoint  of  NO  oxidation. 
The  NO  depletion  rate  of  Eq.  (3. 25)  is  given  by 


d (NO) 
dt 


-ka  (NO)  (03) 


Integrating  the  rate, 


or 


NO 


I 

<NO>o 


d (NO) 
(NO) 


/ - ka  <°3>  dt 
0 


(3.25) 


Now  if  the  exhaust  mixing  is  extensive  enough  so  that  (O  ) » (NO)  , the  ozone  concentra- 

O 

tion  will  be  constant,  and  the  half-time  is 


In  2 

Tl/2  - ka  (03> 
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Figure  3-8  Half  Time  for  NO  Oxidation  by  NO  + NO  + 09  — ■ NO  + NO 

& La  Li 
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For  20-km  altitude  conditions,  Johnston  (Ref.  3-27)  gives  k 

-1  -1  12  -3  a 

molecule  -sec  and  (O^)  ~ 10  molecule-cm 


10 


-14.40 


3 

cm  - 


Therefore 


rl/2  = sec  ~ 2 m^n 

It  is  also  instructive  to  compute  the  half-time  in  the  event  that  (0o)  = (NO)  with  no 

3 o o 

diffusional  replenishment  of  03  as  the  reaction  proceeds . 

In  this  case,  the  half-time  expression  is 

ti/2  ■ v^r0 

which,  when  evaluated  for  the  same  initial  conditions  as  before,  yields 

Tjyg  = 252  sec  ~ 4 min 

Thus,  it  is  clear  that  the  oxidation  of  NO  by  ozone  is  rapid  enough  to  be 
significant  in  the  wake  regime.  In  the  near  wake  the  rate  will  be  limited  by  the  ozone 
supply.  The  O^/NO  approach  to  steady  state  in  the  irradiated  environment  is  treated 
more  fully  in  Section  4. 

3.4.12  Oxidation  of  NO  to  UNO  in  Wake 

Li  O 

In  the  mixing  region,  when  a sufficiently  low  temperature  is  reached  for  nitric  acid  to 
be  stable,  NO  can  be  oxidized  by  the  reaction 

k 

N02  + OH  i HNOs  (3.26a) 
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Although  it  would  appear  from  foregoing  calculations  that  the  hot  exhaust  will  contain 

very  little  NO  as  NO  , this  is  unconfirmed  by  test  data  and,  therefore,  reaction  (3.26a) 

must  be  considered.  However,  even  assuming  all  NO  is  emitted  as  NO  , the  follow- 

X z 

ing  argument  will  demonstrate  that  NO  removal  by  reaction  (3.26a)  is  still  not  of  an  im- 

X 

portant  magnitude. 

As  discussed  in  Section  3.4.10  the  reaction 

k, 

OH  + OH  + M -2  H O + M (3.26b) 

& Li 

is  fast.  If  it  is  sufficiently  so,  OH  will  be  totally  removed  before  reaction  (a)  can 
proceed  to  a significant  degree.  This  is  confirmed  by  the  following  analysis  of  the 
reaction  system  (3.26a)  plus  (3.26b). 

The  rate  of  OH  disappearance  from  reactions  (3.26a)  and  (3.26b)  is  given  by: 

= -2  \ [OH]2  [M]  - ka  [N02]  [OH] 


approximating  [NO^]  by  an  average  value  [NO^]  , (justified  in  the  light  of  results) 
an  integration  produces  the  relation 


[OH] 


ka  [OH]0  [NO,] 
+ 21^  [OH]  q [M] 


-ka  [no2]  t 


1 - exp  ( - k& 


where  subscript  zero  denotes  initial  conditions. 


This  expression,  when  substituted  for  OH  in  the  NO^  rate  equation 


d [N02] 
dt 


- \ [no2]  [oh] 
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yields 


1 d [N02] 

[no2]  dt 


- ka  [OH]0  [NO,]  exp  (-  ka  LN02]  t) 
ka  LN02]  + 21^  [OH]0  [M]  [l  - exp  (-  ka  [N02]  t)] 


Another  integration  then  gives: 


[NOJ 


[NO 


2'  0 


_ OL 


a 


+ Tift  t ' 6XP  (’  ka  [N°2]  ‘). 
2 


where 


a = 


kfetM] 


Taking  the  limit  as  t — > gives  the  final  quantity  of  N02  in  the  wake  as 


tN°2k 


= [NO,] 


Q! 


2J0  ] [OH] ( 

a + 


[nct„] 


This  expression  may  be  evaluated  by  assuming  [NO.J  = [N02]q  as^a  first  approxima- 
tion, calculating  [NO  ] on  this  basis,  and  then  taking  [NO*  ] = — — ^ — 

L Z 00  Z~*  Z 

and  repeating  the  calculation. 


Results  of  the  calculations  are  presented  in  Figure  3-9  in  terms  of  the  fraction  of  N02 

oxidized  to  HNO  , 1 - [NO  ] /[NO  ] , as  a function  of  the  initial  ratio  [NO  ] / 

o Z oo  Z~*  Q Z O 


[°H]0  for  altitudes  of  interest. 
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FRACTION  N00  CONVERTED  TO  HNO 


[N02]0/[OH]0 


Figure  3-9  Fraction  of  N02  Oxidized  to  HN03  in  Infinite  Time 
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In  generating  Figure  3-9,  the  rate  constant  was  evaluated  using  the  parameters 
given  in  Table  3-1.  For  k the  recent  expression  cited  by  Crutzen  (Ref.  3-28)  was 

cl 

used: 


k 

a 


6. 3 x 1012  exp 


p?) 


2.4  x 10 


13 


M 


13 


2.4  x 10~  [M]  + 1013/ ( 1 + 


(i  + 5m.) 


Here  his  expression  has  been  converted  to  mole,  cm  , sec  units. 


By  using  these  rate  data,  the  three  solid  lines  on  Figure  3-9  can  be  computed.  The 

maximum  extent  of  NO  oxidation  is  only  6 percent  and  occurs  at  25  km  when  the  ratio 

“ 3 

[no,L/[oh1  is  reduced  to  10  . Thus  it  appears  that  a very  conservative  upper 

z o o 

limit  on  the  possible  extent  of  NO  depletion  by  the  nitric  acid  reaction  could  be  placed 

X 

at  10  percent. 


To  test  the  consequences  of  the  uncertainties  (about  a factor  of  2)  in  the  rate  constants, 

the  a parameter  was  doubled  (equivalent  to  doubling  k^  or  halving  k^),  producing 

the  dotted  25-km  line  shown  in  Figure  3-9.  The  corresponding  maximum  figure  is  then 

increased  to  11  percent  when  [N02]  /[OH]o  is  10  . In  the  unlikely  event  that  a is 

multiplied  by  4,  the  maximum  goes  up  to  20  percent.  However,  since  only  a fraction  of 

emitted  NO  is  actually  NO  , the  above  conclusion  should  still  be  valid. 

X z 


In  present  engines,  [N O^ J and  [OH]  are  of  comparable  orders  of  magnitude  in  the 
exhaust  jet.  Measures  to  reduce  NO  formation  in  combustion  should  not  affect  the  level 
of  OH,  so  the  condition  [NO(>]q/[OH] q < 10  implies  more  than  a 1,000-fold  reduc- 

tion of  NO  emissions  over  present  engines.  In  that  event  there  is  no  need  to  predict 

X 

the  fraction  of  NO  oxidized  to  nitric  acid  because  there  is  no  NO  problem. 

x x 
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3.4,13  Oxides  of  Sulfur 


At  the  combustor  exit,  both  thermodynamic  and  kinetic  considerations  predict  the  well- 
known  result  that  virtually  all  sulfur  is  present  at  SO  . In  the  nozzle  and  jet , the 
kinetics  of  the  sulfur  oxides  will  be  dominated  by  three  reactions  (Ref.  3-29).  These 
are  the  homogeneous  formation  of  SO 

O 

k 

S02  + O + M -£  SOs  + M (3.27a) 

and  its  reduction  by  H or  O atoms 


S03  + H (or  O)  - S02  + OH  (or  02)  (3.27b) 

This  mechanism  is  analogous  to  the  NO-NO  cycle  in  the  nozzle  and  jet  discussed 

Li 

previously. 

A steady-state  analysis  can  be  applied  as  before,  with  S03  replacing  N02  as  the 
steady-state  species.  The  analysis  gives  the  expression: 

(SOg)  ka  (O)  (M) 

(SOT)  = k [(0)  + (H)] 

ss  a 

15 

Recent  published  values  (Ref.  3-29)  for  k are  in  the  range  k = 1.  5 x 10  to 

1 „16  , 5 , -2  -1.  a a 

10  (cm  -mole  -sec  ) 

For  1^  the  expression  k^  = 6.  5 x 1014  exp  (-10800/RT)  (cm3-mole_1-Sec_1)  has 
been  used.  (Ref.  3-29). 
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,15 


Using  k^  above,  = 4 x 10  , and  nozzle  exit  properties  (T  = 954  K),  the  steady- 

state  relationship  yields 


(SO,)  _4 

= 6 x 10  (nozzle  exit) 


(S02) 


ss 


Therefore,  it  can  be  concluded  that  at  the  nozzle  exit  there  is  only  a negligible 
quantity  of  SC>3 . 

This  is  much  less  than  the  1 to  2 percent  SOQ/SO  ratio  commonly  observed  in  cooler 
combustion  systems  such  as  flue  gases  (Ref.  3-29).  The  explanation  lies  in  the 
10.8-kcal  activation  energy  for  the  SOn  destruction  reaction,  e.g.  , Eq.  (3.27b).  Under 
lower  temperature  conditions,  the  SO^  removal  would  be  greatly  slowed.  For 
example,  if  T = 500°K,  but  the  same  radical  concentrations  are  retained,  then 

(SO  ) 

(soT)  “ °'1  <T  " 500°K> 

^ss 


For  such  conditions  the  assumption  of  steady  state  might  be  invalid. 

It  is  instructive  to  consider  just  the  formation  of  SO^  by  Eq.  (3.27a)  and  neglect  its 
destructions.  If  the  (O)  remains  constant,  this  reaction  has  a half  time 

In  2 

1/2  ka  (O)  (M) 


which  corresponds  to  t = 1 sec.  As  shown  in  Section  3.4.10  the  half  time  for  O atom 

1/ 2 

removal  by  0 + 0+  M-O  + M isa  great  deal  less  - 0(10  ) sec.  This  suggests 

2a  o 

that  O will  not  be  available  in  the  diluted  low -temperature  wake  long  enough  for  a non- 
negligible  fraction  of  SO^  to  be  oxidized  to  SO^. 
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In  the  presence  of  water  vapor,  SO  rapidly  forms  sulfuric  acid  aerosols.  Possible 

O 

sulfur  aerosol  reactions  in  the  wake  will  be  discussed  more  fully  in  Section  4.2,  but  it 
appears  from  the  above  discussion  that  SO  , and  hence  acid  aerosol,  will  be  absent 

O 

from  the  wake. 

3.4.14  Comparison  of  Kinetic  Model  With  Ground  Test  Data 

Test  data  obtained  from  the  wake  of  a J-85  engine  at  sea-level  static  maximum 
A/B  conditions  (Ref.  3-25)  indicate  substantial  CO  oxidation  in  the  exhaust 
jet.  This  has  been  confirmed  by  calculations  with  the  present  kinetic  model  using 
plume  temperature  profiles  given  for  the  J-85  under  similar  conditions.  A comparison 
between  the  test  data  and  computed  CO  concentrations  is  given  in  Figure  3-10.  The 
kinetic  model  calculations  were  carried  out  only  for  a short  distance  beyond  the 
nozzle  exit.  Exact  agreement  should  not  be  expected  since  the  temperature  profile 
used  was  not  confirmed  to  be  compatible  with  the  test.  In  addition,  the  test  data  have 
not  been  corrected  for  dilution  of  the  plume  by  surrounding  air. 

The  rapid  oxidation  of  CO  behind  the  exit  plane  of  a J-85  engine  operating  at  sea  level 
is  quite  a different  result  from  the  slow  change  in  CO  mole  fraction  with  distance 
behind  a GE-4  engine  in  supersonic  flight  in  the  stratosphere.  The  kinetic  model  is 
able  to  account  for  and  explain  the  differences.  In  both  cases  the  model  predicts  almost 
no  conversion  of  NO  to  N09,  however,  as  discussed  previously. 
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Figure  3-10  Comparison  of  Kinetic  Model  With  Ground  Test  Data 


3.4.15  Essential  Exhaust  Reactions 


In  the  light  of  the  complete  analysis,  many  of  the  reactions  in  Table  3-1  can  be 
eliminated  from  consideration.  The  deletions  can  be  divided  into  three  catagories 
First,  some  species  are  present  only  in  negligibly  small  concentrations.  Providing 
this  condition  is  not  the  result  of  a dynamic  balance  which  can  shift  (e.g.  , by  the 
removal  of  a destructive  radical),  these  species  can  be  neglected.  Second,  reactions 
that  are  very  slow  in  both  forward  and  reverse  directions  can  be  eliminated.  Third, 
frequently  one  direction  of  reaction  can  be  neglected  in  comparison  to  the  other. 

In  the  kinetic  calculations  presented  in  this  report  the  reverse  reactions  have  been 
included  automatically,  with  the  reverse  rate  constants  calculated  from  the  forward 
rate  constants  and  thermodynamic  data.  In  the  wake  model  program  under  develop- 
ment, forward  and  reverse  reactions  will  be  treated  as  completely  separate  reactions. 

On  the  basis  of  the  first  two  catagories,  all  reactions  involving  the  species  N and 

N^O  can  be  omitted.  For  all  practical  purposes  N atoms  are  not  present  in  the 

engine  exhaust.  Figure  3-4  and  Table  3-4  show  that  the  N90  is  virtually  frozen 

^ -3 

at  0. 1 ppm  because  of  slow  reactions;  this  is  a tiny  fraction  (10  ) of  the  total  NO  , 

X 

and  hence  N^O  may  be  neglected. 

As  exemplified  by  Table  3-4,  nearly  all  three  body  reactions  fall  into  the  third 
category.  Only  the  recombination  directions  of  these  reactions  need  to  be  considered; 
the  direct  thermal  disassociation  of  all  species  except  HNO^  can  be  neglected. 

A list  of  the  reactions  deemed  essential  for  engine  exhaust  nozzle  and  jet  kinetic 
calculations  is  given  in  Table  3-7.  Forward  and  reverse  reactions  are  tabulated 
separately  and  the  rate  constant  parameters  are  listed  in  both  molar  and  molecular 
units  for  convenience.  The  indicated  references  can  be  found  in  the  reaction  rate 
data  reference  list  (section  3.7)  as  before.  The  stated  uncertainties  in  the  rate 
constants  are  also  included  under  ±log  k column.  No  entry  appears  in  this  column  if 
no  confidence  level  was  given  by  the  evaluator.  In  such  a case,  the  probable  un- 
certainty can  usually  be  taken  to  be  large. 
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ESSENTIAL  EXHAUST  REACTIONS 
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♦References  in  this  Table  appear  as  a special  reference  list  in  Section  3. 
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Table  3-7  may  not  constitute  a completely  sufficient  list  since  unknown  (e.g.  , 
hydrocarbon)  reactions  undoubtedly  are  possible  in  addition  to  the  ones  shown  (e.g.  , 
the  methane  oxidation  sequence).  The  addition  of  other  HC  reactions  would  not  cause 
any  foreseeable  change  of  such  a magnitude  as  to  have  an  environmental  impact.  On 
the  contrary,  the  higher  hydrocarbons  tend  to  have  more  reactive  groups  such  as 
reactive  (a  ) hydrogen,  tending  to  speed  up  termination  reactions.  Thus,  the  previous 
conclusions  attributed  to  rapid  removal  of  inorganic  radicals  by  unburned  hydrocarbons 
are  apt  to  be  even  more  strongly  reinforced  by  inclusion  of  oxidation  mechanisms  for 
the  higher  weight  unburned  hydrocarbon  fractions.  Of  course,  "unknown"  reactions 
need  not  be  limited  to  those  involving  organic  species. 

3.4. 16  Sensitivity  of  Results  to  Input  Data 

The  objective  of  a sensitivity  analysis  is  to  identify  those  factors  to  which  computed 
results  are  most  sensitive.  In  the  present  case,  the  computed  results  will  be  the  final 
wake  concentrations  of  species  with  a potentially  significant  climatic  impact.  Obviously, 
the  stated  objective  cannot  be  achieved  fully  until  the  model  is  developed  and  exercised 
during  the  next  phases  of  the  program.  However,  some  discussion  is  warranted  there, 
because  identification  of  the  essential  features  of  the  model  has  been  dependent  through- 
out on  the  data  that  have  been  consulted. 

The  two  categories  of  results  for  which  sensitivity  analyses  are  needed  are: 

• Dependence  of  model  features  on  uncertainties  in  the  data  underlying  the 
model  formulation 

• Sensitivity  of  results  computed  by  the  model  to  input  data  uncertainties 

In  either  of  the  above  cases,  a systematic  investigation  may  be  divided  into  two 
parts:  (1)  sensitivity  to  uncertainties  in  rate  constants,  and  (2)  sensitivity  to  input 
species  concentrations. 

The  simplified  chemical  model  that  is  currently  being  developed  has  evolved  from 
broader  considerations.  The  most  important  factor  that  can  influence  the  nature  of  the 


3-60 


model  is  the  possible  extent  of  conversion  of  NO  to  HNOQ  . This  is  discussed  in 
section  3.4.12,  where  it  is  seen  that  the  maximum  uncertainties  in  radical  recombin- 
ation rate  constants  are  insufficient  to  allow  for  significant  formation  of  HNO3  in  the 
near  field.  Even  for  a worst -on-worst  situation  corresponding  to  maximum  uncertainty 
in  each  rate  constant,  in  opposite  directions,  and  for  an  unrealistically  cooperative 
assumption  that  all  NC^  is  N02  , the  fraction  converted  to  HNO^  only  begins  to 
approach  significent  values  in  terms  of  climatic  impact  analyses.  If  this  had  not 
been  the  case,  i.e.,  if  HNO^  formation  appeared  to  be  extensive,  then  a much  more 
complicated  and  general  model,  cumbersome  and  expensive  to  run,  might  have  been 
necessary. 

The  model  will  "turn  on"  at  the  nozzle  exit  plane.  However,  the  model  characteristics 
will  depend  on  the  various  simplifying  relationships  discovered  in  the  chemistry  farther 
upstream  — i.e. , within  the  nozzle.  For  purposes  of  illustration,  one  important  factor 
mentioned  previously  — the  NO  which  is  converted  to  NO^  — will  be  discussed  in  terms 
of  its  sensitivity  to  uncertainties  in  a specific  reaction.  This  example  is  chosen  to 
illustrate  the  principles  involved  in  performing  the  sensitivity  analysis. 

Table  3-7  lists  the  set  of  essential  reactions  needed  for  the  kinetics  model  and  shows 
the  rate-constant  uncertainties  given  by  the  evaluators.  Sensitivity  to  rate-constant 
uncertainties  may  be  handled  parametrically  as  illustrated  below  for  one  NO  oxidation 
reaction. 

It  is  known  that,  under  certain  conditions,  the  hydroperoxyl  radical  can  rapidly 
oxidize  nitric  oxide,  as  follows: 


k36 

H02  + NO-~  N02  + OH 


The  value  for  k0„  listed  in  Table  3-7  is  a "latest  upper  limit"  estimate.  The 

^l  2 

published  value  for  k3g  = 1.25  x 10  cm  /mole-sec  had  been  used  previously. 
From  the  computer  output,  the  five-fold  increase  in  k^g  results  in  only  a 0.25 
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increase  in  NOg  concentration.  Thus,  NO^  is  not  found  to  be  particularly  sensitive 
to  kQg  under  the  present  conditions.  As  will  be  seen  below,  this  simple  "substitution 
of  values,  run  the  computer"  approach  is  often  wasteful  and  unnecessary. 

The  sensitivities  of  species  concentrations  to  their  input  values  at  the  throat  can  be 
assessed  readily  with  Figures  3-4  to  3-7.  Where  there  are  reactions  tending  to 
drastically  and  rapidly  alter  initial  concentrations  — illustrated  by  the  strong  non- 
equilibrium production  of  radicals  at  10  cm  downstream  from  the  throat  as  discussed 
in  section  3.4.4  for  example  — the  sensitivity  to  initial  concentration  is  very  small. 
The  rapid  destruction  of  N09  seen  in  Figure  3-4  shows  this  — where  an  artificially 
high  initial  concentration  of  N00  was  used.  Such  cases  can  be  identified  easily  by 
rapid  approach  to  steady-state  levels  that  are  very  different  from  the  starting  con- 
centrations. On  the  other  hand,  species  concentrations  that  do  not  change  appreciably 
from  their  starting  values  may  or  may  not  be  highly  sensitive  to  those  values.  Ob- 
viously, any  inert  species  will  have  concentrations  at  any  point  which  are  linearly 
dependent  on  their  initial  values.  Although  the  NO  profile  in  Figure  3-4  is  flat,  it 
certainly  does  not  follow  that  NO  is  inert.  It  behaves  so  only  because  initial  N02 
concentrations  are  so  small.  If  the  intial  NO^/NO  ratio  were  large,  say  one -half 
or  more,  then  there  would  be  a very  rapid  increase  in  NO  from  its  initial  to  its 
steady-state  value. 

To  aid  in  discussion  of  sensitivity  of  a species  to  local  concentrations  of  a particular 
reactant  of  associated  rate  constants,  or  to  uncertainties  in  such  data,  the  HC>2 
radical  mentioned  above  is  again  a convenient  choice.  It  might  be  supposed,  for 
example,  that  sufficiently  high  levels  of  H00  would  raise  the  NO2  levels.  The 
primary  destruction  reaction  tending  to  keep  N02  small  in  Figure  3-4  is  attack 
of  NC>2  by  H atoms,  which  are  present  in  high  concentration.  However,  in  many 
private  discussions  it  was  frequently  asserted  that  high  concentrations  of  hydrogen 
atoms  implied  high  H09  concentrations  because  of  the  following  rapid  reaction: 

H + 02  + M — H02  + M 


3-62 


The  relationship  of  HO^  to  H can  be  assessed  readily  with  the  help  of  Table  3-4, 
which  shows  that  reactions  21,  22,  and  25  define  the  major  formation  and  destruction 
processes  for  HO?  . In  this  event  the  rate  expression  is 


When  H09  is  in  a steady  state,  however,  the  following  equation  may  be  derived 


Thus,  the  H09  concentration  is  independent  of  H . This  condition  obviates  the 
need  for  a parametric  sensitivity  study  based  on  uncertainties  in  the  concentrations 
of  hydrogen  atoms.  Also,  it  is  clear  that  the  above  equation  defines  the  sensitivity 
of  HC^  formation  to  uncertainties  in  the  rate  data  for  reactions  21,  22,  and  25  . 
Hence,  a parametric  study  is  again  unnecessary. 

Considerations  of  the  above  types  have  been  made  throughout  the  work  to  define  the 
model.  It  is  felt  that,  except  for  "unknown"  reactions,  the  sensitivities  of  the 
chemical  model  features  to  rate  data  have  been  treated  adequately.  As  was  indicated, 
if  OH  recombination  and  N02  oxidation  rates  had  been  vastly  different  than  they 
are,  a multidimensional  model  capable  of  handling,  for  example,  cross  stream  and 
radial  diffusion  reactions  of  OH  and  NO^  might  have  been  needed.  Fortunately,  this 
is  not  the  case,  and  a simple  stream  tube  approach  should  suffice. 

Finally,  sensitivity  of  the  output  results  calculated  by  the  model  to  uncertainties 
in  its  input  data,  such  as  measured  exit -plane  exhaust  concentrations,  can  be 
better  defined  in  subsequent  phases  of  the  program,  after  final  model  programming 
and  executions  for  test  cases  of  interest  are  in  progress. 


d [H02] 


dt  k21  ^ ^°2^  (^22  + k25^  ^ tH°2^ 
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3.5  SIGNIFICANCE  OF  UNBURNED  HYDROCARBONS 


3.5.1  General  and  Summary 

Aircraft  turbojet  engines  without  afterburners  yield  extremely  low  concentrations  of 
unburned  hydrocarbons  (Ref.  3-30).  However,  with  afterburning  "reheat"  the  con- 
centrations of  unburned  hydrocarbons  increase  drastically  (Refs.  3-31  through  3-33). 
With  maximum  reheat,  the  unburned  hydrocarbons  may  equal  carbon  monoxide  con- 
centration levels;  i.  e. , hydrocarbon  levels  of  over  1,  000  ppm  may  be  present  in  the 
exhaust.  Even  at  maximum  reheat,  overall  equivalence  ratios  are  still  considerably 
less  than  1,  and  there  is  ample  oxygen  present  to  burn  these  hydrocarbons  in  the  hot 
exhaust  stream.  In  fact,  their  presence  in  the  first  place  cannot  arise  from  thermo- 
chemical limitations,  but  must  depend  on  reaction  and  mixing  rate  limitations  of 
combustion  in  the  afterburner.  As  seen  previously,  the  concentrations  of  reactive 
species  — i.e. , O,  H,  and  OH  — may  be  high  in  this  exhaust  mixture.  It  appears 
unlikely  that  the  requisite  rate  constants  and  detailed  mechanisms  can  be  assembled 
readily  for  meaningful  chemical-kinetics  numerical  analysis  of  this  complex  system. 
However,  it  is  clear  that  concentrations  of  unburned  hydrocarbons  in  the  final  wake 
will  be  equal  to  or  less  than  values  based  on  their  measured  concentrations  at  the  exit 
plane  after  accounting  for  appropriate  dilution  in  the  wake.  The  question  then  arises 
whether  it  is  necessary  to  specify  these  concentrations  with  any  greater  accuracy  than 
an  appropriate  upper  limit  corresponding  to  freezing  in  the  nozzle  and  subsequent 
dilution.  In  the  following  arguments  it  is  shown  that  the  impact  of  the  added  unburned 
hydrocarbons  on  the  chemical  processes  occurring  in  the  wake,  dispersion,  and  trans- 
port regimes  and,  ultimately,  their  climatic  impact,  is  negligible. 

Four  reaction  categories  are  considered,  which  encompass  the  spectrum  of  thermo- 
chemical hydrocarbon  reactions  in  the  stratospheric  flight  environment.  Results  of 
these  considerations  may  be  summarized  for  each  of  the  four  categories  as  follows: 

• Hydrocarbons  do  not  enter  into  long  chain-length  catalytic  cycles  that  can 
affect  ozone.  The  reactions  degrade  the  hydrocarbon  chains,  and  the 
maximum  total  effect  on  ozone  is  negligible. 
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• Although  hydrocarbon  reactions  that  remove  NO  appear  to  be  possible  in  the 

X 

stratosphere,  for  mechanistic  reasons  and  because  the  total  quantities  of 
added  hydrocarbons  are  dwarfed  by  the  ambient  methane,  the  impact  on  both 
the  local  rate  and  final  total  NO  removal  by  the  hydrocarbons  added  from  the 
exhaust  is  negligible. 

• The  destruction  of  oxygen  atoms  by  organics  in  the  exhaust  is  not  competitive 

with  faster  processes.  Thus,  the  impact  of  unburned  hydrocarbons  on  exhaust 

SO  and  NO  is  small.  If  the  concentrations  of  unsaturated  hydrocarbons  in 
xx 

the  exhaust  are  well  below  the  S09  emission  levels,  the  hydrocarbon  influence 
on  particulate  formation  is  negligible. 

• Reactions  with  organic  radicals  that  can  accelerate  conversion  of  NO  to  NO^ 
in  the  stratosphere  are  of  negligible  importance  for  at  least  four  reasons : 

— Competing  rates  are  faster . 

— Processes  are  noncyclic. 

— Stable  products  are  formed. 

— Hydrocarbons  are  degraded. 

3.5.2  Hydrocarbon/Ozone  Catalytic  Chains 

By  analogy  to  NO  and  HO  catalytic  cycles  that  can  conceptually  destroy  ozone,  one 

X X 

is  tempted  to  investigate  RO  cycles,  i.e., 

X 

RO-  + 03  - R02-  + 02 

ro2;  + o — ro-  + o2 
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\ 


R’  + Og  — *■  RO-  + 09 
RO-  + O — R-  + 02 


The  net  effect  in  either  case  is 


° + °3  - 2°2 

Aside  from  difficulties  in  achieving  the  first  reactions  on  each  of  these  cycles , the 
main  point  to  be  made  is  that  R-  is  not  N or  H ; i.e.  , R is  not  stable  and  cannot 
"go  around  the  cycle"  very  many  times.  First  of  all,  there  is  competition  for  RO  by 
molecular  oxygen,  which  is  in  great  excess  over  or  O : 

RO-  + 02  — aldehyde  + peroxyl 

If  RO  is  CHgO  at  room  temperature* 

— 18  8 

CHgO-  + Og  — - CH90  + HOg  , k = 3 x 10  cm  /molecule  sec 

Although  this  rate  constant  is  as  much  as  six  orders  of  magnitude  lower  than  those 
for  some  of  the  competing  reactions,  the  concentration  of  0()  is  at  least  six  orders 
of  magnitude  higher  than  the  competing  reactant  levels.  For  example,  the  oxygen 
atom  and  hydroxyl  radical  concentrations  are  about  11  orders  of  magnitude  lower 
than  molecular  oxygen  levels  at  20-km  altitude. 


*A11  rate  constants  are  taken  from  the  compilation  in  Ref.  3-34  unless  otherwise 
indicated. 
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The  formaldehyde  formed  by  the  above  reaction  is  relatively  stable  and  is  slowly 
converted  to  carbon  dioxide  by  several  reactions : 

CH  O + OH  — CHO-  + H O , k = 1.6  x l(f 11 

z z 

CH  O + hr  — CHO-  + H 

La 

CHO-  + O — CO  + HO  ■ , k > 2 x l(f15 

La  La 

CO  + OH  — C02  + H 

Secondly,  the  reactions  are  energetic  and  lead  to  decomposition  of  the  radicals. 

For  example,  from  thermochemical  estimates  (Ref.  3-3'5)  one  finds  that  the  second 
reaction  in  the  proposed  catalytic  cycle  is  highly  exothermic: 

RO  ’ + O — *•  RO-  + 09  AH  ~ - 60  kcal/mole 

It  is  unlikely  that  all  of  the  reaction  energy  would  be  partitioned  solely  among  translational 
modes  of  the  products.  It  is  well  known  that  very  little  residual  energy  in  the  alcoxy 
radical  will  result  in  its  decomposition.  Ethoxy  radicals,  for  example,  decompose 
rapidly  with  energy  inputs  as  small  as  10  kcal/mole  (Ref.  3-36) 

ch3ch2o  — ch3-  + ch2o 

The  products  would  react  further,  formaldehyde  going  to  C02  as  previously  discussed, 
and  CH3‘  entering  the  degradative  cycle  thus: 

CH3 • + 02  -*■  CH302‘  , k > 2.5  x 10"15 

AH  = -27.3 
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One  may  write  a large  number  of  mechanisms  for  stepwise  degradation  of  RO.  One 
such  sequence  follows,  where  NO  is  also  present: 

X 


RCH20-  + 02  — RCHO  + HO  ■ 
RCHO  + OH  — RC  = 0 + H90  , 

k ~ 

3 x 10  13 

0 

II 

RC  = 0 + 0£  — RC02  , 

k ~ 

2 x 10"15 

O 

II 

RC02  + no  — rco9  + N02  , 

k ~ 

1 x 10‘12 

RCO-  — * R*  + C02  , 

k ~ 

1 x 1012 

Also,  reaction  of  the  alkoxy  radical  with  nitrogen  dioxide  gives  stable  products. 


3.5.3  Hydrocarbon-Assisted  Nitric  Oxide  Oxidation 
Acceleration  of  the  reaction 


02  + NO  — N02 


by  intermediates  such  as 


°2  + R'  -*  R(Y 

NO  + RO  ■ — N02  + RO* 
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will  affect  ozone  through  NO  photolysis  if  the  net  rate  exceeds  the  rate  of  nitric  oxide 

& 

oxidation  by  ozone,  i.e.  , 


NO  + Os  — N02  + 02 

However,  this  is  not  likely  to  be  the  case  as  indicated  in  the  following  methane-derived 
radical  reactions  (photooxidation  schemes  involving  hydrocarbons  in  NO  oxidation 
have  been  discussed  by  many  authors  - Refs.  3-37  and  3-38): 

CH4  + OH  — CH^  + H20  , 

m + o2  + o — o3  + m, 

CH3  + 02^CH3°2-  • 
ch3o2-  + no  — no2  + CH30-, 

N02  + CHsO-  — CH30N02  , 

NO  + 03  - N02  + 02  , 

The  steady-state  concentrations  of  CH309-  should  be  much  lower  than  the  ambient 
ozone  that  is  about  1 ppm  at  20  km.  Comparison  of  the  first  three  of  these  reactions 
illustrates  this,  because  (O)  ~ (OH)  and  (09)  »(CH  ).  Proper  treatment  would 
of  course  include  all  formation  and  destruction  reactions.  The  third  reaction  is 
necessary  to  convert  the  slowly  formed  CH3  to  CHgO^  ; it  seems  clear  that  CH302' 
will  be  present  only  in  small  concentration.  Comparison  of  the  fourth  and  last  of  the 
foregoing  reactions  shows  that  oxidation  of  NO  by  03  should  be  faster  than  the 
hydrocarbon-assisted  oxidation. 

Attack  of  the  radicals  has  been  discussed  previously,  making  their  regeneration 
and  participation  in  cyclic  peroxyl  oxidation  of  NO  unlikely.  Also,  reaction  products 

can  form  other  stable  products,  as  illustrated  by  the  fourth  and  fifth  of  the  foregoing 
reactions  taken  in  sequence. 


k » 

10  14 

k9.  * 

. Hf15 

(M  taken  as 

k « 

i<f15 

k « 

1(T13 

k » 

10-12 

k « 

io“14 
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Finally,  the  energetic  reactions  lead  to  degradation  of  the  radicals,  removing  them 
from  the  system,  as  discussed  previously. 

3.5.4  Reactions  Leading  to  Stable  Hydrocarbon-NO  r Products 

Reactions  of  NO  that  are  of  interest  are  given  below 
x 


R-  + NO  — RNO  (nitrosoalkane) 

R-  + NO  — - RNO  (nitro alkane) 
z z 

RO-  + NO  — RONO  (alkylni trite) 

RO-  + N02  — R0N02  (alky lnitr ate) 

R-  + N02  — ■ RONO  (alky lnitr ite) 

O O 

n ii 

RC02  + N09  — RC02N09  (peroxyacetyl  nitrate,  PAN) 


(3.28) 

(3.29) 

(3.30) 

(3.31) 

(3.32) 

(3.33) 


Reactions  (3.29)  and  (3.31)  give  stable,  chain  ending  products.  All  the  reactions 
result  in  NO  being  tied  up  with  an  organic  moiety.  Most  of  these  reactions,  excluding 
Eq.  (3.33)  have  been  studied  in  depth  by  Phillips  and  Shaw  (Ref.  3-39).  The  products 
in  Eqs.  (3.28),  (3.30),  and  (3.32)  are  not  particularly  stable  to  ultraviolet  (uv) 
photolysis,  whereas  Eqs.  (3.29)  and  (3.31)  products  are  quite  stable  to  uv.  They  do 

o 

not  absorb  much  energy  above  2,000  A and  can  be  considered  to  be  stable,  long-lived 
reaction  products.  Certainly  the  nitroalkanes  should  be  as  stable  as  PAN.  All  these 
reactions  are  exothermic,  based  on  the  thermochemical  data  given  in  the  treatise  by 
Benson  (Ref.  3-35). 

Oxygen  atom  and  hydroxyl  radicals  can  react  with  the  nitroalkanes  and  the  nitrates 

to  form  nitroalkane  radicals.  These  radical  species  are  relatively  unreactive  and 

should  not  contribute  to  catalytic  chain  cycles  with  ozone.  The  chemistry  of  the  alkyl 

nitrate  radical  (RCHON09)  is  quite  similar  to  that  of  the  nitroalkyl  radical  (RCHNO  ). 

z z 


It  is  conceivable  that  high  concentrations  of  hydrocarbons  in  the  exhaust  relative  to 
ambient  levels  might  lead  to  some  NO  capture.  As  shown  previously,  CH-  might 

X o 

be  present  in  the  exhaust  jet,  leading  to  the  following  reaction: 

CH*  + NO  — CH  NO  AH  = -39.6 

O ^ O 

However,  thermochemical  and  photochemical  stabilities  of  nitrosoalkanes  are  not 
particularly  great.  The  corresponding  methyl  radical  reactions  with  NO  would  be 
limited  by  the  small  N09  concentrations  in  the  near  wake. 

Possible  NO^-removing  reactions  of  unburned  hydrocarbons  in  the  wake  itself  are 
illustrated  below  for  radicals  derived  from  reactions  of  methane  with  ambient  oxidizing 
radicals.  The  rate  constants  are  taken  from  Niki  (Ref.  3-34,),  and  the  thermochemical 
estimates  are  based  on  Benson's  tables  (Ref.  3-35): 

CH^  + N02  - CH3N02  , 

CH^  + N09  — CH30-  + NO  , 

CH30-  + NO  — CHgONO  , 

ch3o-  + no2  — ch3ono9  , 
ch3o2*  + no  — no2  + ch3o  , 

In  the  wake,  relatively  high  concentrations  of  hydrocarbons  might  promote  some  of  the 
foregoing  reactions,  as  illustrated  below 

CH4  + OH  — CH3  + H90  , k * 

CH-  + N02  CH3N02  , k « 


10 


-14 


10 


-12 


-12  3 

k = 2.8  x 10  cm  /molecule  sec 
AH  = -33  kcal/mole 

k = 5.4  x 10"12 
AH  = -17 


k = 6.  7 x 10 
AH  = -41 

-12 

k = 2 x 10 
AH  = -41 


-13 


k = 2 x 10 
AH  = -17 


-13 
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The  second  step  is  fast  but  must  compete  with 


CH3  + o2  - ch3o-  , 

The  latter  reaction  is  very  fast  because  of  the  high  O 2 concentration. 

It  is  important  to  note  the  high  ambient  concentrations  of  methane  - i.e.  , about 
5 x 1012  molecule /cm°  or  1 ppm  at  20  km  (Ref.  3-40).  The  total  methane  contribution 
from  an  SST  fleet  is  estimated  to  add  only  0.5  ppb  to  the  background  (Ref.  3-41). 

Thus,  in  the  long  rim,  the  ambient  atmospheric  reactions  should  dwarf  any  effects  of 
the  hydrocarbons  added  from  the  aircraft  exhaust.  If  ambient  methane  can  be  shown 
to  participate  in  reactions  that  cause  a significant  net  NO^  depletion  the  conclusions 
of  H.  S.  Johnston  and  others  regarding  ozone  destruction  would  have  to  be  modified. 

3.5.5  Comparison  of  Roles  of  NO  and  Hydrocarbons  in  Urban  Atmospheres  With 
Their  Possible  Roles  in  Stratosphere 

It  is  well  known  that  nitrogen  oxides  are  a necessary  ingredient  of  photochemical  smog, 
specifically  ozone  production  in  urban  atmospheres.  The  question  may  arise,  there- 
fore, as  to  why  NO^.  , which  helps  produce  ozone  at  sea  level,  should  help  to  destroy 
it  in  the  stratosphere.  Also,  it  has  been  indicated  that  hydrocarbons  are  not  expected 
to  contribute  significantly  to  stratospheric  smog  even  though  it  is  well  known  that  they 
are  of  key  importance  in  urban  smog  formation.  The  underlying  explanation  stems 
from  the  fact  that  the  intense  uv  radiation  in  the  stratosphere  dissociates  molecular 
oxygen  to  form  O atoms,  whereas  at  sea  level  this  effect  is  of  little  significance. 
However , NO  may  be  photolyzed  in  both  environments  thereby  providing  a potential 
source  of  O atoms  at  sea  level  when  NO^  is  added.  At  sea  level,  the  presence  of 
NO  alone  is  insufficient  for  production  of  ozone  because  photolysis  of  N09  also 
produces  nitric  oxide,  a good  oxygen  atom  scavenger.  Since  most  of  the  NO  is 

X 

added  initially  as  nitric  oxide , it  is  oxidized  by  any  O atoms  produced  in  preference 
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to  oxidation  of  molecular  oxygen  to  form  ozone.  It  is  here  that  the  hydrocarbons  added 
to  urban  atmospheres  play  a major  role  in  formation  of  photochemical  smog.  They 
provide  a path  by  which  molecular  oxygen  can  be  used  to  oxidize  the  nitric  oxide  to 
N09  . Once  the  nitric  oxide  has  been  depleted  sufficiently, the  photolysis  of  NO  and, 
consequently,  production  of  O atoms  leads  to  ozone  formation,  since  there  is  now  a 
greatly  reduced  competition  for  the  O atoms  from  nitric  oxide. 

The  details  of  the  hydrocarbon  role  in  helping  molecular  oxygen  to  convert  NO  to 
NOg  are  complicated.  The  key  species  in  the  oxidation  steps  are  organic  peroxy 
radicals.  Also,  carbon  monoxide  can  play  a role  through  production  of  hydrogen  atoms 
during  its  slow  oxidation.  The  hydrogen  atoms  react  with  molecular  oxygen  and  these 
products,  in  turn,  can  oxidize  carbon  monoxide,  liberating  more  hydrogen  atoms. 

The  key  species  in  this  cycle  is  the  hydroxyl  radical,  which  is  currently  believed  to 
be  the  important  intermediate  in  converting  molecular  oxygen  to  the  organic  peroxy 
radicals. 

Figure  3-11  shows  some  typical  chains  whereby  the  hydrocarbon  species  provide  paths 
for  the  ground  level  oxidation  of  nitric  oxide  and  other  species  initiated  by  hydroxyl 
radicals.  It  is  noteworthy  that  when  the  organic  peroxy  radicals  oxidize  NO  , they 
are  converted  to  alkoxy  radicals.  These  decompose  to  stable  carbonyl  — e.g.  , 
formaldehyde  and  a carbon-centered  radical.  The  reactions  continue  until  degradation 
of  the  original  hydrocarbon  molecule  into  stable  oxidized  products  is  complete. 

Thus,  in  polluted  urban  atmospheres,  there  is  no  ozone  until  the  hydrocarbons  have 
provided  a path  that  converts  nitric  oxide  to  NC^.  It  is  the  photolysis  of  the  NC^ 
and  the  lack  of  competition  from  nitric  oxide,  which  has  been  reduced  to  low  con- 
centrations, that  liberates  and  increases  the  levels  of  oxygen  atoms  and,  consequently, 
the  ozone  level  through  reaction  with  molecular  oxygen.  Eventually,  ozone  begins  to 
react  with  unsaturates.  The  ozonolysis  reactions  lead  to  nitrogen-containing  products. 
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Figure  3-11  Photo-Oxidation  of  trans-but-2-ene  in  the  Presence  of  Nitric  Oxide  and 
Nitrogen  Dioxide:  Skeleton  Reaction  Scheme  Involving  Hydroxyl  Radicals 

, (Taken  From  J.  A.  Kerr  et  al.  , Chemistry  in  Britain.  8,  252  (June  1972) 
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The  situation  in  the  stratosphere  is  quite  different  because  the  initial  ozone  and  oxygen 
levels  are  quite  high  and  the  photodissociation  of  molecular  oxygen  tends  to  keep  them 
high.  Thus,  ample  mechanism  for  oxidizing  nitric  oxide  to  N09  is  already  present 
in  the  stratosphere.  The  oxygen  atom  concentrations  are  sufficiently  high  to  assist  in 
conversions  of  N09  back  to  nitric  oxide  via  the  reaction 

no2  + o — no  + o2 

Thus,  the  production  of  O atoms  from  the  photolysis  of  N09  in  the  stratosphere 
affects  the  destruction  of  ozone  through  the  NO  catalytic  cycle.  The  reasons  for 

X 

this  will  be  seen  more  clearly  in  Section  4. 
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3.6  CHEMICAL  REACTIONS  IN  TURBULENT  MEDIUM 


3.6.1  Introduction 

An  important  question  arising  in  this  and  other  programs  is  relative  to  the  effect  of 
turbulence  upon  the  chemistry.  To  gain  some  insight  into  this  problem  without  claiming 
to  have  solved  it,  we  will  review  basic  work  supported  under  the  Lockheed  Independent 
Research  program  in  connection  with  missile  reentry  applications.  Derivations  and 
results  of  this  work  will  be  presented  elsewhere. 

To  clarify  the  discussion,  a brief  description  is  given  on  the  effect  of  turbulence  in 
a fluid  on  a parcel  of  contaminating  fluid  injected  into  it. 

It  has  long  been  recognized  that  turbulence  causes  concentration  fluctuations.  The 
earliest  related  work  on  this  subject,  which  is  quoted  in  the  recent  literature,  is  over 
20  years  old.  It  is  credited  to  the  Russian  Obukhov  (Ref.  3-42)  and  the  American 
Corrsin  (Ref.  3-43).  Their  work  dealt  directly  with  temperature  fluctuations.  The 
methodology  is  not  changed  significantly  when  one  treats  concentration  fluctuations, 
which  may  be  caused  by  temperature  fluctuations . 

Concentration  fluctuations  are  modeled  in  the  following  manner.  A contaminant  fluid 
is  added  to  a background  fluid  undergoing  turbulent  motion.  Turbulence  is  defined  by 
an  irregular  fluctuating  component  in  the  velocity  field  of  the  background  fluid.  The 
additive  fluid  particles  move  with  the  local  fluid  velocity  because  of  collisions  with 
the  fluid  background.  These  particles  thereby  acquire  an  overall  drifting  component 
of  motion,  due  to  the  mean  flow,  and  an  irregular  motion,  or  jitter,  due  to  the 
fluctuating  velocity  field  of  the  background  fluid.  They  behave  as  a tracer  element  or 
dye  spot  of  the  background  motion.  It  is  assumed  that  the  tracer  does  not  affect  the 
background  motion,  which  is  determined  by  some  external  energy  source.  In  this 
sense  the  tracer  is  called  a dynamically  passive  additive.  Another  source  of  con- 
centration fluctuations  is  due  to  a random  distribution  of  the  initial  concentrations. 

Much  of  the  recent  literature  models  the  turbulent  medium  in  the  limit  of  fast  reaction 
by  ordinary  differential  equations  in  time,  subject  to  these  random  initial  conditions. 

We  also  call  such  models  turbulent,  despite  the  absence  of  the  basic  velocity  con- 
tribution to  real  turbulence. 
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The  tracing  is  not  precise  because  of  diffusion,  which  smears  the  additive  concentration 
and  thereby  dilutes  it.  Since  no  reactions  take  place  in  this  model,  the  total  number  of 
each  smeared  additive  species  is  conserved.  The  model  is  called  the  conserved  passive 
additive  model. 

Experimental  verification  of  this  effect  is  performed  most  easily  when  the  passive 
additive  is  simply  a hot  gas  added  to  the  background  fluid.  The  temperature  fluctuation 
is  then  observable  directly. 

It  is  conceptually  important  to  distinguish  turbulent  diffusion  from  other  basic  processes 
that  appear  to  have  a formal  similarity.  For  example,  there  is  an  analog  between 
turbulent  diffusion  of  a single  fluid  particle  and  Brownian  motion.  When  the  Lagrangian 
viewpoint  is  used  to  follow  the  fluid  particle  in  a homogeneous  turbulent  medium,  a 
long-time  behavior  for  the  mean  distance  traveled  by  the  particle  is  derived  as  pro- 
portional to  the  first  power  of  time  (Ref.  3-44).  This  dependence  is  identical  with  the 
Einstein  result  for  the  mean  displacement  of  a Brownian  particle  (Ref.  3-45).  A Brownian 
particle  such  as  in  an  aerosol  is  larger  in  mass  than  the  molecules  that  surround  it. 

The  molecular  background  bombards  the  Brownian  particle  in  a temporally  random 
manner,  causing  it  to  suffer  small  angle  deflections.  These  deflections  are  thus 
statistical  in  origin  and  are  unimportant  for  the  large  fluid  parcels  considered  in  the 
turbulence  problem.  Thus,  Brownian  motion  per  se  is  not  the  process  underlying 
turbulent  diffusion. 

In  the  work  discussed  here,  the  Eulerian  or  field  viewpoint  is  used  throughout.  The 
mean  displacement  of  a particle  does  not  appear  in  this  formalism.  This  displacement 
is  not  to  be  confused  with  a characteristic  "local”  length  scale  defined  in  terms  of  one 
of  the  partial  differential  equations  governing  the  dynamics  of  the  turbulence . One 
such,  for  example,  is  the  concentration  microscale  (Ref.  3-46)  which  is  defined  in  a 
manner  analogous  to  the  Taylor  microscale,  or  dissipation  length,  for  velocity 
fluctuation  (Refs.  3-47  and  3-48).  Finally,  a calculation  describing  the  evolution  of 
a contaminant  spot  size  with  time,  based  on  the  Eulerian  viewpoint,  shows  that  when 
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strain  is  present  the  spot  area  grows  much  faster  than  the  first  power  of  time 

(Ref.  3-46).  It  is  this  time  dependence  that  characterized  the  homogeneous  turbulence 

result  as  well  as  Einstein's  Brownian  motion  result. 

By  1958  Corrsin  (Ref.  3-49)  had  included  the  case  of  a single  chemical  reactant  in 
the  preceding  model  of  the  dynamically  passive  additive.  Formally,  this  change 
requires  only  a simple  modification  of  the  equations  governing  the  passive  additive 
model. 

More  recently,  O'Brien,  1971  (Ref.  3-50)  considered  the  problem  of  two  reactants 
undergoing  a second- order  reaction.  A significant  result  of  his  work  (Refs.  3-50  and 
3-51)  is  that  the  concentration  fluctuations  can  be  large  compared  with  mean 
concentration  values . 

In  a quiescent  gas  it  is  the  mean  concentrations  that  are  measured  directly.  The 
equations  governing  the  time  and  space  behavior  of  the  concentration  (often  assumed 
to  be  spatially  uniform)  are  in  this  case  governed  by  the  mean  quantities  alone.  When 
turbulence  is  present,  the  time  and  space  behavior  of  the  mean  quantities  depends  on  the 
correlations  between  the  fluctuating  components  of  the  concentration. 

Later  work  by  du  P.  Donaldson  and  Hilst  (Refs.  3-52  and  3-53)  includes  a discussion 
of  the  possibility  of  fluctuations  in  the  reaction  rate  constant  affecting  evolution  of  the 
mean  concentrations.  Although  we  do  not  agree  with  the  equations  of  Ref.  3-52,  the 
concept  must  be  examined  because  of  the  exponential  dependence  of  the  backward 
reaction  upon  temperature. 

The  idea  is  as  follows.  The  reaction  rate  constant  is  a function  of  temperature  — to  a 
significant  degree  in  the  bond-breaking  direction.  The  temperature  of  a gas  is  a 
macroscopic  quantity  defined  by  an  average  over  a one-particle  distribution  for  the 
molecular  velocities.  (It  is  presumed  that  a large  number  of  molecules  are  involved 
in  this  averaging  process. ) Then  the  reaction  rate  constant  is  itself  a macroscopically 
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defined  quantity.  Because  of  the  energy  put  into  the  fluid  system  (e.g. , in  the  form  of 
a shearing  flow),  large  local  temperature  fluctuations  can  exist  in  the  fluid.  These  can, 
in  turn,  cause  the  reaction  rate  constant  k to  vary  from  point  to  point.  It  may  be 
convenient  to  express  this  variation  in  k(T)  as  a part  due  to  the  mean  temperature 
T and  a part  due  to  the  fluctuations  in  temperature  AT  . 


Then 


k(T)  = k (T  + AT) 


k (T) 


+ ^ 
3T 


_ AT 
T = T 


= k + k' 


where 


k = k (T)  , 


k’ 


3k| 

aTl  T=f 


AT 


At  this  point,  we  rename  k(T)  as  k^^  to  emphasize  that,  although  it  contains  the 

ordinary  reaction  rate  constant  k(T)  , which  appears  in  the  chemical  kinetics  literature, 

k,  , is  not  this  familiar  rate  constant,  and  should  not  be  confused  with  it.  Care 
turb  ’ 

should  be  exercised  in  the  literature  to  avoid  confusion. 


3.6.2  Discussion  of  Previous  Work 


The  problem  of  chemical  reaction  of  two  species  in  a turbulent  flow  is  posed  as  a set 
of  partial  differential  equations  subject  to  appropriate  boundary  conditions.  In 
turbulent  flow  one  measures  the  average  values  of  quantities  that  have  a fluctuating 
component.  In  general,  the  probability  distribution  function  for  the  fluctuating  field 
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of  velocity,  pressure,  and  concentration  is  unknown.  Thus,  the  averaged  quantities 
replace  the  precise  dependent  variables  that  are  solutions  of  the  set  of  partial  dif- 
ferential equations.  Unfortunately,  because  of  the  nonlinearity  of  the  convective  time 
derivative  and,  with  second-order  chemical  reactions,  the  nonlinear  coupling  of  the 
concentration  variables,  the  number  of  averaged  equations  is  less  than  the  number  of 
new  unknown  dependent  variables  that  are  the  averaged  quantities . The  problem  of 
closing  the  set  of  equations  (known  as  the  closure  problem)  has  beset  the  theory  of 
the  turbulence  since  its  beginning. 

The  traditional  method  of  resolving  the  closure  problem  is  to  relate  averages  of 
products  of  n-fluctuating  quantities  to  averages  of  products  of  (n-l)-fluctuating  quantities 
by  means  of  "mixing  length"  parameters  that  are  determined  experimentally.  Donaldson 
and  Hilst  (Ref.  3-53)  use  this  approach.  O'Brien  (Ref.  3-50),  on  the  other  hand, 
assumes  a closure  hypothesis  for  the  second-order  reaction  term  that  contains  no 
arbitrary  parameters  with  which  to  fit  the  data.  It  appears,  however,  that  in  treating 
the  diffusion  terms  for  the  full  diffusion-reaction  problem,  O'Brien  must  introduce 
the  traditional  mixing  length  approximation;  however,  he  does  not  discuss  this  problem. 
For  his  closure  scheme,  he  claims  the  advantage  that  it  guarantees  the  required 
positive  definiteness  of  the  mean-square  concentration  fluctuations. 

O'Brien  has  achieved  some  numerical  solutions  for  the  mean  concentration  and  the 
root-mean-square  (rms)  fluctuations  in  concentration  for  the  case  in  which  diffusion 
effects  are  negligible.  He  finds  the  rms  fluctuations  can  be  large.  In  his  work  no 
consideration  is  given  to  the  possibility  of  thermal  fluctuation  contributions  to  the 
reaction  rate  constant.  Donaldson  and  Hilst  discuss  briefly  the  problem  of  fluctuation 
contributions  to  the  reaction  rate  constant,  but  their  formulation  does  not  take  into 
account  the  fact  that  the  backward  and  forward  reaction  constants  have  significantly 
different  temperature  dependencies.  They  have  not  yet  obtained  solutions  for  the 
concentration  variable. 


3-80 


3.6.3  Characteristic  Times 


It  is  possible  to  estimate  the  importance  of  diffusion  relative  to  chemical  reaction  in 
determining  the  behavior  with  time  for  the  concentration  of  one  of  two  species  reacting 
in  the  presence  of  turbulence.  This  crude  estimate  is  based  on  an  equation  describing 
the  evolution  in  time  of  the  single  space  point  correlation  for  the  fluctuating  compo- 
nents of  the  two  species  . Here  the  subscripts  denote  the  two  reactants; 

the  bar  denotes  an  averaging  process,  and  the  primed  quantities  are  the  fluctuating 
parts  of  the  respective  concentration  variables.  This  equation  can  be  written 
symbolically  in  the  form 
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(3.34) 


where  the  first  term  on  the  right-hand  side  (rhs)  of  the  equation  denotes  the  diffusion 
terms,  and  the  second  term  denotes  the  chemistry  terms.  The  diffusion  terms  are 
proportional  to  the  diffusion  coefficient  D , while  the  chemical  reaction  terms  are 
proportional  to  the  reaction  coefficients  k.  (i  = 1,2). 


The  size  of  the  diffusion  term  in  this  equation  is  next  compared  with  that  of  the 
chemistry  term.  For  this  purpose  it  is  convenient  to  use  Eqs.  (15)  and  (20)  of 
Ref.  3-53  for  the  rhs  of  Eq.  (3.34).  An  approximation  is  sought  for  the  rhs  to  reduce 
Eq.  (3.34)  to  the  form 
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CiC2 


+ T 


> Ci  °2 


(3.35a) 


Integrating,  one  obtains 

cjc£(t)  = CVC’T~=T)  e t/T  (3.35b) 

-1  -1  -1  -1 

where  T = (r , + t ) , r is  the  rate  of  decay  of  the  correlation  due  to  diffusion, 

d c d 

and  t is  the  rate  of  decay  due  to  the  chemistry, 
c 
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The  relative  importance  of  the  two  terms  in  Eq.  (3.35a)  is  evidently  determined  by  the 
relative  magnitude  of  the  two  rates  or,  equivalently,  by  the  ratio  N (Ref.  3-54) 


(3.36) 


N » 1 implies  a diffusion  limited  reaction,  and  N « 1 corresponds  to  a reaction 
controlled  by  the  chemistry,  because  the  slowest  process  limits  the  rate  at  which  a 
reaction  can  go. 


The  validity  of  approximating  Eq.  (3.34)  by  Eq.  (3.35a)  depends  on  the  justification  for 
discarding  certain  terms  in  the  exact  equation.  In  the  case  of  the  chemistry  terms, 
this  approximation  amounts  to  a prescription  for  discarding  all  the  self-correlative 
terms  for  the  concentration  fluctuations.  This  approximation  scheme  is  not  founded 
on  physical  grounds.  It  is  to  be  justified  a posteriori  by  comparison  with  the  results 
of  a more  precise  calculation.  Its  usefulness  resides  in  the  fact  that  it  leads  to  a 
tractable  way  of  obtaining  some  needed  numbers. 


Thus,  one  obtains  the  approximate  equality 


*2°1  +klC2>°iC2  " Tclci°2 


(3.37) 


where  r = (k  C + k C ) 1,  k = k (M  /M  ) and  M.  is  the  molecular  weight  of 
C 2 1 1 2 2 12  1 1 

species  i . 


The  expression  above  for  t is  identical  with  the  r obtained  in  Ref.  3-54.  In  the 

c c 

present  work,  t was  obtained  by  approximating  the  chemistry  term  of  Eq.  (3.34) 

6 C’  C’  /<5t  , in  the  time  rate  of  change  of  C'  C'  . The  authors  of  Ref.  3-54  obtained 

r as  an  approximation  to  the  time  rate  of  change  of  the  quantity  C C . Thus, 

C 12 
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the  Tc  as  derived  by  these  authors  is  not  appropriate  for  comparison  of  the  chemistry 
term  with  the  diffusion  term  of  Eq.  (3.34). 


imation  to  Eq.  (3.  35b).  This  approximation  agrees  with  Eq.  (3.35b)  for  M1  ^ M9 

and  for  mean  concentrations  such  that  (a)  C9  or  (b)  where  one  species  is  dominant. 

say,  C « C . For  either  case  (a)  or  (b)  the  ratio  R of  Ref.  3-53  (R  = t'  /r  ) is 
2 1 c d 

closely  related  to  N . For  M = Mn,  k = k . Then  r'  ^ r and  N as  R-l 

12  12  c c 

The  characteristic  time  for  mixing  due  to  diffusion  is  obtained  from  the  same 
equation  as  that  appearing  in  Refs . 3-53  and  3-54 . Discarding  the  term  with  a 
positive  coefficient,  one  obtains 


where  the  microscale  A is  called  the  dissipation  length.  It  is  chosen  identical  with 
the  microscale  (Refs.  3-46  and  3-47)  one  would  obtain  for  a passive  additive  for  the 
case  where  only  a single  species  is  present. 

One  finds,  therefore,  that  (Ref.  3-45) 


In  another  article  (Ref.  3-53 ) these  authors  use  a different  estimate  for  r without 

— -1  c 

attempting  to  justify  it  [ = (k^C^)  ] . This  expression  is  a nonsymmetrical  approx- 


t  . = A2/2D 
d 


(3.38) 


N = A 


(3.39) 


2D 


To  evaluate  this  expression  one  must  reduce  A to  directly  measurable  quantities. 
For  pure  mixing  (i.e. , without  reaction)  Corrsin  (Ref.  3-49)  has  shown 


where  is  the  Taylor  microscale  for  turbulence,  a is  the  Schmidt  number 
{<j  = v / D),  v is  the  kinetic  viscosity,  and  D is  the  diffusion  coefficient. 


For  molecules  of  equal  mass  (M^  = M^) 
kinetic  theory  shows  D = v (Ref.  3-53). 

A « 


and  approximately  equal  size,  classical 
Thus 


Now  At  is  related  to  the  integral  or  outer  scale  A of  turbulence  by  the  relation 
(Ref.  3-46) 


At/A  =*  (15)1/2 


Re 


-1/2 

A 


where  the  turbulence  Reynolds  number  for  the  outer  scale  is  given  by 


Re 


A 


_ ju'2/  A 


and 


p)1/2  i 


is  the  observed  rms  fluctuation  velocity.  Thus 


_ 1/2 

x’2  = m A • A~2 


(3.40) 


Equation  (3.40)  is,  essentially  a result  quoted  in  Ref.  3-53.  The  expression  following 
Eq.  (21)  of  this  reference  differs  from  Eq.  (3.40)  only  by  the  numerical  factor 
which  they  taken  as  0.05.  Since  these  results  are  based  on  order  of  magnitude  approxima- 
tions, the  present  result  can  be  said  to  agree  with  their  estimate  for  A . The  numerical 
value  chosen  in  Ref.  3-53  for  quantities  on  the  rhs  of  Eq.  (3.40)  are  in  reasonable 
agreement  with  other  estimates  of  atmospheric  numbers  (Ref.  3-46). 
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In  Ref.  3-53  molecules  of  approximately  equal  size  and  mass  ( v = D)  are  assumed. 
Their  choice  for  the  outer  scale  is  A = 10  m , whereas  Ref.  3-46  uses  A = 100  m 
for  essentially  the  same  numerical  values  for  the  other  parameters.  Examination 
of  the  U.S.  Standard  Atmosphere,  1962,  shows  that  the  Ref.  3-53  choice  for  kinematic 
viscosity  and  atmospheric  density  corresponds  to  an  altitude  of  3 km. 

Thus.  Eq.  (3.39)  reduces  to 


where  k^  - k^  - k . This  equation  is  useful  in  making  estimates  as  to  whether  or  not 
a reaction  is  diffusion  limited. 

3.6.4  Application  to  Jets  and  Wakes 

It  is  appropriate  to  apply  Eq.  (3.41)  to  chemical  reactions  occurring  in  a high- 
temperature  supersonic  exhaust  jet.  Qualitatively,  the  following  remarks  can  be 
made:  the  reaction  rate  constant  will  tend  to  increase  with  temperature  above 
atmospheric  values,  the  outer  scale,  corresponding  to  energy  containing  eddies,  will 
be  of  the  order  of  the  diameter  of  the  exhaust  jet.  At  typical  stratospheric  flight 
altitudes  this  length  is  much  less  than  the  atmospheric  outer  scale.  Equation  (3.41) 
shows  that  these  two  effects  tend  to  offset  each  other.  Thus,  a reaction  that  is  diffusion 
limited  in  the  atmosphere  may  tend  to  remain  so  in  an  engine  exhaust.  Of  course,  a 
quantitative  statement  depends  on  the  specific  numbers. 

A turbulent  jet  is  considered  with  mean  exhaust  velocity  of  1,  600  m/sec  and  jet 
diameter  of  1 m,  which  is  set  equal  to  A . A nominal  value  of  turbulence  intensity 
for  the  rms  velocity  (u'^j^^/(U)  = 0.01  is  assumed  near  the  exit  plane  of  the  nozzle, 
within  the  cone  of  constant  properties. 


(3.41) 
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This  number  is  not  inconsistent  with  subsonic  data  (Refs.  3-56  and  3-57).  Equation 
(3.41)  is  applied  to  the  supersonic  case,  recognizing  that  the  theory  of  this  equation 
was  based  on  an  incompressible  fluid  model.  For  this  case,  N reduces  to 
N-  k (Cx  +C2). 

The  following  thermal  reactions  will  be  considered,  using  the  reaction  rate  constants 
of  Section  3 and  exit  plane  concentrations  appropriate  to  the  flow.  N can  then  be 
calculated  readily: 


Reactions 

k (ppm-sec)  1 

C1  (ppm) 

C 2 (ppm) 

N 

1 

2 

CO  + 

OH  — C02 

+ 

H 

0.23 

820 

197 

2.3  x 103 

H + 

X 

O 

t 

C^I 

O 

+ 

O 

0.022 

380 

5.5  x 104 

1.2  xlO4 

H + 

N02  — NO 

+ 

OH 

19.0 

380 

3.8  x 10“3 

7.2  x 104 

It  can  be  seen  that  all  three  reactions  are  diffusion  limited.  The  first  reaction  at 
these  near -jet  conditions  has  about  ten  times  the  value  of  N as  does  this  reaction 
under  the  atmospheric  conditions  of  Refs.  3-53  and  3-54.  Thus,  it  remains  diffusion 
limited  — in  agreement  with  the  preceding  qualitative  remarks.  Note  that  in  the  far 
field  behind  the  jet,  aircraft- induced  fluctuations  have  diminished  to  atmospheric 
levels.  (Here  results  like  those  of  Donaldson  and  Hilst,  Ref.  3-53  , Table  1, 
appropriately  modified  to  stratospheric  altitudes,  apply.) 

A criticism  of  Eq.  (3.41)  and  of  the  results  of  Refs.  3-53  and  3-54  arises  from  use 
of  a diffusion  coefficient  corresponding  to  reactants  of  equal  mass.  Because  of  its 
small  mass,  hydrogen  for  example  can  be  expected  to  diffuse  more  rapidly  than  a 
heavier  molecule. 

As  a possible  candidate  for  a chemically  limited  process,  a reaction  is  considered  that 
occurs  in  the  wake  of  a Mach  2 aircraft.  Consider  the  reaction 

NO  + Os  - N02  + 02 
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where  k = 5 x 10  (ppm-sec)  , Cq3  = 1 ppm,  and  CNQ  = 5 ppm.  The  fluctuation 
intensity  is  estimated  as  follows.  Measurements  on  a two-dimensional  Mach  3 wake 
are  performed  in  Ref.  3-49 . Longitudinal  velocity  fluctuations  relative  to  the  mean 
velocity  values  are  about  7 percent  in  the  turbulent  wake.  It  is  assumed  that  the 
mean  flow  velocity  is  given  by  U(x)  ^ U^q(j  /10  =*  60  m/sec.  Thus,  the  rms  velocity 
fluctuation  is  approximately  40  m/sec.  An  outer  scale  of  about  0.3  m arising  from 
the  aircraft  turbulent  boundary  layer  is  appropriate  to  this  problem.  Then  N =* 

_3 

4 x 10  . This  regime  is  clearly  chemistry  limited. 

3.6.5  Conclusions 

Turbulence  produces  measurable  effects  on  a passive  additive  such  as  the  concentration 
of  a chemical  injected  into  a turbulent  medium.  Such  effects  have  been  observed  in 
measurements  made  on  temperature  fluctuations  of  a hot  gas.  The  importance  of  such 
turbulence  effects  to  the  time  and  space  behavior  of  mean  concentration  of  the  chemical 
components  is  uncertain.  Two  kinds  of  effects  are  envisioned.  O'Brien's  work 
(Ref.  3-51)  indicates  that  the  concentration  fluctuations  themselves  are  large  in  the 
presence  of  chemical  reactions.  They  may,  therefore,  have  a large  influence  on  the 
mean  concentrations  one  wishes  to  predict.  Further,  effective  turbulent  reaction 
rates  may  differ  significantly  from  quiescent  laboratory  values.  Even  in  a well-mixed 
gas,  turbulence  may  speed  up  or  slow  down  reactions  relative  to  the  usual  rates. 

It  has  been  shown  that  diffusion-limited  reactions  may  exist  in  a near -jet  field,  and 
chemically  limited  reactions  may  occur  in  a wake.  One  cannot  yet  say  with  confidence 
how  these  results  apply  to  the  supersonic  jet  aircraft.  This  statement  may  appear 
surprising  in  view  of  the  following  result.  In  a one-component  irreversible  second- 
order  reaction  modeled  in  the  limit  of  fast  reaction  by  a purely  time  dependent  equation 
for  the  concentration,  the  reaction  is  slowed  by  the  mean  square  random  fluctuations 
of  the  concentration.  For  a two-component  system,  however,  the  correlation  function 
for  the  two  components  can  take  on  negative  as  well  as  positive  values  as  it  evolves  in 
time.  Thus,  without  numerical  analysis  of  specific  conditions  it  is  impossible  to 
predict  whether  or  not  turbulence  speeds  up  or  slows  down  the  reactions. 
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The  preceding  review  was  based  on  and  supported  by  ongoing  Independent  Research 
at  LMSC.  Results  of  this  work  will  be  presented  in  a forthcoming  report.  Although 
there  has  been  progress,  it  is  felt  that  work  in  this  field  has  not  yet  reached  the  stage 
where  attempts  to  obtain  solutions  for  practical  jet  engine  exhaust  kinetics  problems 
are  warranted. 
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Section  4 

PHOTOCHEMISTRY  IN  THE  WAKE  REGIME 


4. 1 TIME-DEPENDENT  CATALYTIC  REDUCTION  OF  STRATOSPHERIC  OZONE  IN  A 
THIN  STRIP 


Following  injection  of  NO  contained  in  an  aircraft  wake  into  the  stratosphere,  several 
reactions  begin  to  take  place.  The  two  which  most  affect  ambient  species  concentrations 
are  the  oxidation  of  nitric  oxide  by  ozone,  which  depletes  ambient  ozone,  and  the  photol- 
ysis of  nitrogen  dioxide,  which  increases  oxygen  atom  concentrations.  Near  the  ex- 
haust jet  region  where  the  local  NO  concentration  is  higher  than  that  of  ozone,  the 
reaction  is  limited  by  ozone,  and  there  is  little  nitrogen  dioxide  formed  relative  to 
the  nitric  oxide  level.  Thus,  there  is  an  initial  period  of  adjustment,  during  which 
the  wake  is  diluted,  where  two  things  occur:  (1)  the  NOx  concentrations  decrease 
to  below  the  ambient  Og  levels,  and  (2)  the  NO/NO2  ratio  increases  toward  its 
proper  "photochemical  steady  state"  value.  As  has  been  shown  in  Section  3,  this 
period  of  adjustment  is  slow  with  regard  to  the  wake  time  frame  of  a few  minutes. 
However,  it  is  relatively  fast  in  comparison  with  elapsed  times  of  several  hours  and 
beyond,  so  that  a convenient  starting  point  for  analysis  of  the  progress  of  photochemical 
catalytic  effects  on  ambient  ozone  is  suggested,  namely,  to  start  with  a uniformly  mixed 
thin  layer  or  ribbon  of  NO  , irradiated  from  above,  and  at  concentrations  below  that 

X 

of  the  ambient  ozone  in  the  strip.  This  is  the  "ozone  in  a box"  problem,  or,  rather, 
"ozone  in  an  irradiated  ribbon.  " 

Interdiffusion  with  the  surrounding  environment  does  proceed  while  the  photochemical 
reactions  are  taking  place.  However,  it  is  instructive  to  examine  the  maximum 
effects  of  the  catalytic  cycle  for  the  case  where  diffusion  times  are  long  compared 
with  reaction  times.  The  analysis  can  be  performed  in  a relatively  simple  manner 
for  this  artificial  "closed  system. " The  results  will  be  useful  to  show  trends  and 
relationships  among  the  variables,  and  will  approach  quantitative  accuracy  for  those 
limited  conditions  under  which  the  closed  system  assumption  becomes  more  nearly 
satisfied. 
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An  analysis  has  been  carried  out  to  determine  the  ozone  concentration  as  a function 
of  time  for  injection  of  a given  concentration  of  NO  in  a thin  strip  at  various  altitudes 
in  the  stratosphere.  The  analysis  makes  use  of  several  justifiable  simplifications  in 
order  to  avoid  the  complications  of  integrating  the  full  set  of  differential  equations 
governing  stratospheric  photochemistry  in  the  presence  of  oxides  of  nitrogen. 

4.1.1  Reactions 

The  important  chemical  reactions  governing  catalytic  reduction  of  ozone  by  oxides  of 
nitrogen  have  been  given  by  Johnston  (Ref.  4-1)  and  are  as  follows: 


(a) 

°2  + hv 

2 0 

(b) 

O + O2  + M 

°3  + 

M 

(c) 

O3  + h^  — 

O2  + 

0 

(d) 
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NO  + 
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NO2  + hy  — 
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NO  + NO  + 02  — 
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+ no2 

(j) 

O + NO  + M — 

no2 

+ M 

(k) 

no2  + o3  - 

NO3 

+ O2 

(1) 

NO3  + hy 

NO  + 

o2 

The  rate  expressions  and  rate  constants  for  these  reactions  have  also  been  given 
by  Johnston.  Reaction  (1)  is  much  faster  than  (k)  so  that  in  daylight  N03  is  in 
steady  state  and  reaction  (1)  can  be  eliminated. 

Also,  at  20  km  and  above  reactions  (i),  (j),  and  (k)  a re  relatively  slow  compared  to 
(f),  (g),  and  (h)  and  do  not  greatly  effect  either  the  rate  of  ozone  destruction  or  the 
steady-state  ozone  levels  (Refs.  4-1,  4-2).  With  these  simplifications,  the  above  set 
is  reduced  to  the  eight  reactions  (a)  through  (h). 
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For  the  low  altitude  (20  to  25  km)  regions  of  the  stratosphere  of  interest  here,  the 
ozone  destroying  reactions  (d)  and  (e)  are  very  much  slower  than  the  catalytic  de- 
struction reactions  (f)  and  (g)  and  are  therefore  neglected  in  the  present  analysis. 
Thus,  the  set  of  governing  reactions  is  reduced  to  the  following  six: 


(a) 

02  + hr 
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2 0 

(b) 

O + 02  + M 
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°3  +M 
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O3  + hy 
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°2+  O 

(f) 

NO  + 03 
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N°2  + 02 

(g) 

0 + no2 
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NO  + 02 

(h) 

N02  + hy 

— 

NO  + O 

4.1.2  Steady-State  Simplifications 

In  complex  chemical  kinetic  analyses,  the  steady-state  assumption  is  often  applied 
to  minor  species  with  fast  reaction  rates.  Basically,  what  is  assumed  is  that  such 
a species  can,  by  virtue  of  its  fast  production  and  destruction  rates,  instantly  adjust 
its  concentration  level  to  any  changes  of  major  species.  The  explicit  dependence 
of  concentration  upon  time  for  a steady-state  species  is  thus  eliminated  in  favor  of 
an  implicit  dependence  through  the  major  species. 

In-the  above  reaction  scheme,  it  is  found  that  the  O atoms  can  rapidly  adjust  to 
changes  in  the  ozone  or  NOx  levels  and  thus  qualify  as  a steady-state  species.  Also, 
the  NO  and  NO2  may  be  taken  to  be  in  steady  state  provided  that  [NO]  « [O^]  . 
Thus,  the  further  simplification  that  both  O and  NO/NO2  are  in  steady  state  is 
made  here. 

The  steady-state  assumption  has  been  confirmed  by  numerical  integration  of  the 

reaction  set  (a)  through  (h)  using  the  rates  given  by  Johnston.  It  has  been  found 

that  [O]  and  [N0]/[N02]  reach  their  steady-state  values  before  significant  O3 

depletion  occurs,  provided  that  [NO  ] « [ Oq]  • 

x <J 
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4. 1.  3 Governing  Equations 


The  set  of  reactions  (a)  through  (h)  leads  to  the  following  rate  expressions 

d[  03  ] 

— J = k,,  1 02  ] [ O 1 l M ] - kf  [ NO  ] [ 03  ] - jc  [ 03  J 

= 2 ia[  °2  I - kb  1 °2  ) 1 ° H M J + Jc[  °3  I 
kg  [ O J [ N02  J + jh  [ N02  J 

d [ N02  ] 

^ = kj [ NO  J [ 03  ] - kg[  O ] [ N02  ] - jh[  N02  ] 

The  steady  state  condition  on  O and  N09  gives 

d[  O j = d [ NO 2 ] = Q 

dt  dt 

Setting  these  two  rates  equal  to  zero  and  adding 

d[  O]  d[  NQ2  ] 
dt  dt 

gives 

kb  l °2  ] [ O ] [ M j - jc  [ 03  ] - kf  [ NO  ] [ 03  ] = 2 ja  [ 02  ] - 

2 kg  [ O ] [ N02 

Substituting  this  result  into  the  ozone  rate  gives: 

d[  03  J 

ST-  = 2 Ja  1 °2  1 - 2kgfO]  [N02] 


(4.1) 
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With  the  exception  of  the  small  contribution  from  reaction  (e)  which  has  been  neglected 
here,  this  expression  is  the  same  as  given  by  Johnston.  (See  Ref.  4-1. ) 

It  now  remains  to  write  the  above  rate  as  a function  of  ozone  only  so  that  a simple 
numerical  integration  will  give  ozone  as  a function  of  time. 

From 


d[  Q] 

dt 


0 


(4.2) 


[Oj  = 


jc  1 °3  1 + Jh  t N02  ] 

kb  [ M ] [ °2  ] 


and  from 


Here 


V°i 


d[  NQ2  ] 
dt 


[N02]  = 


« (kf[03]  + jh) 


0 using  [ NO  ] = [ NOx  ] - [ N02  ] 

kf  [ °3  ] + kg  [ 0 ] + jh 

so  that 


[ no2  ] 


kf[Os]  [NOx] 
kf  f °3  1 + ^h 


(4.3) 


Substitution  of  Eq.  (4.3)  into  (4.2),  and  Eqs.  (4.2)  and  (4.3)  into  (4.  1)  gives  the  final 
expression  for  the  above  rate 


dt°3] 

dt 


F(L  > j, 


[ 02  ] , [ M ] , [ Os  ] [ NOx  ]) 
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For  a given  altitude,  all  parameters  are  constant  except  [Og]  , and  the  above 
expression  can  be  numerically  integrated  to  give  [O3]  as  a function  of  time 


[ Or 


d[  Os  ] 


dt 


Oc 


4.  1.  4 Results 


The  curves  of  [0Q]  as  a function  of  time  for  various  [NO  ] levels  and  altitudes  are 

presented  in  Figures  4-1  and  4-2.  The  rate  constants  and  atmospheric  properties 

for  these  calculations  have  been  taken  from  Johnston  and  are  given  in  Table  4-1. 

The  photolysis  rate  constants  are  24-hour  averages  as  computed  by  Johnston  for 

a = [NO  ] /[M]  = 0.  The  values  for  a = 0 have  been  used  since,  for  a thin  strip 
x 

of  NOx  , the  incident  radiation  is  that  of  an  unperturbed  stratosphere.  Some  results 
for  daytime  12-hr  average  radiation  are  also  shown.  The  effect  of  the  addition  of  NO 
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Table  4-1 

PARAMETERS  FOR  PHOTOCHEMICAL  CALCULATIONS 


Parameters 

20  km  Altitude 

25  km  Altitude 

ja 

Hf14-49 

i<r12-98 

I—1 

0 

1 

CO 

to 

cn 

00 

lo'32-63 

jc 

kf 

10-5.73 

i<f3'72 

lo'15-10 

i<f15-01 

k 

lo"14-40 

10-14.34 

g 

Jh 

Hf2-35 

Hf2-35 

log10  IMj 

18.27 

17.93 
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Figure  4-2  Ozone  Reduction  at  20  km 


is  to  reduce  the  concentration  in  the  narrow  strip  to  a new  steady-state  value. 
Higher  concentrations  of  NO  reduce  the  ozone  faster  and  to  lower  steady-state  con- 
centrations.  The  dashed  curves  of  "validity"  simply  marks  the  concentration  regions 
where  ozone  falls  below  the  NO  levels,  and  the  analysis  becomes  invalid. 

The  main  conclusion  to  be  drawn  from  Figures  4-1  and  4-2  is  that  the  ozone  level  is 
not  reduced  in  the  wake  time  frame  O (10  sec).  The  characteristic  times  for  the 
catalytic  reduction  cycle  to  have  a significant  effect  are  much  longer.  It  would  take 
very  large  changes  in  the  rate  constants  to  invalidate  this  conclusion. 

The  ozone  half-times  (time  for  [O^  ] to  be  reduced  to  one-half  its  initial  value)  are 
presented  in  Figure  4-3  for  20-  and  25-km  altitudes.  The  half-times  have  been 
computed  by  two  different  methods.  The  curves  marked  ambient  [O]  have  been  taken 
from  the  half-time  data  given  by  Johnston  and  are  based  on  initial  steady-state  [NO  ] 
levels.  The  curves  marked  steady-state  [O]  are  the  half-times  taken  directly  from 
the  ozone  versus  time  curves  of  Figures  4-1  and  4-2. 

The  half-times  for  steady-state  [O]  are  considerably  shorter  than  those  for  ambient 
[O]  , especially  at  higher  [NO  ] levels.  The  reason  for  this  is  that  the  initial  effect 

X 

of  adding  NO^  is  to  raise  the  steady-state  level  of  O atoms  by  virtue  of  the 
jh  [N02]  term.  This  new  steady-state  [ O]  level  is  reached  quite  rapidly  before 
any  depletion  of  O^  begins.  The  effect  of  this  higher  [ O]  is,  of  course,  to  speed 
up  the  catalytic  destruction  of  ozone  — hence,  the  shorter  half-times.  It  is  felt  that  the 
steady-state  [O]  half-times  are  the  appropriate  ones  for  the  present  discussion. 

It  should  also  be  noted  that  the  steady-state  [O]  half-times  presented  here  closely 
agree  with  those  calculated  by  Johnston's  method  if  the  initial  steady-state  [O]  level 
is  used  in  Johnston's  half-time  relationship.  This  may  seem  fortuitous  because  the 
half-time  relationship  given  by  Johnston  is  for  constant  [O]  and  constant  [NO] , 
while  these  quantities  vary  with  time  in  the  present  analysis.  However,  Johnston's 
half-time  depends  upon  the  product  [ O]  [NO]  which  does  tend  to  remain  constant 
during  the  initial  ozone  depletion  since  [O]  is  decreasing  while  [NO]  is  increasing. 
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Figure  4-3  O /NO  Catalytic  Half-Times 


4.2  PARTICULATE  FORMATION 


4.2.1  General  Considerations 

In  the  H,  C,  N,  O,  S,  h v system,  there  are  many  possible  reactions  that  can  lead  to 
formation  of  particulates.  These  are  reviewed  below. 

It  is  concluded  that  many  of  these  reactions  must  be  regarded  as  "unknown"  in  light  of 
present  knowledge.  Except  for  these  "unknown  areas"  it  is  not  anticipated  that 
heterogeneous  chemical  effects  attributable  to  the  presence  of  particulates  in  the  wake 
regime  will  significantly  modify  the  results  that  have  been  obtained  on  the  basis  of 
homogeneous  gas-phase  kinetics.  This  is  largely  a consequence  of  the  relatively  short 
residence  times  in  the  wake  — i.  e. , less  than  20  min. 

The  possible  significance  of  particulates  and  aerosols  in  stratospheric  chemistry  has 
been  discussed  briefly  by  Harrison  (Ref.  4-3).  The  two  types  of  mechanism  through 
which  particulates  can  affect  the  chemistry  are  (1)  adsorption,  which  simply  removes 
species  from  the  gas  phase  and  (2)  surface  catalysis  through  which  the  rates  of 
homogeneous  gas  processes  can  be  greatly  accelerated.  Harrison  concludes  that 
heterogeneous  chemistry  may  affect  the  stratospheric  kinetics  significantly  because 
characteristic  times  for  gas  diffusion  to  an  aerosol's  surface  may  be  comparable  to 
or  shorter  than  those  for  many  important  homogeneous  reactions.  He  further  concludes 
that  catalytic  effects  appear  to  be  potentially  more  significant  than  removal  by 
adsorption. 

Harrison's  remarks  are  based  on  ambient  stratospheric  concentrations  of  particulates. 
Therefore,  the  competition  between  homogeneous  and  heterogeneous  kinetics  may  be 
expected  to  become  even  more  important  if  additional  particulates  are  injected  in  the 
stratosphere  from  jet  regime  exhaust.  A start  has  been  made  at  treating  the  whole 
subject  of  stratospheric  aerosols  in  the  context  of  their  potential  climatic  modifications 
(Ref.  4-4). 
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The  present  discussion  focuses  on  photochemical  development  of  new  aerosols  within 
the  aircraft  wake.  Particulates  added  to  the  flow  at  the  engine  exit  plane,  such  as 
solid  carbon,  and  formation  of  new  aerosols  at  late  times  in  the  stratosphere  will  be 
covered  by  other  portions  of  the  CIAP  program. 

Formation  and  growth  of  new  aerosols  in  the  wake,  regime,  as  in  any  regime,  depends 
on  an  excess  vapor  pressure  condition,  i.e. , there  must  be  super  saturation  with 
respect  to  the  species  which  is  to  be  condensed.  In  other  words,  the  equilibrium 
vapor  pressure  of  species  in  condensation  with  the  condensed  aerosol  must  be  exceeded. 
There  is  a body  of  literature  that  refers  to  the  growth  of  aerosols  in  "unsaturated" 
conditions.  This  refers  to  the  fact  that  new  species  with  higher  molecular  weights  are 
formed  in  the  irradiated  gas,  and  that  the  vapor  pressure  of  these  species  with  respect 
to  the  various  groups  it  contains  is  quite  low.  Thus,  in  an  environment  that  is  highly 
unsaturated  with  regard  to  water  — i.  e. , when  the  vapor  pressure  of  water  is  far  below 
the  saturation  vapor  pressure  — photochemical  formation  of  new  species  with  low  vapor 
pressure  may  still  be  possible.  In  this  case,  the  growth  of  aerosols  from  the  vapor 
phase  is  still  possible,  skin  to  photochemical  synthesis,  and  this  type  of  process  is 
sometimes  referred  to  as  particle  formation  under  unsaturated  conditions  (Ref.  4-5). 

For  consideration  of  normal  aerosol  formation  under  saturated  or  supersaturated 
conditions,  condensation  of  water  vapor  from  the  engine  exhaust  is  the  most  obvious 
potential  process. 

Here  ther  are  ample  condensation  nuclei  to  serve  as  condensation  sites  if  condensation 
favoring  conditions  exist.  These  nuclei  include  ambient  particulate  species  taken  into 
the  engine,  soot  and  ions  (exhaust  plasma  levels  correspond  to  approximately 
10  ions/cm  ).  On  the  other  hand,  the  rates  of  condensation  in  relation  to  the  time 
scales  of  interest  are  highly  uncertain.  It  is  noteworthy  that  the  nucleation  and  growth 
of  water  particles  to  visible  size  within  the  wake  regime  is  unlikely  because  visible 
contrails  are  almost  never  observed  in  stratospheric  flight. 
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A simple  exercise  will  demonstrate  that  a saturated  condition  for  H O vapor  exists 

& 

only  briefly  in  the  wake.  Assuming  complete  homogeneity,  a wake  energy  balance 
yields  the  following  relation  between  water  concentration  and  temperature : 


[H2OJw  “ ’ 


T ) 
a 


- T ) 
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where:  subscript  w denotes  the  wake,  subscript  a the  ambient,  and  subscript  e the 
engine  exhaust  conditions;  Z is  the  ratio  of  the  specific  heat  of  the  exhaust  products 
to  that  of  air.  Because  of  its  very  low  2-ppm  level  (Ref.  4-6)  the  ambient  stratospheric 
water  vapor  can  be  neglected. 

The  above  relationship  has  been  plotted  in  Figure  4-4  for  the  GE-4  engine  with  after- 

5 

burning  — T = 950 °K,  [H  O]  = 10  ppm,  and  Z = 1.  The  curve  for  cruise  operation 

0 Lt  0 

would  not  be  greatly  different  because  the  effect  of  the  lower  water  concentration  would 
be  cancelled  by  a lower  exhaust  temperature.  Also  included  in  Figure  4-4  are  curves 
representing  saturated  conditions  at  15,  20,  and  25  km  altitude. 


The  figure  shows  there  is  no  condensation  to  be  expected  at  25  km,  but  there  is  a brief 
period  when  slightly  supersaturated  conditions  exist  in  a 20-km  wake.  This  period 
encompasses  the  region  between  2300  ppm  and  550  ppm  water  concentration,  which 
implies  the  range  of  air  dilution  factors  from  45  to  180.  At  15  km  the  wake  becomes 
saturated  at  a dilution  factor  of  25,  and  the  condition  continues  until  a dilution  of 
550  is  reached.  The  time  scale  for  the  wake  to  grow  through  this  range  of  dilution 
is  of  the  order  of  1 to  10  sec.  Therefore,  at  lower  altitudes  the  extent  of  water 
condensation  in  the  wake  will  depend  on  the  rates  of  nucleation  and  droplet  growth. 

There  should  be  ample  nuclei  in  the  engine  exhaust.  The  lack  of  contrail  formation 
in  the  stratosphere  must  therefore  be  due  to  the  failure  of  the  particles  to  grow  to 
visible  size  during  the  brief  period  of  saturated  conditions. 

In  reviewing  the  many  possible  reactions  capable  of  forming  particulates,  the  objective 
has  been  to  identify  those  reactions  in  the  wake  of  a GE-4  engine  operating  at  20-km 
altitude  which,  in  the  process  of  particulate  formation,  will  lead  within  20  minutes 
to  significant  changes  (e.g.  , 10  percent)  in  the  appropriate  exhaust  species  composition. 
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Figure  4-4  Water  Condensation  Regimes  in  the  Wake 
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Additionally,  reactions  with  half  times  of  up  to  about  a week  are  identified  to  aid 
in  analysis  of  the  diffusion  regime.  For  initiation  of  the  latter,  a dilution  of  the 
exhaust  by  a factor  of  3,000  is  assumed.  To  provide  a basis  for  kinetic  analysis,  an 
estimate  of  the  engine  exhaust  characteristics  must  be  made.  Additional  information 
is  required  relative  to  the  ambient  chemical  species,  particulates,  and  light  in- 
tensity and  spectrum. 

From  an  analysis  of  aircraft  engines  without  afterburners,  it  has  been  found  that  the 
concentrations  of  unburned  hydrocarbons  are  extremely  low  (Ref.  4-7).  With  after- 
burners in  operation,  however,  the  hydrocarbon  concentration  generally  increases 
drastically  (Refs.  4-8  to 4-10).  The  concentration  of  hydrocarbons  may  be  as  much 
as  10  percent  of  the  carbon  monoxide  formed.  The  maximum  concentration  of  sulfur 
in  the  JP-5  fuel  is  specified  as  0.1  weight  percent.  With  an  overall  equivalence  ratio 
of  0.8,  a maximum  sulfur  dioxide  content  of  about  100  ppm  is  calculated.  Other  exhaust 
species  concentrations  are  taken  from  the  previous  sections.  These  are  summarized 
in  Table  4-2.  It  should  be  noted  that  the  nitric  oxide  and  SO^  concentrations  are 
considerably  higher  than  those  employed  in  atmospheric  pollution  studies  which  have 
been  the  principal  source  of  photochemical  kinetic  data  for  this  analysis. 

Table  4-2 


ENGINE  EXHAUST  CONCENTRATIONS 


Concentration 

Component 

Mol  Fraction 

PPM  by  Volume 

n2 

0.747 

h2o 

0. 102 

C02 

0.093 

o2 

0.  055 

CO 

0.  001 

h,h2,no,o 

400 

OH 

200 

Hydrocarbon 

250 

so2 

<100 
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However,  in  those  cases  where  information  is  available  for  relatively  high  concentrations, 
the  results  are  in  general  accord  with  those  obtained  at  lower  concentrations. 

Approximate  concentrations  (Ref.  4-11)  of  principal  ambient  species  are  summarized  in 
Table  4-3.  The  estimated  concentration  of  ambient  oxygen  atoms  is  from  Johnston 
(Ref.  4-1). 

The  ambient  particulate  characteristic  at  20  km  are  described  by  Junge  (see  Ref.  4-12). 

3 

A number  density  of  0.1  particle /cm  found  above  midwestern  United  States  probably 
represents  a world-wide  condition  with  little  variation  in  time  and  latitude. 


Table  4-3 

AMBIENT  SPECIES  AT  20  KM 


Component 

3 

Log  Concentration,  Molecule/Cm 

N2 

18 

°2 

17.  3 

CH 

13 

°3 

12 

0 

6 

A majority  of  the  particles  at  this  altitude  fall  in  the  0. 1-1  pm  range.  Molecular  sized 
particles  will  be  attached  very  rapidly  to  the  established  background  population. 
Particles  smaller  than  0.  01  pm  radius  will  be  eliminated  by  coagulation  in  times 
that  are  short  compared  to  established  residence  times  for  stratospheric  air.  They 
may,  however,  materialize  as  a result  of  injection  by  the  engine  exhaust.  Larger 
particles  will  be  increasingly  subject  to  the  effect  of  gravity.  Specific  analysis  of 
particles  recovered  from  20  km  has  revealed  a particle  radius  variation  from  0.  23 
to  1.48  pm  (see  Ref.  4-13).  Silicon,  iron,  and  sulfur  were  detected  with  the  sulfur 
content  of  particles  indicated  to  vary  from  5 to  71  percent. 
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Since  many  of  the  reactions  which  might  lead  to  the  formation  of  particulates  are 
photochemical  in  nature,  the  light  intensity  at  20-km  altitude  is  relevant  to  the 
problem.  Attenuation  of  the  solar  spectrum  will  result  principally  from  Rayleigh 
scattering,  aerosol  scattering,  and  ozone  absorption.  The  relative  importance  of 
each  will  vary  as  a function  of  the  wave  length.  Tabulations  of  extinction  coefficients 
and  optical  thickness  (Ref.  4-14)  plus  subsequent  modifications  (Ref.  4-15)  were  used 
as  a basis  for  the  relative  intensities  given  in  Table  4-4.  The  intensities  are  ex- 
pressed as  relative  amounts  reaching  20  km  from  outer  space  and  relative  amounts 
reaching  sea  level  from  20  km.  Attenuation  coefficients  are  included  to  show  the 
relative  importance  of  each  of  the  mechanisms. 

Laboratory  light  intensities,  when  given  in  the  aerosol  literature,  are  usually  related 
to  sunlight  at  sea  level.  The  increased  intensities  at  20  km  must  therefore  be  taken 

o o 

into  account.  From  Table  4-4,  this  factor  may  vary  from  2 at  4,  000  A to  16  at  3,  000  A. 
Short  wave  lengths  are  absorbed  by  ozone. 

The  effect  of  increased  light  intensity  is  often  roughly  compensated  for  by  the  effect  of 

reduced  temperature  at  20  km.  For  example,  in  the  laboratory  photo -oxidation  of 

nitric  oxide  in  the  presence  of  trans -2 -butene,  Altshuller  and  Bufalini  found  that  an 

increase  of  20°C  caused  the  same  increase  in  the  rate  of  oxidation  as  a 2-fold  increase 

in  light  intensity  (see  Ref.  4-16).  Since  the  temperature  at  20  km  is  roughly  80°C  below 

4 

laboratory  temperatures,  a reduction  in  rate  of  2 , or  16,  is  indicated  for  a process 

o 

of  the  same  activation  energy.  At  wave  lengths  above  3000  A , therefore,  the  net  effect 
of  the  increased  altitude  may  actually  be  a reduction  in  the  rate  observed  under  labor- 
atory conditions. 

The  reactions  of  the  engine  exhaust  species  with  each  other  and  with  ambient  species 
to  form  particulates  may  be  grouped  as  sulfur  dioxide  reactions,  nitrogen  oxide 
reactions,  and  combinations  of  both.  Finally,  consideration  is  given  to  heterogeneous 
reactions  between  particulates  and  gaseous  species. 
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Table  4-4 


ATMOSPHERIC  ATTENUATION 


Wave  length 
A 

3 -1 

Attenuation  Coefficients,  10  km 

Relative  Intensities 

Rayleigh 

Aerosol 

Ozone 

I km/lon  km 
o 20 

*20  km/I~ 

6000 

0.6 

1.4 

2.2 

0.73 

0.96 

4000 

3. 1 

1.9 

0 

0.  52 

0.97 

3400 

6.2 

2.3 

1. 1 

0.35 

0.94 

3200 

8.  0 

2.4 

14.7 

0.25 

0.78 

3000 

11 

2.  5 

166 

0.  06 

0.11 

2800 

14 

2.6 

1740 

3.2  x 10~7 

1. 12  x 10'10 

4.2.2  Sulfur  Dioxide  Reactions 

As  discussed  previously  in  section  3.4.  13  the  oxidation  of  SOQ  by  O atoms  in  the  wake 
is  negligible.  In  the  rapidly  cooling  mixing  region,  function  of  S03  is  forestalled  by 
the  faster  removal  of  O atoms  by  O + O^  + M — O^  + M . 

S02  Photolysis  (dry).  Widely  varying  rates  of  photolysis  of  S02  in  air  have  been 

published.  Hall  found  that  S00  (56  to  230  mm)  and  09  (5  to  200  mm)  reacted  in  sun- 

" —6  —1 

light  with  a first-order  rate  constant  of  8 x 10  min  (see  Ref.  4-17).  Sethi  et  al. 

-3  ° 

claimed  a quantum  yield  of  2 x 10  over  the  wave  length  range  of  2800  to  4200  A (see 

Ref.  4-18).  Renzetti  and  Doyle  (Ref.  4-19)  found  a photolysis  rate  of  0.45  percent/ 
min.  Urone  et  al.  (Ref.  4-20)  cite  a rate  0. 1 percent/hr  in  the  equivalent  of  a noon- 
day sun.  In  the  0. 14  to  1.  0 ppm  sulfur -dioxide  range,  Cox  and  Penkett  (Ref.  4-21 ) 

-3  -1 

found  a first-order  rate  constant  of  2 x 10  hr  , essentially  independent  of  SC>2 
concentration.  Gerhard  and  Johnstone  (Ref.  4-22)  worked  in  the  2950  to  3650  A 
range  with  SC>2  (5  to  30  ppm)  in  air  and  found  a first-order  conversion  of  0. 1 to  0.2 
percent/hr  in  intense  sunlight  (3x  sun).  The  short-term  (20  min)  conversion  to  S03 
will  be  insignificant  based  on  most  of  the  above  rates.  However,  in  one  week  the 
conversion  should  exceed  10  percent. 
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SC>2  Photolysis  (moist).  The  effect  of  moisture  on  the  oxidation  reaction  rate  is  debat- 
able. Water  vapor  was  found  to  have  no  effect  on  the  rate  of  less  than  0. 1 percent/hr 
except  at  100  percent  relative  humidity,  where  the  rate  was  slightly  greater  (Ref.  4-22). 

Gartrell  et  al.  found  (Ref.  4-23)  that  the  oxidation  of  S09  increased  from  1/6  percent/ 
min  to  about  2 percent/min  in  passing  from  low  to  high  humidity. 

The  conclusion  is  drawn  that  the  somewhat  moist  S02  oxidation  probably  falls  in  the 
same  category  as  the  dry  reaction  — of  minor  significance  for  20  min,  but  significant 
over  a 1-wk  period.  The  significant  aspect  is  that  the  subsequent  formation  of  aerosol 
is  always  fast  compared  to  the  oxidation  rate.  Thus,  SOQ  + H90  — -H9S0  is  fast  — 
i.e.  , H2SC>4(£)  droplet  formation  is  much  faster  than  the  2%/min  SOg  maximum  for- 
mation rate. 

SO2  + Ozone.  This  reaction  is  slower  than  the  irradiation  oxidation  (Ref.  4-24)  and  is 
also  considered  to  be  too  slow  to  be  significant  at  atmospheric  concentrations  (see  Ref. 
4-25).  It  therefore  would  probably  not  be  significant  in  20  min  in  the  stratosphere.  It 
is  possible,  however,  that  a significant  amount  could  react  in  a week,  because  of  the 
much  higher  ambient  concentration  of  ozone  relative  to  oxygen  atoms. 

— 32  6 2 

SC>2  + O + Third  Body.  A third-body  reaction  rate  constant  of  8 x 10  “ cm  molecule 

r~l -4 

sec  was  published  by  Cadle  (see  Ref.  4-26).  This  leads  to  an  initial  rate  of  6 x 10 

percent/min  with  oxygen  atoms  in  the  stratosphere,  but  the  rate  is  eight  orders  of 

magnitude  greater  in  the  exhaust  due  to  the  vastly  higher  concentration  of  O-atoms. 

15  6 ~2  — 1 

A rate  constant  of  1.  5 x 10  cm  mole  sec  at  784 °K  has  been  calculated  by  Levy 
et  al.  (see  Ref.  4-27).  This  is  actually  somewhat  smaller  than  other  values  reported 
in  the  same  paragraph  for  the  reaction  at  300 °K  and  an  order  of  magnitude  lower  than 
the  estimate  of  Webster  and  Walsh  (Ref.  4-28)  for  784°K.  With  respect  to  the  reaction 
rate  near  the  exhaust  jet,  high  O-atom  levels  and  oxidation  rates  have  been  discussed 
previously. 

SO2  + Ozone  + Olefins.  High  concentrations  of  ozone  and  olefins  have  been  observed 
(Ref.  4-29)  to  produce  a dense  smoke,  particularly  with  cyclic  olefins,  and  the 
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addition  of  S02  did  not  significantly  affect  the  rate  of  aerosol  formation  or  the  amount 

obtained.  However,  the  aerosol  appeared  rapidly  and  the  maximum  amount  formed 

within  5 min.  Leighton  has  tabulated  (Ref.  4-30)  second-order  rate  constants  for 

ozone-olefin  reactions.  Most  are  in  the  2000  to  20,000  1 mole-1  sec-1  range  with  a 
5 -1  -1 

maximum  of  10  1 mole  sec  for  trans-2-butene.  If  the  latter  rate  is  applied 

_y 

to  the  exhaust  (Oq  = 1.6  x 10  mole/1  and  assuming  10  percent  of  the  exhaust 
^ -8 

hydrocarbon  is  trans -2 -butene  = 1.8  x 10  mole/1)  it  yields  an  initial  maximum 
rate  of  1.  0 percent/m  in.  Both  the  rate  constant  and  olefin  concentration  are  likely 
to  be  much  lower.  In  addition,  a volume  dilution  factor  of  about  3000  in  the  wake  is 
anticipated. 

Kinetic  data  are  also  given  by  Cox  and  Penkett  (Ref.  4-25)  for  low  concentrations. 

It  is  suggested  that  the  mechanism  involves  the  formation  of  a peroxide  zwitterion 
during  the  ozone-olefin  reaction.  Although  the  S09  did  not  have  any  significant  effect 
on  the  rate  of  disappearance  of  ozone  or  hydrocarbon,  large  amounts  of  aerosol  were 
produced  when  S02  was  added  to  the  mixture.  At  ozone  and  olefin  concentrations 
of  0.  05  ppm,  the  oxidation  rate  of  0. 1 ppm  S02  is  estimated  to  be  about  3 percent/hr 
for  cis-2-pentene  and  0.4  percent/hr  for  propene. 

Based  on  the  foregoing  discussion,  it  is  concluded  that  this  reaction  will  not  be 
significant  within  20  min,  but  that  the  olefin  concentration  may  be  reduced  by  reaction 
within  a 1-wk  period.  However,  if  the  rate  of  S09  removal  is  proportional  to  the 
olefin  concentration,  dilution  by  more  than  three  orders  of  magnitude  will  reduce 
the  rate  to  insignificance  for  a 1-wk  period. 

o 

Photolytic  S09  + Organic.  Irradiation  at  2400  to  3340  A of  saturated  and  unsaturated 
hydrocarbons  in  the  presence  of  S0o  has  been  found  by  Dainton  and  Ivin  to  produce 
sulfinic  acids  (see  Ref.  4-31).  To  the  extent  that  the  aerosol  is  organic  it  may  ex- 
perience instability  by  evaporation  or  oxidation  (see  Ref.  4-32).  Harkins  and  Nicksic 
(Ref.  4-33 ) found  that  S09  consistently  formed  an  aerosol  in  contact  with  hydrocarbons, 
that  heavy  hydrocarbons  produced  aerosols  without  irradiation,  and  that  very  little 
S09  actually  oxidized.  Another  interesting  observation  was  that  an  inverse  correlation 
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exists  between  chamber  temperature  and  aerosol  formation.  Johnston  and  Dev  Jain 
(Ref.  4-34)  irradiated  50  mm  each  of  n-butane  and  S09  in  400  mm  air  and  produced 
heavy  fogs  within  2 min. 

Although  aerosols  containing  sulfinic  acids  are  produced,  the  rate  of  aerosol  formation 
is  negligible  when  the  reactants  are  at  concentrations  of  100  ppm  or  less  (see  Ref. 

4-35).  Kopcynski  and  Altshuller  also  found  that  saturated  hydrocarbons  did  not  con- 
tribute significantly  to  aerosol  production  when  compared  to  olefins.  They  concluded 
that  the  rate  of  photochemical  reaction  of  paraffins  or  olefins  with  SO^  was  very  slow 
in  the  100  to  1000  ppm  range.  It  is  therefore  concluded  that  the  reaction  will  be  insig- 
nificant in  the  wake  regime.  There  is  insufficient  information  to  indicate  its  possible 
importance  in  the  diffusion  and  transport  regime. 

4.  2.  3 Nitrogen  Oxide  Reactions 

The  exhaust  is  estimated  to  contain  300  ppm  NO  with  almost  no  N02  present  initially. 

The  relative  amounts  of  NO  and  N09  are  very  significant  in  the  reaction  rates  of  organic 
compounds.  This  influence  is  believed  to  be  dependent  on  the  fact  that  NO  is  more 
rapidly  oxidized  than  hydrocarbons  by  reactive  oxygen  species.  In  general,  the 
presence  of  NO  retards  the  organic  oxidation  processes. 

Photolytic  NO  + Organic.  Methane,  ethane,  propane,  and  2,  2 -dim ethyl  propane  are 
relatively  unreactive  in  the  photo -oxidation  of  NO.  However,  NO  photo -oxidation  rates 
appear  to  correlate  with  aerosol  formation  (see  Ref.  4-36).  Di-substituted  internal 
olefins  (2,  3 dimethyl-2-enes)  have  the  highest  NO  photo -oxidation  rate.  For  example, 

2,  3 dimethyl -2 -butene  (1  ppm)  and  NO  (0.4  ppm)  provided  a photo -oxidation  rate  of 
0.  059  ppm  min  \ The  majority  of  some  40  other  olefins  evaluated  by  Glasson  and 
Tuesday  reacted  in  the  range  of  0.  002  to  0.012  ppm  min  . Stephens  and  Schuck 
(Ref.  4-37)  irradiated  NO  plus  auto  exhaust  and  obtained  aerosols  after  30  to  60  minutes. 

It  is  concluded  that  certain  olefins  could  react  significantly  within  a 20-min  period; 
most  will  not.  Since  the  special  olefins  will  constitute  only  a very  minor  part  of  the 
total  exhaust  hydrocarbons,  the  total  hydrocarbon  content  will  not  change  significantly 
during  this  period. 
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Photolytic  NO  + Olefin  + CO.  According  to  Westberg  et  al.  , the  presence  of  CO 
markedly  accelerates  the  disappearance  of  olefin,  the  conversion  of  NO  to  N09,  and 
the  appearance  of  ozone  (Ref.  4-38).  However,  from  the  published  curves,  the 
presence  of  100  ppm  CO  had  little  effect  in  20  min  on  the  photo -oxidation  of  NO  (1.  5 ppm) 
in  the  presence  of  isobutene  (3  ppm)  and  70  percent  relative  humidity.  The  mechanism 
is  dependent  upon  a hydroxyl  carrier,  as  follows: 

OH  + CO  — C02  + H 

H + 02  + (M)  — OOH  + (M) 

OOH  + NO  — - OH  + N02 

The  more  rapid  reactions  of  ozone  and  O-atoms  with  NO  would  probably  override 
this  reaction. 

The  reaction  does  not  appear  to  be  of  concern  in  the  20-min  picture,  since  dilution 
by  3000  will  reduce  the  CO  concentration  from  1000  to  less  than  1 ppm.  From  the 
data  in  Ref.  4-38,  it  appears  that  changes  in  the  olefin  and  NO  concentrations  could 
be  significant  in  the  diffusion  regime  although  in  this  case  the  CO  would  tend  to 
reduce  the  amount  of  N02  formed. 

Photolytic  NO2  + Organic.  The  photolysis  of  N02 -hydrocarbon  mixtures  in  air  has 
been  examined  by  Schuck  and  Doyle  (see  Ref.  4-39  ).  The  maximum  rate  constant 
for  3-methyl  pentane  (6ppm)  and  N09  (2ppm)  corresponded  to  0. 15  percent/min. 

The  olefins  produced  much  greater  rates.  For  example,  trans -2 -butene  (~  3 ppm) 
and  N09  (~  1 ppm)  yielded  a rate  corresponding  to  13.  6 percent/min.  Other  olefins 
were  in  the  1.5  to  3.5  percent/min  category. 

From  these  data,  it  is  concluded  that  the  total  hydrocarbons  will  not  be  significantly 
reduced  in  the  20-min  period  when  the  3000-fold  volume  dilution  is  taken  into  account. 
However,  the  occurrence  of  significant  photolytic  reactions  with  olefins,  even  with 
dilution,  is  conceivable  over  a one-week  period. 
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Inorganic  Nitrates  and  Nitrites.  Nitric  acid  is  not  sufficiently  hygroscopic  to  form 
aerosols  in  the  stratosphere  (see  Ref.  4-40).  There  are  many  possible  reactions  for 
formation  and  destruction  of  nitrous  acid,  for  example: 


M + NO  + OH 
N02  + OH 

no2  + h2o 

etc. 

H + HN02 
OH  + HN02 
etc. 


— HN02  + M 

— hno2  + O 

— HN02  + OH 

— h2  + N02 

— h2o  + N02 


The  first  reaction  in  each  set  utilizes  readily  available  species,  i.e.  , NO,  OH,  and 
H from  the  exhaust  jet. 

Rate  data  are  not  available  for  some  of  these  reactions  (Ref.  4-41),  but  it  seems 
unlikely  that  aerosol  formation  will  result,  even  if  HN02  is  present.  In  moist  air, 
nuclei  may  form  by  hydration  of  HN02  or  HN03  , but  the  aerosol  formation  rate 
is  very  slow  according  to  Vohra  (see  Ref.  4-42). 

Ammonium  nitrate  could  easily  form  particulates,  for  example: 

nh3  + hno3  — nh4no3 

However,  as  discussed  in  Section  3,  very  little  HNC>3  is  anticipated  in  the  exhaust 
and  wake,  and  there  is  no  evidence  that  significant  quantities  of  NH3  may  be  expected. 

In  general,  inorganic  nitrate  particulate  formation  in  the  wake  is  not  anticipated. 

4.2.4  Combinations  of  Nitrogen  and  Sulfur  Oxide  Reactions 

SO2-NO.  The  combination  of  S02  and  NO  produces  no  aerosol  according  to  Prager  et  al. 
(see  Ref.  4-43 ). 
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Photolytic  SO2  + NO  + Organic.  Haagen-Smit  (Ref.  4-29)  found  that  the  irradiation 
of  cyclohexene,  NO,  and  SC>2  did  not  produce  any  aerosol  in  the  first  half  hour. 

In  the  irradiation  of  propylene  (2  ppm),  NO  (0.5  ppm),  and  S09  (1.2  ppm)  Altshuller 
et  al.  (Ref.  4-44)  found  that  the  rate  of  consumption  of  propylene  was  independent 
of  the  amount  of  S02.  With  essentially  the  same  concentration  of  reactants,  Cox  and 
Penkett  obtained  a 7 to  10  percent  conversion  of  S09  to  sulfate  aerosol  in  a 2 -hr 
irradiation  (see  Ref.  4-21).  With  cis-2-pentene  (0.98  ppm),  NO  (0.73  ppm),  and 
S02  (0.  07  to  0.25  ppm),  the  same  experimenters  obtained  an  S02  reaction  rate  of 
0.  64  percent/min.  In  the  latter  case,  the  chamber  intensity  was  only  one-fourth 
that  of  sunlight.  Aerosol  appeared  after  100  min  (Ref.  4-45)  in  the  system  1-butene 
(4  ppm),  NO  (1  ppm),  and  S02  (0.72  ppm).  Ozone  does  not  appear  until  almost  all  the 
NO  is  converted  to  N02  (see  Ref.  4-43). 

Since  the  induction  period  for  particulate  formation  with  5 ppm  NO  approaches  20  min 
a significant  concentration  change  is  not  anticipated  within  this  time  interval.  Over 
the  period  of  1 wk  however,  certain  olefins  may  undergo  significant  changes  although 
a 3000-fold  dilution  will  tend  to  offset  this  possibility.  Except  for  independence  of  SO2 
concentration,  not  enough  details  of  the  mechanism  are  known  to  predict  the  effect 
on  NO  and  S02  concentrations. 

Photolytic  S02  + N09.  The  irradiation  of  S02  and  N02  has  been  studied  by  Jaffee  and 

Klein  (see  Ref.  4-46).  They  obtained  an  overall  rate  constant  which  was  a function 

of  the  third  body  concentration.  For  the  conditions  in  the  exhaust  and  in  the  ambient 

10  2 -2 

stratosphere,  k = 1.4  x 10  1 mole  sec.  The  overall  rate  at  the  exit  plane  of  the 

exhaust  is 

Rate  = 1.4  x 1010  [O]  [S02]  [M]  60  = 5. 1 x 10~5  moll”1  min"1 
4 

or  1.8  x 10  percent/min,  a rate  reasonably  close  to  the  earlier  estimates  for  the 
three-body  SOg-O  reaction.  Gerhard  and  Johnstone  (Ref.  4-22)  irradiated  S02  (5  to 
30  ppm)  and  N09  (2  ppm),  and  found  no  effect  of  salt  or  moisture  on  the  rate. 
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Renzetti  and  Doyle  found  (Ref.  4-19 ) that  N02  enhanced  aerosol  formation  whereas 
NO  suppressed  it. 

This  reaction  may  be  presumed  similar  to  the  three-body  S09-0  reaction  previously 
discussed. 

Photolytic  S09  + N09  + Organic.  The  addition  of  12  ppm  n-hexane  to  S02  (16  ppm), 
N02  (15  to  410  ppm),  and  air  containing  0 to  100  percent  relative  humidity  showed  that 
the  reaction  rate  was  independent  of  moisture  content  except  at  100  percent  humidity 
(see  Ref.  4-20).  The  highest  rate  observed,  0.4  percent/min  was  found  at  the  highest 
moisture  and  N02  levels.  In  most  cases  less  than  10  percent  of  the  S02  reacted  over 
extended  periods  of  time.  Aerosol  was  not  produced  (Ref.  4-43)  with  n-butane  or 
2 -methyl  pentane. 

A considerable  amount  of  work  has  been  devoted  to  examining  the  increased  aerosol 
production  rate  obtained  by  adding  olefins  to  the  N09-S09  system.  All  types  of  olefins 
produced  particulates.  The  induction  period  for  conversion  of  NO  to  N02  was  again 
noted  during  the  irradiation  of  engine  exhaust  (see  Ref.  4-39).  Diolefins  were  found 
(Ref.  4-43)  to  increase  the  reaction  rate  less  than  straight- chain  mono-olefins.  Cyclic 
olefins  formed  particulates  more  rapidly  than  straight-chain  olefins  (see  Ref.  4-29). 
The  reaction  was  completed  in  0.  5 hr  or  less,  according  to  Haagen-Smit. 

The  photolysis  of  1 -hexene  (1.  5 ppm),  with  1 ppm  each  of  S02  and  N02  in  air,  removed 
S02  at  up  to  15  percent/hr  (see  Ref.  4-17).  Endow  et  al.  found  that  sulfuric  acid  was 
a major  component  of  the  resulting  aerosols  and  that  aerosols  formed  at  10  to  20 
percent  relative  humidity  contained  more  inorganic  nitrogen  than  those  formed  at 
50  percent  relative  humidity  (see  Ref.  4-47).  Goetz  and  Pueschel  (Ref.  4-48) 
determined  that  the  growth  rate  of  the  aerosol  was  proportional  to  the  hydrocarbon 
concentration.  They  also  found  that  the  sequence  of  addition  of  reactants  influenced  the 
aerosol.  When  Stevenson  et  al.  added  less  than  1 ppm  SC>2  to  systems  containing 
10  to  20  ppm  unsaturated  hydrocarbons  and  3 ppm  N02»  the  photolytic  aerosol  pro- 
duction increased  (see  Ref.  4-49). 
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It  is  concluded  that  the  occurrence  of  a significant  degree  of  hydrocarbon  reaction 
with  SO?  and  N02  within  20  min  would  not  be  expected;  further,  little  or  no  N02  is 
initially  present.  The  importance  of  the  reaction  will  depend  on  the  speed  with  which 
NO  is  converted  to  N09.  Certain  olefins  may  react  substantially  over  a 1-wk  period. 
This  indication  is  based  upon  the  observed  rates  in  the  S09-free  system  and  the 
increased  aerosol  production  observed  when  SC>2  was  added  to  the  system.  Here  again, 
the  overall  hydrocarbon  concentration  would  not  decrease  significantly,  since  a major 
portion  of  the  composition  is  saturated. 

4.  2.  5 Heterogeneous  Reactions 

The  potential  effect  of  recently-formed  or  existing  particulates  on  chemical  reactions 
forming  aerosols  may  be  explored.  First,  it  has  already  been  shown  that  the  ambient 
concentration  of  particles  at  20  km  altitude  is  extremely  low  (100  particles/liter)  and 
is  therefore  not  likely  to  significantly  affect  the  reactions  of  the  exhaust  species.  The 
potentially  catalytic  nuclei  which  result  from  the  exhaust  itself  may  be  classified  as 
aqueous  sulfate  systems  or  basic  oxides. 

Aqueous  Sulfate  Systems.  The  rate  of  oxidation  of  S09  in  30  percent  sulfuric  acid  (Ref. 
4-22)  is  said  to  be  negligible  compared  to  the  rate  of  photochemical  oxidation.  In 
the  presence  of  ammonia,  however,  Van  der  Hemel  and  Mason  (ammonia-S02-water- 
air  system)  found  (Ref.  4-50)  that  oxidation  of  S09  in  cloud  droplets  may  be  more 
efficient  than  the  photochemical  conversion. 

Scott  and  Hobbs  (Ref.  4-51)  observed  that  the  concentration  of  sulfate  in  a droplet 
after  24  hr  is  about  two  orders  of  magnitude  greater  when  ammonia  is  present  than 
when  it  is  absent.  By  assuming  that  the  rate  of  formation  of  sulfate  is  proportional 
to  the  sulfite  ion  concentration  and  considering  the  several  equilibria  involved,  the 
following  rate  expression  is  developed: 


= K [SOs  ] 


where  H+  is  related  to  P 


NH 


3 


4-26 


Information  on  stratospheric  ammonia  was  not  found.  An  average  sea-level  atmospheric 

_9 

concentration  of  7 x 10  atm,  together  with  the  estimated  exhaust  S00  partial  pres- 

-6  ^ 
sure  of  4.7  x 10  atm,  would  yield  an  initial  conversion  rate  of  1580  percent/min. 

Certain  salts,  notably  those  of  iron,  copper,  and  manganese,  are  known  to  catalyze 
the  conversion  of  S02  to  sulfate.  Johnstone  and  Coughanowr  (Ref.  4-52)  estimated 
that  1 ppm  S09  oxidized  at  1 percent/min.  Manganese  sulfate  produced  the  highest 
rate  among  the  salts  tested.  Johnstone  and  Moll  found,  when  manganese  or  ferrous 
sulfates  were  present  in  the  250-ppm  S02  - moist  air  system,  that  about  1 percent/ 
min  of  sulfuric  acid  formed  in  fogs  (Ref.  4-53). 

It  has  not  been  established  to  what  extent  ammonia,  iron,  or  manganese  are  present 
either  in  the  engine  exhaust  or  in  the  ambient  stratosphere.  At  this  time,  it  can  only 
be  stated  that  if  any  or  all  are  present,  significant  reductions  in  S02  concentration 
could  occur  in  the  dispersion  regime.  In  the  wake  region,  diffusion  to  the  catalytic 
nuclei  would  probably  limit  the  SO2  concentration  reduction  rates. 

Basic  Oxides.  Although  sodium  chloride  showed  relatively  little  effectiveness  in 

oxidizing  S09  (Refs.  4-20,  4-53),  certain  basic  oxides  have  demonstrated  great 

efficiency  in  scavenging  S09.  Forty  percent  of  the  S09  (1  ppm)  in  air  (46  percent 

" 8 z 

relative  humidity)  was  quickly  absorbed  by  1.8  x 10  nuclei  per  liter.  Oxidized 
lead  aerosols  removed  S02  completely  in  less  than  5 min  according  to  Smith  (See 
Ref.  4-54).  Dried  hydrous  iron-oxide  particles  were  successful  in  removing  S09 
within  seconds  (Ref.  4-20),  but  this  result  was  obtained  with  high  powder/SO£  ratios. 
Oxides  of  aluminum,  calcium,  and  lead  also  yielded  high  absorption  rates. 

It  is  unlikely  that  these  oxides  would  be  found  in  either  the  wake  or  dispersion  areas  in 
sufficient  concentration  to  significantly  alter  the  S02  concentration. 
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Section  5 

RADIATIVE  COOLING  AND  RADIATIVE  INDUCED  BUOYANCY 


In  this  section  we  address  ourselves  to  two  questions.  One,  does  the  energy  loss  rate 
due  to  radiative  emission  of  the  hot  exhaust  gases  at  the  exit  temperature  have  a 
significant  effect  on  the  initial  wake  cooling?  Second,  does  radiative  heating  or  cooling 
in  the  long-time  dispersion  wake  result  in  a slight  temperature  change  (of  a few  degrees) 
which  would  induce  buoyancy  effects  ? 

In  answering  the  first  question,  the  proper  perspective  is  obtained  by  recognizing  that 
the  initial  wake  mixing  with  cool  air  results  in  a 500  °K  temperature  drop  in  a few 
milliseconds.  At  typical  exit  temperatures  around  2000 °K  the  primary  source  of 
radiative  emission  is  the  intense  4.3-/i  band  of  CO^.  An  upper  estimate  to  the  radiative 
cooling  rate  was  determined  by  assuming  an  optically  thin  model,  i.e.  , one  in  which 
all  energy  emitted  by  the  CO  molecules  escapes  from  the  wake. 

The  energy  loss  per  unit  mass  of  C09  is 

(Q  ) = 4tt  ol  B (T)  / ^ \ 

' rad  CQ  band  c oq  ^gm  CQ  sec  / 

2 2 

when  o^anc[  is  the  integrated  band  absorption  for  the  4.3 -p  C02  band,  and  B^  (T)  is 
the  Planck  function  evaluated  at  the  band  center  frequency  (cuq  = 2,350  cm  1)  and  at  the 
exhaust  gas  temperature.  The  energy  loss  due  to  C09  is  shared  by  all  of  the  exhaust 
species.  The  energy  loss  per  unit  mass  of  exhaust  gases  is 

^rad^  , CCO  ^rad^co  (gm-sec) 

exhaust  2 C02 
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where  C is  the  mass  fraction  of  COo  species  for  which  we  assumed  a value  of 

ou2 

0.25.  The  temperature  rate  of  change  is  simply  (Q  ,)  divided  by  the  specific 

^aa  exhaust 

heat  of  the  high-temperature  exhaust  gases.  At  an  exhaust  temperature  of  2000°K, 
radiative  cooling  will  produce  a rate  of  temperature  decrease  of  about  13°K  per  milli- 
second. At  an  exhaust  temperature  of  1000 °K,  the  cooling  rate  has  dropped  by  a factor 
of  six.  It  is  clear  from  this  straightforward  upper-bound  estimate  that  the  radiation 
cooling  rates  are  inconsequential  when  compared  with  the  temperature  decrease  due  to 
gas  dynamic  effects. 

The  question  of  radiative  cooling  or  heating  in  the  long-time  dispersion  wake  regime  is 
much  more  difficult  to  assess  . Considered  on  a global  basis,  a parcel  of  air  in  the 
middle  and  upper  stratosphere  (altitude  ~ 30  km)  cools  by  infrared  radiative  trans- 
port. This  cooling  results  primarily  from  energy  lost  to  space  via  emission  in  the 

CO  15— /it  band.  * In  the  lower  stratosphere,  cooling  due  to  15 -p  band  CO  emission 
Z z 

disappears  as  the  parcel  of  air  can  no  longer  "see"  the  cool  space  background. 

Rather,  at  these  lower  altitudes  a parcel  of  air  is  surrounded  by  an  atmosphere  that 
is  nearly  isothermal  over  the  mean  free  path  for  photons  emitted  in  the  15 -p  band.  A 
similar  situation  holds  for  COQ  amounts  due  to  exhaust  gases  at  the  SST  wake  altitude. 

The  situation  with  respect  to  water  vapor  in  the  SST  exhaust  is  quite  different.  There 

will  be  an  imbalance  of  H O molecules  above  the  atmospheric  ambient  level  due  to  the 

SST  wake.  In  the  long-time  dispersion  wake  the  CO  level  will  be  only  slightly  per- 

-4  1 

turbed  from  the  ambient  level  of  about  3 x 10  mass  fraction.  The  H^O  imbalance 
will  produce  two  counterposed  effects:  (1)  absorption  of  solar  radiation  in  the  2.  7-p 
band  and  (2)  emission  in  the  H^O  rotation  band  which  starts  at  approximately  20  /n. 

An  examination  of  these  two  effects  shows  that  the  far-infrared  cooling  dominates. 

A highly  accurate  determination  of  the  amount  of  cooling  would  necessitate  a detailed 


*W.  R.  Kuhn  and  J.  London,  "Infrared  Radiative  Cooling  in  the  Middle  Atmosphere 
(30-110  km)."  J.  Atmos.  Sci,  Vol.  26,  p.  189  - 204,  1969. 
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numerical  solution  of  the  equation  of  radiative  transfer  over  a number  of  wavelength 
intervals  sufficient  to  account  for  the  spectral  dependence  of  the  H^O  rotational  band. 
Rather  than  attempt  such  a laborious  calculation,  the  following  approximate  analysis 
based  on  an  optically  thin  model  was  performed.  We  assume  that  photons  emitted 
upward,  i.e. , over  a 2tt  steradian  view  factor  represented  by  emission  to  free  space, 
participate  in  the  cooling.  Then  the  cooling  rate  becomes 


(Qrad^  = 277  / B^>T)  dco 


H2° 


erg 

gm  sec 


It  is  inadequate  to  make  the  approximation  that  B(oj,  T)  ~ BWo(T)  , since  the  pure 
rotational  band  extends  over  a large  wave-number  range.  Rather,  we  approximate 
the  spectral  absorption  coefficient  variation  as 


aband  /co  \ - (cu/T)^ 

/ 7r  r \T/e 


with 


r = 


4 kT  B 

t 

he 


1/2 


-1 


and  where  B is  the  rotation  constant  which  we  take  as  30  cm  . Performing  the 
e 

above  integration  at  a wake  temperature  of  220 °K  we  arrive  at  a cooling  rate  of 


«rad> 


= 2.7  x 10 


8 


erg 


H2° 


gmH2osec 


Again,  the  energy  loss  is  shared  by  all  of  the  exhaust  species.  Hence  the  energy 
loss  per  unit  mass  of  exhaust  gases  is 
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(Qrad} 

cooling 


<QradH20  ’ Cfi2° 


where  is  the  mass  fraction  of  water  vapor.  The  mass  fraction  (i.e.  , mixing 

ratio)  of  water  vapor  in  the  ambient  atmosphere  is  about  10  6 . In  a time  on  the 

order  of  a few  minutes,  the  wake  is  dispersed  to  an  extent  that  the  wake  water  vapor 

-5 

level  is  only  one  order  of  magnitude  above  ambient  so  that  C„  _ = 10  . The 

H2° 

radiative  cooling  rates  are  sufficiently  low  that  negligible  cooling  occurs  over  the  short 
time  interval  where  C^q  exceeds  10  . At  a level  of  water  vapor  of  10_o  (equal 

to  10  times  ambient)  the  cooling  rate  yields  a rate  of  temperature  drop  of 


AT 

At 


(Qrad} 

cooling 

heat  capacity  of  air 


1 °K/hr 


Because  of  the  optically  thin  model,  the  above  rate  of  temperature  decay  is  a con- 
servative upper  bound.  From  an  examination  of  the  strengths  of  the  lines  in  HgO 
rotational  band,*  we  estimate  that  only  those  lines  which  are  located  at  wavelengths 
shorter  than  about  25  p would  emit  to  free  space.  This  would  lower  the  cooling  rate 
by  about  one  order  of  magnitude.  Furthermore,  the  mixing  of  the  wake  gases  with 
the  atmosphere  will  reduce  the  wake  water  vapor  concentration  to  further  reduce 
the  cooling  rate. 


The  rate  of  temperature  drop  does  not  persist  indefinitely.  As  soon  as  the  wake  drops 
below  ambient  temperature,  absorption  in  the  15-p  band  of  CC>2  occurs.  The  final 
temperature  reached  is  such  that  the  absorption  by  C09  balances  the  loss  of  HgO, 
i.  e.  , 


(Qrad)  = (Qrad) 

heating  a cooling 

C02  h20 


*See  Table  5.5  of  R.  M.  Goody,  Atmospheric  Radiation,  Oxford  Press  New  York,  1964. 
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The  CC>2  15  -p  band  heating  rate  is 


(Qrad),  .. 

heating 

C02 


47t  a 


band 


[B 


co 


(Tarab) 


Bco  ^Twake^ 


which  yields  a final  temperature  about  2 percent  below  ambient  based  on  the  upper- 
bound  optically  thin  model  and  if  the  water  vapor  concentration  remained  at  10  times 
ambient.  It  is  estimated  that  self-absorption  in  the  water  vapor  bands  and  wake 
mixing  will  reduce  the  final  temperature  to  about  0.2  percent  below  ambient. 
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Section  6 

EXHAUST  PARTICULATE  TRANSPORT 


The  need  to  include  the  particulate  fallout  process  as  part  of  the  wake  model  has  been 
investigated.  If  hours  are  required  for  the  particles  to  fall  an  appreciable  distance, 
say  a thousand  feet,  then  this  process  need  not  be  included.  In  the  study,  particle 
sizes  ranging  from  one  micron  in  diameter  to  a tenth  of  a micron  were  considered. 

This  range  of  particle  sizes  is  in  agreement  with  that  measured  in  the  Naval  Air 
Propulsion  Center.*  The  sinking  velocity  of  the  particles  was  determined  from  Newton's 
Second  Law  where,  in  this  case,  the  vertical  force  component  was  the  difference  be- 
tween the  gravitational  force  and  the  drag  force: 


Because  of  the  relative  size  of  the  particle  as  compared  to  the  mean  free  path  at 
typical  SST  cruise  altitudes,  a free-molecular  drag  force  was  assumed: 


where  it  can  be  shown  to  the  first  order  of  accuracy  that  the  drag  coefficient  for  a 
sphere  in  the  free  molecular  flow  regime  is-.** 


mg  - D = m (dV/dt) 


C 


*" Study  of  Altitude  and  Mach  Number  Effects  on  Exhaust  Gas  Emissions  of  an 
Afterburning  Turbofan  Engine,"  NAPTC-ATD-212,  December  1971. 


**Princeton  Series  on  High  Speed  Aerodynamics  and  Jet  Propulsion,  Volume  HI: 
Fundamentals  of  Gas  Dynamics,  p.  704. 
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In  these  expressions  p is  the  ambient  density;  V is  the  vertical  component  of  the 
particle  velocity;  d , the  particle  diameter;  R , the  gas  constant;  and  T , the 
ambient  temperature.  The  parameter  S is  defined  by 


V2  RT 

For  known  ambient  conditions,  the  drag  force  becomes 

D = V27 r RT  p (Vd2)  = (constant)  (Vd2) 

Substituting  this  expression  for  D into  Newton's  law  yields  a first  order  differential 
equation  for  the  velocity  V . Assuming  as  initial  conditions  zero  velocity  and  zero 
distance  traveled,  the  first  order  differential  equation  is  readily  integrated  yielding 

v(t)  = Tg  (1  - e t//j) 

and 


x(t)  = rg  [t  + t (e  t//‘  - 1)] 
where  the  characteristic  time  r is  defined  as 

_ m _ m 

(||- V27rRT"p)  d2  (constant)  (d2) 

For  carbon  particulate  matter  deposited  at  65,000  feet,  these  expressions  show, 
Figure  6-1,  that  2.  5 weeks  are  required  for  the  one  micron  particle  to  fall  a 
thousand  feet,  and  176  days  (approximately  six  months)  are  required  for  the  one- 
tenth  micron  particle  to  fall  the  same  distance.  From  this  it  can  be  concluded  that 
the  transport  of  particulate  matter  in  the  engine  exhaust  jet  and  in  the  vortex  field 
need  not  be  treated  separately  from  the  gas  transport  since  the  particles  will  follow 
the  gas  motions  closely. 


6-2 


6-3 


Figure  6-1  Particulate  Transport  Times 


Section  7 

INTERFACE  BETWEEN  WAKE  MODEL  AND  DIFFUSION  AND  TRANSPORT  MODEL 

In  this  section  the  fluid  mechanical  and  chemical  interface  between  the  wake  model 
and  the  diffusion  and  transport  model  under  development  by  the  Lawrence  Livermore 
Laboratory  is  discussed.  The  interface  defines  the  time  interval  after  passage  of  the 
aircraft  over  which  the  wake  model  will  predict  the  fluid  mechanical  and  chemical 
behavior  of  the  aircraft  wake.  The  interface  will  also  include  definition  of  the  output 
quantities  of  the  wake  model  which  will  be  used  as  input  quantities  to  the  diffusion  and 
transport  model.  Also  included  in  this  section  is  a discussion  of  requirements  imposed 
on  the  wake  model  by  the  basic  objectives  of  the  diffusion  and  transport  model. 

7.1  INTERFACE  DEFINITION 

The  basic  definition  of  the  wake  model  and  the  diffusion  and  transport  model  interface 
is  that  the  wake  model  must  describe  the  aircraft  wake  for  a sufficient  distance  behind 
the  aircraft  until  the  aircraft-induced  perturbations  in  the  atmosphere  become  small 
relative  to  naturally  occuring  perturbations  in  the  atmosphere.  In  other  words  the  wake 
model  must  extend  a sufficient  distance  behind  the  aircraft  until  the  aircraft-induced 
perturbations  no  longer  affect  the  fluid  mechanics  and  chemical  kinetics  of  the  wake. 

The  fluid  mechanical  interface  can  be  best  defined  by  considering  the  wake  growth  data 
and  theoretical  estimates  shown  in  Figure  2-19.  The  most  persistent  dispersion  mecha- 
nisms are  those  associated  with  the  wake  height  — namely,  wind  shear  and  gravitational 
collapse.  Wind  shear  is  dependent  on  the  height  since  the  difference  in  advection  be- 
tween the  top  and  bottom  of  the  wake  volume  is  a function  of  the  wake  height,  vertical 
wind  shear,  and  angle  between  the  wind  vector  and  the  wake  axis.  Gravitational 
collapse  is  height  dependent,  since  the  driving  potential  for  collapsing  the  wake  is  a 
function  of  the  density  difference  between  the  wake  fluid  and  the  atmospheric  air,  which 
in  turn  is  directly  dependent  on  the  wake  height. 
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The  interface  between  the  wake  regime  and  the  transport  and  dispersion  regime,  can 
then  be  defined  by  the  point  in  time  when  the  wake  height  is  no  longer  affected  by  the 
aircraft  flow  field.  This  point  is  when  the  gravitational  wake  collapse  is  completed, 
which  is  at  a wake  age  of  about  20  min  for  the  B-52  and  about  10  to  15  min  for  the 
Convair  990.  At  the  present  time,  a method  for  scaling  the  point  at  which  the  gravi- 
tational collapse  is  completed  has  not  been  developed.  It  is  hoped  that  a scaling  method 
will  be  developed  in  the  near  future.  The  fluid  mechanical  interface  can  only  be  defined 
in  terms  of  the  physical  phenomenon  at  this  time. 

The  output  quantities  from  the  wake  model  will  be  the  width  and  height  of  the  wake  and 
the  concentration  of  the  engine  exhaust  emission.  Although  the  chemical  interface 
has  not  been  firmly  established,  it  is  expected  that  the  species  concentrations  to  be 
defined  at  the  model  interface  will  include  NO,  N09,  0^,  OH,  S02>  SO^,  and, 
perhaps,  hydrocarbons.  The  solid  particulate  concentration  will  also  be  included 
which,  in  the  present  chemical  model,  is  assumed  to  be  inert. 

7.2  DIFFUSION  AND  TRANSPORT  MODEL  REQUIREMENTS  ON  WAKE  MODEL 

The  requirements  imposed  on  the  wake  model  are  dictated  by  the  impact  errors  in 
the  wake  model  outputs  will  have  on  the  diffusion  and  chemistry  in  the  long-term 
wake,  as  predicted  by  the  diffusion  and  transport  model.  For  present  purposes,  it 
is  assumed  that  the  time  frame  of  the  diffusion  and  transport  model  is  of  the  order  of 
days.  If  longer  time  periods  are  considered,  the  requirements  on  the  wake  model 
are  considerably  relaxed.  The  requirements  on  the  wake  model  are  discussed  in 
relationship  to  the  NO  /Or,  photochemistry,  wake  diffusion,  and  particulate  formation 

X o 

in  the  diffusion  and  transport  regime. 

The  half-times  for  catalytic  reduction  of  ozone  by  the  oxides  of  nitrogen  according 
to  the  analysis  discussed  in  Section  4 (Photochemistry  in  the  Wake  Regime)  are 
presented  in  Figure  7-1.  Notethatthe  half-times  are  strongly  dependent  on  the  NO 
concentrations.  The  NO  concentration  is  given  by  the  total  mass  of  NO  in  the 

X X 

wake  and  the  total  volume  over  which  it  is  distributed  (assuming  no  NO  sinks).  The 

X 

total  quantity  of  NO  in  the  wake  should  therefore  be  given  as  accurately  as  possible 

X 
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Figure  7-1  Half-Times  for  Catalytic  Reduction  of  by  NO, 


even  though  significant  NO  /O  photochemical  reactions  may  not  occur  in  the  diffusion 

X u 

and  transport  regime  as  indicated  by  Figure  7-1.  The  NO  requirement  is  probably 

X 

the  most  important  requirement  on  the  wake  model. 


The  dependence  of  long-term  wake  diffusion  in  the  diffusion  and  transport  regime  on 
the  initial  wake  size  is  illustrated  in  Figure  7-2.  These  solutions  to  the  diffusion 
equations  were  obtained  from  Ref.  7-1  where  the  diffusion  coefficient  is  taken  to  be 
of  the  form 


e1/3  CT(t,4/3 


2 3 

where  e is  the  dissipation  rate  (cm  /sec  ) and  a is  the  standard  deviation  of  the  con- 
centration. For  present  purposes , the  standard  deviation  a is  equated  to  the  wake 
width.  At  early  times,  there  is  some  effect  of  the  initial  wake  size  on  the  wake  width 
but  at  times  of  the  order  of  a day,  the  wake  width  is  weakly  dependent  on  the  initial 
wake  size.  These  results  indicate  that  the  initial  wake  size  is  not  a critical  input 
parameter  to  the  diffusion  and  transport  model. 

The  chemistry  related  to  the  formation  of  particulates  was  discussed  in  detail  in 
Section  4 (Photochemistry  in  the  Wake  Regime).  The  results  presented  in  Section  4 
indicate  that  sulphuric  acid  particulates  may  be  formed  in  the  wake  regime  due  to 
the  relatively  high  concentration  of  oxygen  atoms.  If  it  is  desirable  to  trace  the 
formation  of  particulates  in  detail,  it  will  be  necessary  to  include  the  percentage 
conversion  of  sulphur  dioxide  to  sulphuric  acid  as  a wake  model  output  quantity.  If 
only  the  long-term  formation  of  sulphuric  acid  particulates  is  of  interest,  the  relatively 
small  degree  of  conversion  in  the  wake  can  be  ignored. 

The  time  rate  of  diffusion  of  engine  effluents  into  ambient  air  as  a function  of  final  wake 
concentration  distribution  and  diffusion  coefficient  was  investigated  to  determine  the 
necessity  of  including  concentration  gradients  as  a wake  output  parameter.  The  dis- 
tributions considered  were  those  shown  in  Figure  7-3.  To  obtain  the  Gaussian,  it  has 
been  assumed  that  all  of  the  aircraft-generated  turbulent  mixing  has  died  out,  leaving  a 
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Figure  7-2  Effect  of  Initial  Wake  Size  on  Wake  Width  in  Diffusion  and  Transport  Regime 
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cylindrically  diffusing  line  trail.  This,  in  turn,  has  been  made  equivalent  to  a line 
source  of  heat  where  a certain  quantity  of  heat  is  deposited  instantaneously  per  unit 
length.  For  the  line  source  of  heat,  the  radial  distribution  of  temperature  is  given  by 


where 

c = specific  heat 

T = temperature 

Qm  = source  strength 
p = density 

k = coefficient  of  conductivity 

t = time 

The  equivalent  for  the  diffusion  of  species  is 


= original  mass  deposited  per  unit  length 
D = diffusion  coefficient 

For  convenience  in  the  computation,  the  Gaussian  was  normalized,  i.e.: 


where 


Cb  = mass  fraction  of  the  ith  species 


max 
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where 


M 


m 


max 


4 7T  p D t 


and 


]/2TDt 


The  species  mass  fractions  defining  the  height  of  the  three  step  distributions  are 
the  mean  of  the  Gaussian  as  obtained  by  integrating  the  area  under  the  curve  to  a 
radius  where  the  concentration  is,  respectively,  10,  1,  and  0.1  percent  of  the  initial 
centerline  value.  When  studying  the  relative  diffusion  rates  of  the  Gaussian  and  step, 
the  two  were  in  all  cases  terminated  at  the  same  radius. 

To  compute  the  time  rate  of  change  of  the  distributions,  the  species  conservation 
equation  for  diffusion  was  solved  numerically.  In  cylindrical  coordinates,  this 
expression  is  given  as 


For  simplicity,  a two-component  gas  mixture  was  assumed.  The  distributions 
shown  in  all  figures  are  for  nitric  oxide, whereas  its  complement  was  a species  with 
characteristics  the  same  as  that  of  air.  At  all  times,  the  sum  of  the  species  mass 
fractions  was  equal  to  one. 

If  the  time  required  for  the  Gaussian  and  step  distributions  to  converge  is  short, 
say  in  terms  of  hours  rather  than  several  days,  then  the  diffusion  and  transport 
regime  can  be  assumed  to  be  independent  of  the  final  wake  concentration  gradients. 
In  the  first  case  studied,  the  time  rates  of  decay  of  Gaussian  (terminated  at  a radius 
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where  the  concentration  was  1 percent  of  the  centerline  value)  and  step  distribution 

whose  areas  were  equal  to  that  of  the  truncated  Gaussian  were  computed  for  three 

3 4 5 2 

constant  diffusion  coefficients:  10  , 10  , and  10  cm  /sec.  Figure  7-4  shows  the 

3 2 

time  variation  for  the  ratio  of  the  centerline  mass  fractions.  For  10  cm  /sec,  little 

more  than  4 hr  were  required  for  the  centerline  values  of  the  two  distributions  to 

decay  to  within  20  percent  of  each  other,  while  30  min  and  3 min,  respectively,  were 

required  for  the  diffusion  coefficients  of  10  and  10  cm“/sec.  A diffusion  coefficient 
5 2 3 

of  10  cm  /sec  is  probably  a lower  limit  for  the  diffusion  and  transport  regime. 

The  effects  of  locating  the  "edge"  of  the  Gaussian  at  other  than  1 percent  were 

3 2 

studied  for  a diffusion  coefficient  of  10  cm  /sec.  The  results  of  the  computations  are 
shown  in  Figure  7-5.  For  the  case  where  the  centerline  values  of  the  Gaussian  and  the 
step  distribution  of  equal  area  were  initially  the  most  divergent,  a 9-hr  period  was 
required  for  the  two  to  converge  to  within  20  percent  of  each  other.  For  the  other 
two  cases,  4 hr  or  less  were  required  for  this  same  degree  of  convergence.  Thus, 
it  seems  that  wake  concentration  gradients  diffuse  rather  quickly,  within  a matter  of 
hours  for  very  low  values  of  the  diffusion  coefficient  and  within  a matter  of  minutes 
for  more  realistic  values  of  the  diffusion  coefficient,  regardless  of  the  diversity  of  the 
distributions.  Hence,  concentration  gradients  need  not  be  included  in  the  description 
of  the  wake  model  output. 

7.3  REFERENCE  FOR  SECTION  7 

7-1  J.  J.  Walton,  Scale  Dependent  Diffusion,  VCRL-74134,  Lawrence  Livermore 
Laboratory,  25  Aug  1972 
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Distribution  — Centerline  Mass  Fraction 
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Figure  7-5  Relative  Decay  of  Centerline  Mass  Fractures  of  Three  Gaussian 
and  Step  Distributions 


